Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,903,968 T→G 6.5% F75V (TTC→GTC)  nirB → nitrite reductase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,903,9680TG 6.5% 188.3 / 3.4 61F75V (TTC→GTC) nirBnitrite reductase large subunit
Reads supporting (aligned to +/- strand):  ref base T (30/27);  new base G (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGTCGAACAACTGGTCAGTCGCGGCGCCCTGCAGCGCTTCGA‑GGTGCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTGGCCACCGGCTCCTATCCCTTCGTGCCGC  >  NC_002947/1903820‑1904112
                                                                                                                                                      |                                                                                                                                                
cGTCGAACAACTGGTCAGTCGCGGCGCCCTGCAGCGCTTCGA‑GGTGCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGAtt                                                                                                                                                  <  1:45541/149‑1 (MQ=255)
       caacTGGTCAGTCGCGGCGCCCTGCAGCGCTTCGAGGGTGCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTAc                                                                                                                                             >  1:53542/1‑148 (MQ=255)
       caacTGGTCAGTCGCGGCGCCCTGCAGCGCTTCGA‑GGTGCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTAc                                                                                                                                             <  2:53542/147‑1 (MQ=255)
          cTGGTCAGTCGCGGCGCCCTGCAGCGCTTCGA‑GGTGCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGc                                                                                                                                        <  2:72320/149‑1 (MQ=255)
           tGGTCAGTCGCGGCGCCCTGCAGCGCTTCGA‑GGTGCGGGTATTCGGCGAAGAGCGCCAGACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGc                                                                                                                                        >  1:283820/1‑149 (MQ=255)
                aGTCGCGGCGCCCTGCAGCGCTTCGA‑GGTGCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATgg                                                                                                                                  >  4:70165/1‑149 (MQ=255)
                  tCGCGGCGCCCTGCAGCGCTTCGA‑GGTGCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCg                                                                                                                                <  5:180332/149‑1 (MQ=255)
                  tCGCGGCGCCCTGCAGCGCTTCGA‑GGTGCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCg                                                                                                                                >  5:67822/1‑149 (MQ=255)
                       gcgcCCTGCAGCGCTTCGA‑GGTGCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATAGCGtgcac                                                                                                                           <  1:234218/149‑1 (MQ=255)
                            cTGCAGCGCTTCGA‑GGTGCGGGTATTCGGCGAAGAGCTCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGtgcacctgca                                                                                                                      >  2:348215/1‑149 (MQ=255)
                            cTGCAGCGCTTCGA‑GGTGCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGtgcacctgca                                                                                                                      <  5:27357/149‑1 (MQ=255)
                              gCAGCGCTTCGA‑GGTGCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGtgcacctgcacc                                                                                                                    >  1:205160/1‑149 (MQ=255)
                                 gcgcTTCGA‑GGTGCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTgg                                                                                                                 <  3:70165/149‑1 (MQ=255)
                                     ttCGA‑GGTGCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGt                                                                                                               >  8:28410/1‑147 (MQ=255)
                                       cGA‑GGTGCGGGTATTCGGCGAAGAGCTCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAg                                                                                                            <  1:348215/148‑1 (MQ=255)
                                         a‑GGTGCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGcc                                                                                                         <  2:205160/149‑1 (MQ=255)
                                            gTGCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGt                                                                                                       <  6:273525/149‑1 (MQ=255)
                                              gCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGc                                                                                                     >  2:210815/1‑149 (MQ=255)
                                               cGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGc                                                                                                     <  7:28410/148‑1 (MQ=255)
                                                gggTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTg                                                                                                   >  2:20366/1‑149 (MQ=255)
                                                gggTAGTCGGCGAAGAGCGCCA‑ACTCGCGTATGACCGCGTGCACCTGTCCGAATACTGCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTg                                                                                                   <  3:275767/149‑1 (MQ=255)
                                                     ttCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTg                                                                                                   >  3:320465/1‑144 (MQ=255)
                                                       cGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATc                                                                                             <  6:67822/148‑1 (MQ=255)
                                                              gagCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACcgc                                                                                       >  2:145628/1‑147 (MQ=255)
                                                                         gcgcgTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGaa                                                                           <  2:283820/149‑1 (MQ=255)
                                                                             gTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAgtcg                                                                       >  6:116552/1‑149 (MQ=255)
                                                                               aTGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAgtcgt                                                                      >  7:319087/1‑148 (MQ=255)
                                                                               aTGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAgtcgt                                                                      >  6:133066/1‑148 (MQ=255)
                                                                                tGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATGTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCa                                                                    >  1:326341/1‑149 (MQ=255)
                                                                                tGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATGTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCa                                                                    >  1:331687/1‑149 (MQ=255)
                                                                                   ccGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTg                                                                 >  2:81237/1‑149 (MQ=255)
                                                                                              tGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGc                                                      >  3:115336/1‑149 (MQ=255)
                                                                                                tCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGAt                                                                                              >  2:123766/1‑107 (MQ=255)
                                                                                                tCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTa                                                    >  7:27016/1‑149 (MQ=255)
                                                                                                tCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTa                                                    <  5:116552/149‑1 (MQ=255)
                                                                                                  cGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGcgccgc                                                                              <  3:208547/121‑1 (MQ=255)
                                                                                                  cGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGcgccgc                                                                              >  4:208547/1‑121 (MQ=255)
                                                                                                   gAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATgg                                                 <  1:20366/149‑1 (MQ=255)
                                                                                                   gAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATgg                                                 >  5:82078/1‑149 (MQ=255)
                                                                                                   gAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGAGTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATgg                                                 <  3:85532/149‑1 (MQ=255)
                                                                                                                gcAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTg                                    <  1:81237/149‑1 (MQ=255)
                                                                                                                       tgCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTg                              <  1:210815/148‑1 (MQ=255)
                                                                                                                          cGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAg                                                                          >  4:214775/1‑101 (MQ=255)
                                                                                                                          cGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAg                                                                          >  4:219109/1‑101 (MQ=255)
                                                                                                                          cGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAg                                                                          <  3:219109/101‑1 (MQ=255)
                                                                                                                          cGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAg                                                                          <  3:214775/101‑1 (MQ=255)
                                                                                                                          cGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTGGCCa                          >  4:150482/1‑149 (MQ=255)
                                                                                                                             aaCCCTGGCCCTGTTTGAGCACGAATTCTCCGATGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTGGCCACCg                       <  8:319087/149‑1 (MQ=255)
                                                                                                                              aCCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTGGCCACCgg                      >  3:230247/1‑149 (MQ=255)
                                                                                                                                ccTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTGGCCACCGGCt                    >  4:319199/1‑149 (MQ=255)
                                                                                                                                 cTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTgg                                                                                                  >  1:64562/1‑70 (MQ=255)
                                                                                                                                 cTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTgg                                                                                                  <  2:64562/70‑1 (MQ=255)
                                                                                                                                    gCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTGGCCACCGGCTCCTa                >  8:12246/1‑149 (MQ=255)
                                                                                                                                     cccTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTGGCCACCGGCTCCTAt               <  4:230247/149‑1 (MQ=255)
                                                                                                                                     cccTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTGGCCACCGGCTCCTAt               <  3:319199/149‑1 (MQ=255)
                                                                                                                                      ccTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTGGCCACCGGCTCCTATc              >  4:297900/1‑149 (MQ=255)
                                                                                                                                               gCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTGGCCACCGGCTCCTATCCCTTCGt       <  6:127098/147‑1 (MQ=255)
                                                                                                                                               gCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTGGCCACCGGCTCCTATCCCTTCGt       >  5:127098/1‑147 (MQ=255)
                                                                                                                                               gCACGATGTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTGGCCACCGGCTCCTATCCCTTCGTgc     <  2:331687/149‑1 (MQ=255)
                                                                                                                                               gCACGATGTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTGGCCACCGGCTCCTATCCCTTCGTgc     <  2:326341/149‑1 (MQ=255)
                                                                                                                                                cACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTGGCCACCGGCTCCTATCCCTTCGTgcc    >  2:18064/1‑149 (MQ=255)
                                                                                                                                                  cGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTGGCCACCGGCTCCTATCCCTTCGTgccgc  <  4:320465/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
CGTCGAACAACTGGTCAGTCGCGGCGCCCTGCAGCGCTTCGA‑GGTGCGGGTATTCGGCGAAGAGCGCCA‑ACGCGCGTATGACCGCGTGCACCTGTCCGAATACTTCGGCGGCAGCTGTGCCGAAACCCTGGCCCTGTGTGAGCACGATTTCTACGAGGCCAATGGCGTGCACCTGCACCTGGGTGAGGCCGTGCTGGAGATCGACCGCGAGCGCCGCGAAGTCGTCACTGCCGCAGGTCGCTATGGCTACGACCAACTGGTGCTGGCCACCGGCTCCTATCCCTTCGTGCCGC  >  NC_002947/1903820‑1904112

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑