Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,549,5950CA21.7% 119.9 / 13.9 69T732N (ACC→AAC) fdnGformate dehydrogenase‑N, alpha subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base C (30/24);  new base A (2/13);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 7.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCCGCACGAACCGCAATCTGACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTGCTGTACAACC  >  NC_000913/1549450‑1549741
                                                                                                                                                 |                                                                                                                                                  
cagcCGCACGAACCGCAAGCTGACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGc                                                                                                                                                 <  1:632778/149‑1 (MQ=255)
     gCACGAACCGCAATCTGACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCAtctt                                                                                                                                            >  2:191629/1‑149 (MQ=255)
     gCACGAACCGCAATCTGACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTGATTGCGACGAAAGGTCAGTAGCTGAGTAGCTATGCGCATCTGCGTGATGACGGTCCAACCGCAtctt                                                                                                                                            >  1:284691/1‑149 (MQ=255)
      cACGAACCGCAATCTGACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCAtcttc                                                                                                                                           >  2:926757/1‑149 (MQ=255)
        cGAACCGCAATCTGACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGGCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATcttctt                                                                                                                                         <  1:123789/149‑1 (MQ=255)
          aaCCGCAATCTGACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGc                                                                                                                                       <  1:104034/149‑1 (MQ=255)
          aaCCGCAATCCGACGAAGTGGCTAAAGAGCACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGc                                                                                                                                       <  1:108244/149‑1 (MQ=255)
           aCCGCAATCTGACGCAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTAGGACGCTAATGGCGTGCGTATTGCGAAGAAAGGGCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTATTGCt                                                                                                                                      <  1:957161/149‑1 (MQ=255)
           aCCGCAATCAGACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCTAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCt                                                                                                                                      <  1:903883/149‑1 (MQ=255)
             cGCAATCTGACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCGGCGTGATGCCGGTACAACCGCATCTTCTTGCTgg                                                                                                                                    >  1:523120/1‑149 (MQ=255)
              gCAATCTCACGAAGTGTCTAAAGAGAACAACGTCTACGGGCTGGAAGATCTCTATGTAGCTAATGTCGTGCTTATTGCGAATAAAGTTCAGTTTATGAGTATCTTTGCTCATCTGCGTGATGACGGTACAAACTCATCTTATTGCTGGa                                                                                                                                   <  2:154633/149‑1 (MQ=255)
                   cTGACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGCTGACGGTACAACCGCATCTTCTTGCTGGATCTac                                                                                                                              >  2:972732/1‑149 (MQ=255)
                     gACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTacac                                                                                                                            >  2:280482/1‑149 (MQ=255)
                       cGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTCCGAAGAAACGTCAGTTGATCAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTAATGCTGGATCTACACCg                                                                                                                          >  2:424603/1‑149 (MQ=255)
                         aaGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGt                                                                                                                        <  1:373638/149‑1 (MQ=255)
                               ctcacgcgaACACCGGCTCCGCGCTGGAAGAGCTCTATGACGCTAAGGCCGTGCTTCTTGCGAAGAACGGCCAGTTGCCGAGTAGCTTTGCGCATCTGCGTTATTTTTTTACATATTCTGCGACTAGCTTGATCTACACCGGTAGCTgg                                                                                                                  <  2:75402/142‑1 (MQ=255)
                                tAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGTTCTACACCGGTAGCTGGa                                                                                                                 >  2:850658/1‑149 (MQ=255)
                                 aaaGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTTATGTTTTTAATTATTCTTCTTCTTGCTTGATCTACACCGGTAGCTGGAc                                                                                                                <  2:654023/149‑1 (MQ=255)
                                 aaaGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGAc                                                                                                                >  2:321167/1‑149 (MQ=255)
                                    gagaACAACGGCTAGGCGCTGGAAGATCTCTATGTCGCTAATGGCGTGCTTATTGCGAAGAATGTTCAGTTGCTGCGTAGCTTTGCGCATCTGCGTGATGACGGTAGAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACaga                                                                                                             <  1:126984/149‑1 (MQ=255)
                                    gagaACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTAGCGCATCTGCGTGATGACGGTACAAACGCATCTTATTGCTGGATCTACACCGGTAGATGGACaga                                                                                                             >  1:750902/1‑149 (MQ=255)
                                     agaACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACagag                                                                                                            >  1:1281142/1‑149 (MQ=255)
                                       aacaacGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGcc                                                                                                          >  1:678699/1‑148 (MQ=255)
                                        acaacGGCTACGCGCAGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGAGTTATTTTTTTAATTATTCTTCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAg                                                                                                         <  2:367736/149‑1 (MQ=255)
                                         ccacGGCTACGCGCTGGAAGATCGCTTTGTCGCTAATGGCGGGCTGACTGAGAAGAAAGGTCAGGTGCTGAGTAGCTTTGCGCATCTGCGTTTTTTTTTTAATTATTCTTCTTCTTGCGTGATCTACAGCGGTAGCTGGACAGAGCAgg                                                                                                        <  2:998796/147‑1 (MQ=255)
                                         caacGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGATTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTACTTGCTGGATCTACACCCGTAGCTGGACATAGCAgg                                                                                                        >  2:309229/1‑149 (MQ=255)
                                          aacGGCTCCGCGCTGGAAGATCTCTATGACGCAAATGGCGCGCTTAATGCGAAGAAAGGTCAGTGGCAGAGTAGCTTTGCGCATCTGCGTTTTTTTTTTAATTATTCTTCTTCTTGCGTGATCTACACCGGTAGCTGGACAGAGCAggg                                                                                                       <  2:864429/149‑1 (MQ=255)
                                               cTACGCGCTGGAAGCTCTCTATGCCGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAAcc                                                                                                  <  2:1045119/149‑1 (MQ=255)
                                                            gATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGCGCAGGGCAACCAGATGGCTAACcg                                                                                     >  1:1230903/1‑149 (MQ=255)
                                                              tctcTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACcgcg                                                                                   >  2:814271/1‑149 (MQ=255)
                                                               gtctATGACTCTAATGGCGGGCTTTTTGCGAAGAATGTTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATTTCTGTACAACTGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGa                                                                                  <  1:716718/148‑1 (MQ=255)
                                                               ctctATGACGCTAATGGTGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGa                                                                                  <  1:690083/149‑1 (MQ=255)
                                                               ctctATGACGCTAATGGCGTGCTTAATGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGTAGTTAATTAAGCATCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGa                                                                                  <  2:382013/149‑1 (MQ=255)
                                                                   aTGACCCTAATGGCGTGCTCATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGGATCTGATTTTTTTTTTTAATTATTCTTCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAAc                                                                              <  2:290174/149‑1 (MQ=255)
                                                                    tGACGCTAATGGCGTGCTTATTGAGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACAACGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACt                                                                             >  2:658832/1‑149 (MQ=255)
                                                                     gACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTg                                                                                                                   <  2:1031851/110‑1 (MQ=255)
                                                                     gACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTg                                                                                                                   >  1:1031851/1‑110 (MQ=255)
                                                                      ccGCTAATGGCGTGCTTATTGCGAAGAAAGGGCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGATTTTAATTATTCTTCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTcc                                                                           <  2:354141/148‑1 (MQ=255)
                                                                        gCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTTTTTACTGTAAAAAATCTTCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGa                                                                         <  2:1152810/149‑1 (MQ=255)
                                                                        gCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGa                                                                         >  1:1068289/1‑149 (MQ=255)
                                                                        gCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGa                                                                         <  1:926757/149‑1 (MQ=255)
                                                                         cTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGAc                                                                        >  2:1162984/1‑149 (MQ=255)
                                                                           aaTGCCGTGCTTATAGCGAAGAAAGGTCAGTTGCTGAGGAGCTTTGCGCATCTGCTTTATTTATTTAATTATTCTGCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGAccc                                                                      <  2:824965/149‑1 (MQ=255)
                                                                              ggCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTc                                                                   <  1:500722/149‑1 (MQ=255)
                                                                                      ttttGCGAAGAAAGGTCAGTTGCGGAGTATCTTTGCGGATCTGCGTGTTGACGGTACAACCGCATCTGCTTGGTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTgg                                                           <  1:1048253/147‑1 (MQ=255)
                                                                                       tttGCGAAGAAAGGTCAGTTGCTGAGTAGGTTTGCGCATCTGCGTGAAGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCACATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTggg                                                          <  1:544083/148‑1 (MQ=255)
                                                                                             atgaaAGGTCAGTTGCTGACTAGCTTTGCGCATCTGCGTGATGCCGGTACAACCGCAGCTTCGTGCTGCATTTCCCCCGGTAGCTGGACAGATCAGGGCAACCACATGGCTAACCGCGAGAACTCCGACAACTCCCGTCTCGGGAATAc                                                    >  2:405569/3‑149 (MQ=255)
                                                                                               gaaAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGc                                                  <  1:280482/149‑1 (MQ=255)
                                                                                                  aGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTACCCGCGATAACTCCGACCCGTCCGGTCTGGGGAATTCGCTgg                                               >  1:1161168/1‑149 (MQ=255)
                                                                                                  aGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTgg                                               >  2:213202/1‑149 (MQ=255)
                                                                                                    gTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAAACGCATCTTCTTGCTGGATCTACACCGGTAGCTTGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACTCTgaga                                             >  1:301775/1‑149 (MQ=255)
                                                                                                     tCAGTTGCTGAGTAGCTTTGCGCGTCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGAt                                            >  2:203802/1‑149 (MQ=255)
                                                                                                     tCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACTACCTCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGAt                                            <  2:750902/149‑1 (MQ=255)
                                                                                                     tCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGAt                                            >  1:939964/1‑149 (MQ=255)
                                                                                                       aGTTGCTGAGTAGCGGTGCGCATCTGCGTGTTTATGTTAATTATTCTTCTTCTTGCGTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATgg                                          <  2:203064/149‑1 (MQ=255)
                                                                                                         ttGCTGAGTAGCTTTGCGCATCTGCGTGATGAGGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGc                                        <  1:119089/149‑1 (MQ=255)
                                                                                                                   gcagtGCGCAGCTGCGGGATGATGTTACATCAGCATCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGTGGAATACGCTGGGATTGGCCTGGGCGTgg                              <  2:707596/145‑1 (MQ=255)
                                                                                                                    cTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCGTGCTGGATCTACAGCGGTAGCTGGACAGAGCAGGGGAACCAGATGGCTAACCGCGATAAATCCGACCCGTCCGGGGGGGGGTTTACGGTGGGCTGGGCCTGGGCGTgga                             >  2:731501/1‑148 (MQ=255)
                                                                                                                         cgcATCTGCGTGATGATGTTACTTCCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTc                        <  2:523120/149‑1 (MQ=255)
                                                                                                                            aTCTGCGTGATGTATGTATAAATGCTTCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAAc                     <  2:933690/149‑1 (MQ=255)
                                                                                                                            aTCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTTACCGCGATAACTCCGCCCCGTCAGGTCTGTGGAATAAGCTGGGATGGGCCTGGGAGTGGCAGcaaaaa                     >  2:225189/1‑144 (MQ=255)
                                                                                                                                gCGTGATGACGGGACAACCGCATCTTCTTGCTGGATCAACCCCGGTAGCGGGACAGAGCACGGCAACCATATGGCTAACCGCGATAACTCCGACCCCGCCGGGCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTc                 <  1:1018042/149‑1 (MQ=255)
                                                                                                                                    gatgaCGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTACGACCCGTCCCGTCTGGGGAATACGCTGGGATGGGCCTCGGCGTGGCCGCTCAACCGTCGCGt             >  1:628646/1‑149 (MQ=255)
                                                                                                                                     atgaCGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTg            >  2:757149/1‑149 (MQ=255)
                                                                                                                                     atgaCGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTg            >  1:1064004/1‑149 (MQ=255)
                                                                                                                                       gaCGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTGCt          >  1:936509/1‑149 (MQ=255)
                                                                                                                                        aCGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTGCTg         <  1:321167/149‑1 (MQ=255)
                                                                                                                                             accaCCGCATCTTCCTGCTGGACCTCCCCCGGTAGCTGGCCAGAGCAGGGCAACCAGAGGGCTAACCGCGCTAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTGCTGTACaa    <  1:154157/146‑1 (MQ=255)
                                                                                                                                               aaCCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACACGTCCGGTCTGGGGAATACGCTGGGAAGGGCCTGGGAGTGGCCGCTCAACCGTCGCGTGCTGTACAAcc  >  1:959056/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CAGCCGCACGAACCGCAATCTGACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTGCTGTACAACC  >  NC_000913/1549450‑1549741

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