Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,537,123 G→C 5.1% intergenic (+52/+51) nadE → / ← PP_4870 NH(3)‑dependent NAD(+) synthetase/azurin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,537,1230GC 5.1% 342.8 / 5.9 109intergenic (+52/+51)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base G (54/48);  new base C (2/4);  total (56/53)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

CATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGAA  >  NC_002947/5536979‑5537272
                                                                                                                                                |                                                                                                                                                      
cATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCACGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGACCCGCGAATACGAGGGTGAACCCACCGACGCa                                                                                                                                                   <  1:1218539‑M2/150‑2 (MQ=255)
       aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCATCTATCGAAAGACCCGGCACAAGCGCGAACTGCCGTTCGCGCTTTGATAAGTGGCTTGCCGGAGCGGGTTCACCCGCGAATACGAGAGTGACCCCACCGACGCattcgcgg                                                                                                                                            >  2:345625‑M2/1‑142 (MQ=255)
       aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCattcgcgg                                                                                                                                            >  2:88958‑M2/1‑142 (MQ=255)
           ggTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCattcgcgggtga                                                                                                                                        <  2:1292669‑M2/150‑13 (MQ=255)
            gTTGCCCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCattggcgggtgaa                                                                                                                                       >  2:1052147‑M2/1‑137 (MQ=255)
            gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTCCCGTTCGCTCATTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCattcgcgggtgaa                                                                                                                                       >  1:163954‑M2/1‑137 (MQ=255)
               gACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCattcgcgggtgaaccc                                                                                                                                    >  1:2191551‑M2/1‑134 (MQ=255)
                 ccAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCattcgcgggtgaacccgc                                                                                                                                  >  1:1292669‑M2/1‑132 (MQ=255)
                  cAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGCACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCattcgcgggtgaacccgct                                                                                                                                 <  1:1453943‑M2/150‑20 (MQ=255)
                          gTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTCCCCCCGCGAATACGATGGTTGACCCACCGACGCattcgcgggtgaacccgctcctgcagg                                                                                                                         >  1:635094‑M2/1‑123 (MQ=255)
                              gACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCCGGAGCGGGTTAACACGCGCCTACGAGCATGATCCCCGCGACGGACTTGCCGGGGAACGCGCTGCTGCAGGGTTc                                                                                                                     <  2:24638/150‑1 (MQ=255)
                                catcatCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCTAACTGCCGTTCGCTCGAAGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGc                                                                                                                   <  1:1780074/150‑1 (MQ=255)
                                   catcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGcc                                                                                                                  >  1:1150392/1‑148 (MQ=255)
                                   catcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGcc                                                                                                                  <  2:1150392/148‑1 (MQ=255)
                                     tcgtcgCAATCTAGCGCAAGAACCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATCAGTGGCCTGCAGGAGCTGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGTGTGAACCCGCTCCTGCAGGGTGCGCCGTCa                                                                                                              <  2:1004159/150‑1 (MQ=38)
                                          gCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCtt                                                                                                         <  2:2152122/150‑1 (MQ=255)
                                                  ccGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACt                                                                                                 >  2:2206607/1‑150 (MQ=255)
                                                    gCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCattcgcgg                                                                                                                                            >  2:1055093‑M2/1‑97 (MQ=255)
                                                    gCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCattcgcgg                                                                                                                                            <  1:1055093‑M2/105‑9 (MQ=255)
                                                       aGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACg                                                                                            >  2:469911/1‑150 (MQ=255)
                                                         aCCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGTCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCAACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAt                                                                                          <  1:2171208/150‑1 (MQ=255)
                                                         aCCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGATGCAGTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAt                                                                                          <  1:2037188/150‑1 (MQ=255)
                                                            cAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCACGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGAc                                                                                      >  2:1327785/1‑150 (MQ=255)
                                                                   cgcgcGACCTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGGTCGGGCTACTCCGCGAATACGAGGGCGAACCCACCCACGCATTCGCTGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGACTTTGATTACTTGACGATGACCTTGCCt                                                                               <  2:428486/149‑1 (MQ=255)
                                                                    gcgcGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCtt                                                                              <  1:787229/150‑1 (MQ=255)
                                                                       cgAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                                           <  1:2206607/150‑1 (MQ=255)
                                                                         aaCTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtcatc                                                                         >  1:1393080/1‑150 (MQ=255)
                                                                            cgCCGTGCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                                      <  1:1887981/149‑1 (MQ=255)
                                                                               cGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCgag                                                                   <  2:2236460/150‑1 (MQ=255)
                                                                                      tCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCattcgcgggtgaacccgctcct                                                                                                                              <  2:1138656‑M2/85‑23 (MQ=255)
                                                                                      tCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCattcgcgggtgaacccgctcct                                                                                                                              >  1:1138656‑M2/1‑63 (MQ=255)
                                                                                      tCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGc                                                            >  2:266007/1‑150 (MQ=255)
                                                                                        cTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCt                                                                                                          <  1:1624657/103‑1 (MQ=255)
                                                                                        cTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                                           >  2:1690544/1‑133 (MQ=255)
                                                                                        cTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                                           <  1:1690544/133‑1 (MQ=255)
                                                                                         ttGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACtt                                                                                                >  2:2106462/1‑112 (MQ=255)
                                                                                         ttGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACtt                                                                                                <  1:2106462/112‑1 (MQ=255)
                                                                                          ggATAAGTCGCGTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAGCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCgg                                                        <  2:1218539/149‑1 (MQ=255)
                                                                                                 tGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTg                                                 >  2:2327967/1‑150 (MQ=255)
                                                                                                    ccTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCa                                               <  2:2274318/149‑1 (MQ=255)
                                                                                                    ccTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCa                                               >  1:2274318/1‑149 (MQ=255)
                                                                                                     cTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCaga                                             <  1:266007/150‑1 (MQ=255)
                                                                                                       gatctaGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaag                                           <  1:1094142/145‑1 (MQ=255)
                                                                                                         aGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                                      >  1:942715/1‑121 (MQ=14)
                                                                                                         aGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                                      <  2:942715/121‑1 (MQ=14)
                                                                                                          ggAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCGGGGAAAGTGCAGaagaag                                        <  1:80739/150‑1 (MQ=255)
                                                                                                           gAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGTCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                                      <  1:1886881/119‑1 (MQ=11)
                                                                                                           gAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                                      >  2:1886881/1‑119 (MQ=14)
                                                                                                           gAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTTAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGc                                       >  2:1413943/1‑150 (MQ=255)
                                                                                                            aGCGGGTTCACCCGAGAATAACAGGGTAAGCCCAACTACGCATTTACGGGTGTATCTTCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCt                                      <  2:256063/150‑1 (MQ=255)
                                                                                                             gagggTTCCACCGCGAATACGAGGTGGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGACTTTGATTAATCAGTCATCACTCGACCCGGTATCATCGGGATAGTGCCCGGGAAGGTGCAGaagaagatc                                     >  2:1094142/3‑147 (MQ=16)
                                                                                                             gCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAt                                                                                          >  1:136044/1‑98 (MQ=255)
                                                                                                             gCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAt                                                                                          <  2:136044/98‑1 (MQ=255)
                                                                                                                ggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAg                                  >  2:1986424/1‑150 (MQ=255)
                                                                                                                ggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAg                                  >  1:1793524/1‑150 (MQ=255)
                                                                                                                 gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGAtgt                                                               >  2:806881/1‑120 (MQ=255)
                                                                                                                 gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGAtgt                                                               <  1:806881/120‑1 (MQ=255)
                                                                                                                 gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAg                                                    >  1:408091/1‑131 (MQ=255)
                                                                                                                 gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTg                                     >  1:1800680/1‑146 (MQ=255)
                                                                                                                 gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTg                                     <  2:1800680/146‑1 (MQ=255)
                                                                                                                    cACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcac                              >  1:2180530/1‑150 (MQ=255)
                                                                                                                      cccGCGAATACGAGGGTGAACCCACCGACGCattcgcgggtgaacccgctcctgcaggccacttatcaaggagcgaacggcagttcgcgcttgtgctgggtcttgcggtaggtggcgacgatgatgtcgaacgctgcctggtcaaccggc                             <  2:1007829‑M2/150‑120 (MQ=255)
                                                                                                                      cccGCGAATACGAGGGTGAACCCACCGACGCattcgcgggtgaacccgctcctgcaggccacttatcaaggagcgaacggcagttcgcgcttgtgctgggtcttgcggtaggtggcgacgatgatgtcgaacgctgcctggtcaaccggc                             <  1:1673764‑M2/150‑120 (MQ=255)
                                                                                                                      cccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTACACTACTTGACGATGACCTTGCATTTAGTTCACTCGACGCATCCCGGGAGGGTGCAGAAGAAGCTGTAGTcacca                            >  1:1103459/1‑150 (MQ=255)
                                                                                                                       gcgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGGGCAGAAGAAGCTGTAGTcaccac                           >  1:2088836/2‑150 (MQ=255)
                                                                                                                       ccGCGAATACGAGGGTGAACCCACCGACGCattcgcgggtgaacccgctcctgcaggccacttatcaaggagcgaacggcagttcgcgcttgtgctgggtcttgcggtaggtggcgacgatgatgtcgaacgctgcctggtcaaccggct                            <  2:3900‑M2/150‑121 (MQ=255)
                                                                                                                       ccGCGAATACGAGGGTGAACCCACCGACGCattcgcgggtgaacccgctcctgcaggccacttatcaaggagcgaacggca                                                                                                 >  1:1828026‑M2/1‑30 (MQ=255)
                                                                                                                       ccGCGAATACGAGGGTGAACCCACCGACGCattcgcgggtgaacccgctcctgcaggccacttatcaaggagcgaacggca                                                                                                 <  2:1828026‑M2/81‑52 (MQ=255)
                                                                                                                       ccGCGAATACGAGGGTGAACCCACCGACGCattcgcgggtgaacccgctcctgcaggccacttatcaagga                                                                                                           <  1:1749568‑M2/71‑42 (MQ=255)
                                                                                                                       ccGCGAATACGAGGGTGAACCCACCGACGCattcgcgggtgaacccgctcctgcaggccacttatcaagga                                                                                                           >  2:1749568‑M2/1‑30 (MQ=255)
                                                                                                                       ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATg                                                                                         <  2:1242393/89‑1 (MQ=255)
                                                                                                                       ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                           >  1:1617942/1‑150 (MQ=255)
                                                                                                                       ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                           >  1:240244/1‑150 (MQ=255)
                                                                                                                        cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGACTTTGATTCCTTGACCCTGACCCTGGCCTTCGCTGAAGGGGAGTGGCCCGGGGAGGTGCAGGGGAAGGGGTAGTGgccacc                          >  2:509988/1‑150 (MQ=16)
                                                                                                                         gcgAATACGAGGGTGAACCCACCGACGCattcgcgggtgaacccgctcctgcaggccacttatcaaggagcgaacggcagttcgcgcttgtgctgggtcttgcggtaggtggcgacgatgatgtcgaacgctgcctggtcaaccggctgc                          <  1:1202077‑M2/150‑123 (MQ=255)
                                                                                                                         gcgAATACGAGGGTGAACCCACCGACGCattcgcgggtgaacccgctcctgcaggccacttatcaaggagcgaacggcagttcgcgcttgtgctgggtcttgcggtaggtggcgacgatgat                                                      <  1:1225340‑M2/122‑95 (MQ=255)
                                                                                                                         gcgAATACGAGGGTGAACCCACCGACGCattcgcgggtgaacccgctcctgcaggccacttatcaaggagcgaacggcagttcgcgcttgtgctgggtcttgcggtaggtggcgacga                                                          >  2:1480138‑M2/1‑28 (MQ=255)
                                                                                                                         gcgAATACGAGGGTGAACCCACCGACGCattcgcgggtgaacccgctcctgcaggccacttatcaaggagcgaacggcagttcgcgcttgtgctgggtcttgcggtaggtggcgacga                                                          <  1:1480138‑M2/118‑91 (MQ=255)
                                                                                                                         gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGt                                                  <  2:1565835/125‑1 (MQ=255)
                                                                                                                         gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGt                                                  >  1:1565835/1‑125 (MQ=255)
                                                                                                                           gAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGa                                                                                        <  1:774598/86‑1 (MQ=255)
                                                                                                                           gAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGa                                                                                        >  2:774598/1‑86 (MQ=255)
                                                                                                                           gAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGc                       <  2:860948/150‑1 (MQ=255)
                                                                                                                             aTACGAGGGTGAACCCACCGACGCattcgcgggtgaacccgctcctgcaggccacttatcaaggagcgaacggcagttcgcgcttgtgctgggtcttgcggtaggtggcgacgatgatgtcgaacgctgcctggtcaaccggctgcccat                      <  2:313318‑M2/150‑127 (MQ=255)
                                                                                                                              tACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCga                                                                    <  2:1013430/102‑1 (MQ=255)
                                                                                                                              tACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCga                                                                    >  1:1013430/1‑102 (MQ=255)
                                                                                                                              tACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTc                    <  2:1393080/150‑1 (MQ=255)
                                                                                                                              tACGAGGGTGAACCCACCCTCGTattcgcgggtgaacccgctcctgcaggccacttatcaaggagcgaacggcagttcgcgcttgtgctgggtcttgcggtaggtggcgacgatgatgtcgaacgctgcctggtcaaccggctgcccatg                     <  1:1297577‑M2/150‑128 (MQ=255)
                                                                                                                                 gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGCGATGTGGCCCTGTAAGATGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg                 >  2:891771/1‑150 (MQ=255)
                                                                                                                                 gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg                 <  2:1453943/150‑1 (MQ=255)
                                                                                                                                  aGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAgg                                                   >  2:1079293/1‑115 (MQ=255)
                                                                                                                                  aGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAgg                                                   <  1:1079293/115‑1 (MQ=255)
                                                                                                                                    ggTGAACCCACCGACGCattcgcgggtgaacccgctcctgcaggccacttatcaaggagcgaacggcagttcgcgcttgtgctgggtcttgcggtaggtggcgacgatgatgtcgaacgctgcctggtcaaccggctgcccatgcaggaa               >  1:313318‑M2/1‑17 (MQ=255)
                                                                                                                                    ggTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCtt              >  2:730659/1‑150 (MQ=255)
                                                                                                                                           ccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCgaga                                                                  >  2:2153843/1‑91 (MQ=255)
                                                                                                                                           ccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCgaga                                                                  <  1:2153843/91‑1 (MQ=255)
                                                                                                                                            cACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCa                                               <  1:1544129/109‑1 (MQ=255)
                                                                                                                                            cACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCa                                               >  2:1544129/1‑109 (MQ=255)
                                                                                                                                             aCCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCACCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTc     >  1:631425/1‑150 (MQ=255)
                                                                                                                                                gACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGaa  >  2:385264/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
CATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGAA  >  NC_002947/5536979‑5537272

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