Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,541,077 G→C 8.0% A505G (GCC→GGC)  narZ ← nitrate reductase 2 (NRZ), alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,541,0770GC 8.0% 299.3 / 15.2 87A505G (GCC→GGC) narZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (41/39);  new base C (4/3);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01

GTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGG  >  NC_000913/1540932‑1541225
                                                                                                                                                 |                                                                                                                                                      
gTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTg                                                                                                                                                     >  4:661644/1‑149 (MQ=255)
  tCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTggg                                                                                                                                                   <  4:177074/149‑1 (MQ=255)
    aTGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGTGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTgggg                                                                                                                                                  >  1:790479/1‑148 (MQ=255)
     tGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGgtg                                                                                                                                                <  3:862493/149‑1 (MQ=255)
     tGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGgt                                                                                                                                                 >  2:750177/1‑148 (MQ=255)
      gTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGgtgt                                                                                                                                               >  4:206809/1‑149 (MQ=255)
      gTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtg                                                                                                                                                <  2:371756/148‑1 (MQ=255)
       tCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGgtgt                                                                                                                                               >  2:535418/1‑148 (MQ=255)
       tCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtga                                                                                                                                               <  2:272301/148‑2 (MQ=255)
          aTGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACgg                                                                                                                                           >  4:800744/1‑149 (MQ=255)
            gtgATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGtt                                                                                                                                         <  3:847322/149‑1 (MQ=255)
             tgATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGttt                                                                                                                                        >  3:228478/1‑149 (MQ=255)
               aTACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTg                                                                                                                                       <  1:13701/148‑1 (MQ=255)
               aTACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTg                                                                                                                                       >  2:409693/1‑148 (MQ=255)
               aTACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTg                                                                                                                                       <  2:289548/148‑1 (MQ=255)
                tACCAGCGGTTAACACCTGCGCCGAGGATAATCATCCAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGAt                                                                                                                                     <  3:800744/149‑1 (MQ=255)
                   cAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTc                                                                                                                                  >  3:658627/1‑149 (MQ=255)
                   cAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTc                                                                                                                                  >  4:76574/1‑149 (MQ=255)
                   cAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTgg                                                                                                                                                    >  1:189379/1‑131 (MQ=255)
                   cAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTgg                                                                                                                                                    <  2:189379/131‑1 (MQ=255)
                     gTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAg                                                                                                                                <  4:73383/149‑1 (MQ=255)
                         ttAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTAc                                                                                                                                             <  1:83321/132‑1 (MQ=255)
                         ttAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTAc                                                                                                                                             >  2:83321/1‑132 (MQ=255)
                          tAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTaa                                                                                                                          <  3:206809/149‑1 (MQ=255)
                          tAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTa                                                                                                                           >  1:874736/1‑148 (MQ=255)
                               ccTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCttt                                                                                                                     <  3:76574/149‑1 (MQ=255)
                               ccTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCtt                                                                                                                      >  1:386333/1‑148 (MQ=255)
                               ccTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCtt                                                                                                                      <  1:750177/148‑1 (MQ=255)
                                cTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCttt                                                                                                                     <  2:845015/148‑1 (MQ=255)
                                          gATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGtttt                                                                                                          >  4:332048/1‑149 (MQ=255)
                                           aTAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGtttt                                                                                                          <  2:386333/148‑1 (MQ=255)
                                           aTAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑c                                                                                                        <  3:831081/149‑1 (MQ=255)
                                             aaTCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CAt                                                                                                      <  3:684141/149‑1 (MQ=255)
                                               tcatcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACACCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATc                                                                                                     <  2:783686/148‑1 (MQ=255)
                                                 atcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCTAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCtt                                                                                                   <  1:596656/148‑1 (MQ=255)
                                                 atcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCtt                                                                                                   <  1:409693/148‑1 (MQ=255)
                                                 atcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCtt                                                                                                   <  1:549044/148‑1 (MQ=255)
                                                   cGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTcc                                                                                                 >  1:757474/1‑148 (MQ=255)
                                                       cgcCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGccc                                                                                            <  4:90401/149‑1 (MQ=255)
                                                            aTGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGa                                                                                        <  2:431503/148‑1 (MQ=255)
                                                            aTGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGAt                                                                                       >  3:139012/1‑149 (MQ=255)
                                                             tGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGAt                                                                                       <  2:757474/148‑1 (MQ=255)
                                                                       gggCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAg                                                                                                                                >  4:510044/1‑99 (MQ=255)
                                                                       gggCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAg                                                                                                                                <  3:510044/99‑1 (MQ=255)
                                                                             tATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGcc                                                                      <  4:228478/149‑1 (MQ=255)
                                                                              aTCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCa                                                                     <  4:139012/149‑1 (MQ=255)
                                                                                 ggCAAATTCACGGGCGATGGTTTCAATATACTGGTGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTa                                                                  >  4:301004/1‑149 (MQ=255)
                                                                                 ggCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGt                                                                                                                                          <  1:123572/79‑1 (MQ=255)
                                                                                 ggCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGt                                                                                                                                          >  2:123572/1‑79 (MQ=255)
                                                                                  gCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTAc                                                                 >  4:713404/1‑149 (MQ=255)
                                                                                   cAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATGCGCCAGTAc                                                                 >  1:7250/1‑148 (MQ=255)
                                                                                      aTTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCaa                                                              <  2:874736/148‑1 (MQ=255)
                                                                                          aCGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCAt                                                         >  3:196505/1‑149 (MQ=255)
                                                                                             ggCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTc                                                                                                                                  <  1:63615/75‑1 (MQ=255)
                                                                                             ggCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTc                                                                                                                                  >  2:63615/1‑75 (MQ=255)
                                                                                             ggCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATaaa                                                      <  4:649495/149‑1 (MQ=255)
                                                                                             ggCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATaaa                                                      >  4:646429/1‑149 (MQ=255)
                                                                                                gATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGc                                                   >  3:24083/1‑149 (MQ=255)
                                                                                                   ggTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGa                                                <  4:481445/149‑1 (MQ=255)
                                                                                                    gTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGac                                               <  3:307782/149‑1 (MQ=255)
                                                                                                          atatACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCgg                                         >  3:148119/1‑149 (MQ=255)
                                                                                                           tataCTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCgg                                         >  2:353426/1‑148 (MQ=255)
                                                                                                           tataCTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGa                                        >  3:143541/1‑149 (MQ=255)
                                                                                                             taCTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACa                                      <  3:431890/149‑1 (MQ=255)
                                                                                                                tGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAgg                                    >  2:55656/1‑148 (MQ=255)
                                                                                                                tGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAgg                                    >  1:216893/1‑148 (MQ=255)
                                                                                                                 ggCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATaa                                                       >  1:367464/1‑128 (MQ=255)
                                                                                                                 ggCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATaa                                                       <  2:367464/128‑1 (MQ=255)
                                                                                                                  gcgcGGCCCGCC‑GTAATTTGCTCCCCCCCGGCTGGGGGGTAAGGTTTTGTTTCAGCCGTAA‑CTTTCGGCATGTTTTCCATC‑TCCAGCCCGGCAATGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGtt                                <  4:713376/149‑1 (MQ=255)
                                                                                                                  gcgcGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGgtgt                                 <  4:658627/149‑1 (MQ=255)
                                                                                                                  gcgcGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGAACACCGGACAGgtg                                  <  2:790479/148‑1 (MQ=255)
                                                                                                                          cGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATc                          >  2:767201/1‑148 (MQ=255)
                                                                                                                           gCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGa                        <  3:661644/149‑1 (MQ=255)
                                                                                                                                tAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTaa                    <  1:487209/148‑1 (MQ=255)
                                                                                                                                  aaccgtacACCCCAGCCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGt                 >  4:906794/8‑149 (MQ=255)
                                                                                                                                   tttGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCTAGGCCGTAATTCGCCAGTACCAAATCATAATCGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTc                >  4:672775/1‑149 (MQ=255)
                                                                                                                                      gtacACCCCAGCCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAg              >  1:624977/4‑148 (MQ=255)
                                                                                                                                       tacACCCCAGCCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGt            >  4:253327/3‑149 (MQ=255)
                                                                                                                                        tCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGtt           >  4:713551/1‑149 (MQ=255)
                                                                                                                                        tCACCCCAGGCTAGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGtt           >  4:567001/1‑149 (MQ=255)
                                                                                                                                         cACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGtt           >  2:282092/1‑148 (MQ=255)
                                                                                                                                         cACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGtt           >  2:806825/1‑148 (MQ=255)
                                                                                                                                          aCCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGtttt         >  4:643805/1‑149 (MQ=255)
                                                                                                                                           ccccAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGttttt        >  4:720911/1‑149 (MQ=255)
                                                                                                                                           ccccAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGtttt         >  1:97983/1‑148 (MQ=255)
                                                                                                                                             ccAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTaa      <  3:713404/149‑1 (MQ=255)
                                                                                                                                             ccAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTa       >  2:10537/1‑148 (MQ=255)
                                                                                                                                                ggCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACg    <  2:216893/148‑1 (MQ=255)
                                                                                                                                                 gCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACggg  <  3:301004/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
GTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAG‑CGTAATCTTTCGCACTGTTTT‑CATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGG  >  NC_000913/1540932‑1541225

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑