Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,951,3920GC29.4% 139.6 / 69.1 85intergenic (‑128/+33)oguA/PP_25858‑oxoguanine deaminase/alpha‑ketoglutaric semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (27/33);  new base C (14/11);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGCACAC  >  NC_002947/2951244‑2951535
                                                                                                                                                    |                                                                                                                                                
tGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGa                                                                                                                                                 <  2:1290590/150‑1 (MQ=255)
tGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGa                                                                                                                                                 <  2:1290629/150‑1 (MQ=255)
 gCAGGGTTTGTAGGTGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGaa                                                                                                                                                <  2:2102006/150‑1 (MQ=255)
 gCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGaa                                                                                                                                                <  1:1050938/150‑1 (MQ=255)
 gCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGaa                                                                                                                                                <  2:1737200/150‑1 (MQ=255)
    gggTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAagag                                                                                                                                             <  2:469905/150‑1 (MQ=39)
     ggTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCaa                                                                                                                                                <  1:887815/146‑1 (MQ=255)
     ggTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCATCTTCGCGGCTCTAGCCGCCaa                                                                                                                                                >  2:887815/1‑146 (MQ=255)
      gTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAAGAGGc                                                                                                                                           <  2:290970/150‑1 (MQ=39)
       tgcgTAGATGACATGGTTCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTTGCGGCTAGAGCCGCGAAGAGGcc                                                                                                                                          >  2:1995194/4‑150 (MQ=255)
              aTGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAAGAGGCCGGCACgg                                                                                                                                   <  2:544550/150‑1 (MQ=38)
               tGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACggg                                                                                                                                  <  2:679117/150‑1 (MQ=255)
                      gTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGc                                                                                                                           >  2:639292/1‑150 (MQ=255)
                        gCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCa                                                                                                                         <  1:612218/150‑1 (MQ=255)
                        gCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGGATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCa                                                                                                                         >  1:1008729/1‑150 (MQ=255)
                           aGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAagag                                                                                                                                             >  1:829964/1‑127 (MQ=38)
                           aGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAagag                                                                                                                                             <  2:829964/127‑1 (MQ=38)
                                      tACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACAcc                                                                                                           <  1:1632920/150‑1 (MQ=255)
                                         gCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCATACCGcc                                                                                                        >  1:1066943/1‑150 (MQ=255)
                                                cGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGgc                                                                                                <  1:1725377/150‑1 (MQ=255)
                                                 gATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATGCAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCaa                                                                                                                                                <  1:784439/102‑1 (MQ=37)
                                                 gATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGc                                                                                                                                           >  1:924622/1‑107 (MQ=255)
                                                 gATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGc                                                                                                                                           <  2:924622/107‑1 (MQ=255)
                                                      ccGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTca                                                                                          >  1:160394/1‑150 (MQ=255)
                                                              aTAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTcc                                                                                  >  1:936604/1‑150 (MQ=255)
                                                              aTAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAaga                                                                                                                                              <  2:600627/91‑1 (MQ=25)
                                                              aTAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAaga                                                                                                                                              >  1:600627/1‑91 (MQ=25)
                                                               tAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCa                                                                                 >  1:2113006/1‑150 (MQ=255)
                                                                aGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAg                                                                                <  1:2037896/150‑1 (MQ=255)
                                                                  gcgcTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCgg                                                                                                 >  1:2085225/1‑131 (MQ=255)
                                                                  gcgcTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCgg                                                                                                 <  2:2085225/131‑1 (MQ=255)
                                                                    gcTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAAc                                                                            <  2:2113006/150‑1 (MQ=255)
                                                                       cGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCgg                                                                         <  1:639292/150‑1 (MQ=255)
                                                                         gTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCgg                                                                                                 >  1:1130910/1‑124 (MQ=255)
                                                                         gTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCgg                                                                                                 <  2:1130910/124‑1 (MQ=255)
                                                                         gTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCaa                                                                                                                                                >  1:444465/1‑78 (MQ=25)
                                                                         gTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCaa                                                                                                                                                <  2:444465/78‑1 (MQ=25)
                                                                              gTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATcc                                                                  >  2:2143791/1‑150 (MQ=255)
                                                                               tAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCAGGCTGGCGTGGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCAAACCATCCAGCAACCGGGCAATccc                                                                 >  1:664241/1‑150 (MQ=255)
                                                                                  gcaagcTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAg                                                              >  2:849427/1‑150 (MQ=255)
                                                                                               tGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTca                                                                                          >  1:1028126/1‑109 (MQ=255)
                                                                                               tGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTca                                                                                          <  2:1028126/109‑1 (MQ=255)
                                                                                               tGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGttt                                                 >  1:1632461/1‑150 (MQ=255)
                                                                                               tGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGttt                                                 <  2:1893025/150‑1 (MQ=255)
                                                                                               tGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGttt                                                 >  2:98295/1‑150 (MQ=255)
                                                                                               tGAATGTATACAGCTCGAAGTGCCCGTGCCAGCCTCGTGTCGGCTAGAGCTGCGAAGAAGCCGGCACGGGTTGCCGGCATGAAAGAAACCCCAGCCTGCTCCGCAAGAACACGGTAGAGCGACCGGGCAGTCCCCAGCGGGTTAGCGttt                                                 >  2:984305/1‑150 (MQ=255)
                                                                                                  ttGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTc                                              >  2:1868237/1‑150 (MQ=255)
                                                                                                    gTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGTGATAATAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCaa                                            <  1:1165163/150‑1 (MQ=255)
                                                                                                     tataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCATACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAAc                                           <  2:1066943/150‑1 (MQ=255)
                                                                                                     tataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGGTGCTCACACCCTCGAGCAACCTGGGAAGCCCCAGGGGGTTAGCGCTTTTCAAc                                           >  2:1165163/1‑150 (MQ=255)
                                                                                                     tataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCaa                                            >  2:1165203/1‑149 (MQ=255)
                                                                                                     tataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGAGGCCAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCaa                                            <  1:1165203/149‑1 (MQ=255)
                                                                                                         cAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCAt                                       <  1:2143791/150‑1 (MQ=255)
                                                                                                                tgcaacccGTGCCGGCCTCTTCGCGGCTCTAGCCCACAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGgcag                                >  1:1795445/6‑150 (MQ=37)
                                                                                                                  gTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGgcagca                              <  1:1868237/150‑1 (MQ=255)
                                                                                                                  gTGCCCGTGCCGGCCTCTTGGAGATTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCCCACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGgcagca                              <  2:1008729/150‑1 (MQ=255)
                                                                                                                        gTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGt                        <  2:1221117/150‑1 (MQ=255)
                                                                                                                         tGCCGGCCTCTTGGCGGCGAGAGCCGCGAAGAGTCCGGCACGTGTTTCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTc                       <  2:2079283/150‑1 (MQ=255)
                                                                                                                            cGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGt                            >  1:1587917/1‑142 (MQ=255)
                                                                                                                            cGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGt                            <  2:1587917/142‑1 (MQ=255)
                                                                                                                            cGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAgg                    <  2:664241/150‑1 (MQ=255)
                                                                                                                             ggCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAggg                   <  2:936604/150‑1 (MQ=255)
                                                                                                                             ggCCTCTTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAggg                   >  2:947367/1‑150 (MQ=39)
                                                                                                                             ggCCTCTTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAggg                   >  2:1263555/1‑150 (MQ=39)
                                                                                                                              gCCTCTTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGt                  >  1:388824/1‑150 (MQ=39)
                                                                                                                                ctctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGAAg                >  1:941152/1‑150 (MQ=255)
                                                                                                                                ctctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAg                >  2:997163/1‑150 (MQ=255)
                                                                                                                                ctctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAg                >  1:1083389/1‑150 (MQ=39)
                                                                                                                                ctctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAg                >  2:2160517/1‑150 (MQ=39)
                                                                                                                                 tctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGt               <  2:1583245/150‑1 (MQ=255)
                                                                                                                                 tctTCGCGGCTCTAGCCGCCAAGAGGCCTGCACGGGTTGCAGCCATGAAAGATCACACCGCCTcgc                                                                                                   >  1:2148334/1‑66 (MQ=17)
                                                                                                                                 tctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTcgc                                                                                                   <  2:2148334/66‑1 (MQ=17)
                                                                                                                                 tctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGt               >  2:59454/1‑150 (MQ=39)
                                                                                                                                 tctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGt               >  2:2141488/1‑150 (MQ=39)
                                                                                                                                   ttGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTc             <  2:1632461/150‑1 (MQ=255)
                                                                                                                                    tGGCGGCGAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCt            <  1:98295/150‑1 (MQ=255)
                                                                                                                                      gcggcgAGAGCCGTGAAGAGGCCGGTACGGGTTGCAGCCATGAAAGATCATTCCGCTGCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTTAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTgg          <  2:957941/150‑1 (MQ=255)
                                                                                                                                      gcggcTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGTTCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCTGTGTTCAACTCATCCGGCAGCAGTTCGTCAGGGTAGTTCTgg          >  2:1518868/1‑150 (MQ=255)
                                                                                                                                       cggcTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGa         >  2:935100/1‑150 (MQ=255)
                                                                                                                                        ggctctagCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGaa        <  1:2160517/144‑1 (MQ=255)
                                                                                                                                        ggctctagCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGaa        <  2:388824/144‑1 (MQ=255)
                                                                                                                                          cTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGc      >  1:1385190/1‑150 (MQ=255)
                                                                                                                                           tAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAg       >  2:1898104/1‑148 (MQ=255)
                                                                                                                                           tAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAg       <  1:1898104/148‑1 (MQ=255)
                                                                                                                                              agCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGCacac  >  2:1879512/1‑150 (MQ=255)
                                                                                                                                                  gcgAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAgg                    >  1:1792968/1‑128 (MQ=255)
                                                                                                                                                  gcgAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAgg                    <  2:1792968/128‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGCACAC  >  NC_002947/2951244‑2951535

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: