Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,621,920 G→A 5.5% L264L (CTG→CTA) ‡ PP_1421 → sensor histidine kinase TctE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,621,9200GA 5.5% 263.5 / 4.4 90L264L (CTG→CTA) ‡PP_1421sensor histidine kinase TctE
Reads supporting (aligned to +/- strand):  ref base G (45/40);  new base A (3/2);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.01e-01

CTTGCGAGCGTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCGGGTCGAGAATG  >  NC_002947/1621773‑1622065
                                                                                                                                                   |                                                                                                                                                 
cTTGCGAGCGTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTggc                                                                                                                                                 <  2:1472094/150‑1 (MQ=255)
     gAGCGTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGt                                                                                                                                            >  1:1654743/1‑150 (MQ=255)
      aGCGTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGtt                                                                                                                                           <  1:296405/150‑1 (MQ=255)
         gTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCGCAGCGCATAGCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGaa                                                                                                                                        <  2:1928000/150‑1 (MQ=255)
         gTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGaa                                                                                                                                        <  1:847485/150‑1 (MQ=255)
              cGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCAc                                                                                                                                   <  1:1447627/150‑1 (MQ=255)
              cGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCAc                                                                                                                                   >  2:1261128/1‑150 (MQ=255)
               ggtggtGGAGGTGCAGCGCGAACTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACTAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACg                                                                                                                                  >  2:283923/1‑150 (MQ=255)
                gtggtgCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACgt                                                                                                                                 >  1:1191591/1‑150 (MQ=255)
                     gcaggtgcagCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGa                                                                                                                            <  2:540805/150‑1 (MQ=255)
                      caggtgcagCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAg                                                                                                                           <  2:1419352/150‑1 (MQ=255)
                         gtgcagCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGCGAGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCATGTGTCGAGCTg                                                                                                                        >  2:377457/1‑150 (MQ=255)
                          tgcagCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTgg                                                                                                                       >  2:535986/1‑150 (MQ=255)
                           gcagCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTCCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTggg                                                                                                                      >  1:326960/1‑150 (MQ=255)
                            cagCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGc                                                                                                                     <  1:1015340/150‑1 (MQ=255)
                            cagCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGc                                                                                                                     <  2:1345285/150‑1 (MQ=255)
                            cagCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGc                                                                                                                     <  2:1515201/150‑1 (MQ=255)
                            cagCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTGCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGc                                                                                                                     >  2:943908/1‑150 (MQ=255)
                              gcgcGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCtt                                                                                                                   >  1:714233/1‑150 (MQ=255)
                                gcgAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTgc                                                                                                                 <  2:1161791/150‑1 (MQ=255)
                                gcgAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTCGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTgc                                                                                                                 <  2:1299550/150‑1 (MQ=255)
                                    gCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACTCCTCTGGCGGCGGTGATGGCACGTGTCGAGCTGGGCTTGCGCTc                                                                                                             >  2:804200/1‑150 (MQ=255)
                                    gCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCt                                                                                                                         <  1:679619/138‑1 (MQ=255)
                                    gCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCt                                                                                                                         >  2:679619/1‑138 (MQ=255)
                                      tGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCgg                                                                                                           <  1:1209839/150‑1 (MQ=255)
                                        ggCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGGGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGcc                                                                                                         <  2:98817/150‑1 (MQ=255)
                                          cccGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGa                                                                                                       <  1:1281932/150‑1 (MQ=255)
                                           ccGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAg                                                                                                      >  1:413886/1‑150 (MQ=255)
                                            cGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGc                                                                                                     <  2:1654743/150‑1 (MQ=255)
                                                     gcgcCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTgc                                                                                                                 >  2:1458598/1‑129 (MQ=255)
                                                     gcgcCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTgc                                                                                                                 <  1:1458598/129‑1 (MQ=255)
                                                       gcCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTgg                                                                                          >  2:1063375/1‑150 (MQ=255)
                                                        ccTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCtt                                                                                                                   >  1:407827/1‑124 (MQ=255)
                                                        ccTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCtt                                                                                                                   <  2:407827/124‑1 (MQ=255)
                                                        ccTTGAACCACTTTACCGAGCGCTTGCGAGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGCACGCCGCTGGCGGCGTTGCAGGCATGTGTCGATCTAGGCTTGCGCTCGGCCGAGCCGCAGGAGTggc                                                                                         >  2:908395/1‑150 (MQ=255)
                                                         cTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGGACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTggcg                                                                                        >  2:1310757/1‑150 (MQ=255)
                                                         cTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTggcg                                                                                        >  2:413473/1‑150 (MQ=255)
                                                         cTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTggcg                                                                                        >  1:905870/1‑150 (MQ=255)
                                                         cTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCCAGCGGCGttc                                                                                                                                          <  2:428841/100‑2 (MQ=255)
                                                         cTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCCAGCGGCGttc                                                                                                                                          >  1:428841/1‑99 (MQ=255)
                                                          ttGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCg                                                                                                            <  2:360068/129‑1 (MQ=255)
                                                          ttGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCg                                                                                                            >  1:360068/1‑129 (MQ=255)
                                                               ccACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGAc                                                                                  >  1:1352069/1‑150 (MQ=255)
                                                                cACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGAcc                                                                                 >  2:166686/1‑150 (MQ=255)
                                                                 aCTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGACCGCCGCTGGCGGCGGTGAAGGCACGTGTCGAGCTGGGCTTGCGCCCGGCCGAGCCGCAGGAGCGGCGGCAGACCt                                                                                >  2:717733/1‑150 (MQ=255)
                                                                  cTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTc                                                                                                             >  1:1779485/1‑120 (MQ=255)
                                                                  cTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTc                                                                                                             <  2:1779485/120‑1 (MQ=255)
                                                                  cTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCtt                                                                               <  2:714233/150‑1 (MQ=255)
                                                                   tttACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTg                                                                              <  1:1001516/150‑1 (MQ=255)
                                                                   tttACCGAGCGCTTGCGCGGCCAGTCCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGGCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTg                                                                              <  2:326960/150‑1 (MQ=255)
                                                                       ccGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGt                                                                          >  2:1488459/1‑150 (MQ=255)
                                                                        cGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTc                                                                         >  1:15140/1‑150 (MQ=255)
                                                                         gAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCg                                                                        <  2:608990/150‑1 (MQ=255)
                                                                          aGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTg                                                                              >  1:855504/1‑143 (MQ=255)
                                                                          aGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTg                                                                              <  2:855504/143‑1 (MQ=255)
                                                                          aGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCgg                                                                       <  1:166686/150‑1 (MQ=255)
                                                                           gcgcTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGc                                                                      <  2:15140/150‑1 (MQ=255)
                                                                              cTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGa                                                                                   >  1:889159/1‑134 (MQ=255)
                                                                              cTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGa                                                                                   <  2:889159/134‑1 (MQ=255)
                                                                              cTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGc                                                                   >  2:1417819/1‑150 (MQ=255)
                                                                               ttGCGCGGCCAGTTCTAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGcc                                                                  <  2:1352069/150‑1 (MQ=255)
                                                                                tGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGTCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGccc                                                                 <  1:1297931/150‑1 (MQ=255)
                                                                                     ggCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGaa                                                                                                                                        >  2:1388807/1‑74 (MQ=255)
                                                                                     ggCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGaa                                                                                                                                        <  1:1388807/74‑1 (MQ=255)
                                                                                              gAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAgg                                                   >  1:1529905/1‑150 (MQ=255)
                                                                                                gcgcCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCTGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGAAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCt                                                 >  1:1408649/1‑150 (MQ=255)
                                                                                                gcgcCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCt                                                 >  1:1328939/1‑150 (MQ=255)
                                                                                                     aGGCGCAGTTCATTTCCAGCGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCCGGAGAGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGAccc                                            >  2:1148389/1‑150 (MQ=255)
                                                                                                       gcgcAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCAc                                          >  1:1839913/1‑150 (MQ=255)
                                                                                                                  tGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCa                               <  1:1488459/150‑1 (MQ=255)
                                                                                                                  tGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCa                               <  1:1201936/150‑1 (MQ=255)
                                                                                                                   gCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACCGGCTGACCCACCTGGCCAACCAg                              >  2:303088/1‑150 (MQ=255)
                                                                                                                    ccGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGGGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGc                             >  2:1479510/1‑150 (MQ=255)
                                                                                                                            gCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTcc                     >  1:1289900/1‑150 (MQ=255)
                                                                                                                            gCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTcc                     <  1:413473/150‑1 (MQ=255)
                                                                                                                              ccACGAACTGCGAACGCCGCTGGTGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCt                   >  2:1938957/1‑150 (MQ=255)
                                                                                                                               cACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTg                  >  1:293268/1‑150 (MQ=255)
                                                                                                                                aCGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCATGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTgg                 <  1:377457/150‑1 (MQ=255)
                                                                                                                                   aaCTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGgcgc              <  1:1600135/150‑1 (MQ=255)
                                                                                                                                     cTGCGAACGCCGCTGGCGGCGTTGAAGCCCCGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCTGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCCCCCGTTGTCCCTGGCGCgg            >  1:1299144/1‑150 (MQ=255)
                                                                                                                                      tGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGc                                                  <  2:289383/111‑1 (MQ=255)
                                                                                                                                      tGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGc                                                  >  1:289383/1‑111 (MQ=255)
                                                                                                                                       ttcaaCGCCGCCAGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCAc                                                          >  2:1470750/4‑102 (MQ=255)
                                                                                                                                       ttcaaCGCCGCCAGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCAc                                                          <  1:1470750/99‑1 (MQ=255)
                                                                                                                                       ttcaaCGCCGCCAGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCGGGt          >  2:10199/4‑150 (MQ=255)
                                                                                                                                        cGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGAAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCGGGTc         <  2:1408649/150‑1 (MQ=255)
                                                                                                                                         gAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCGGGTCg        >  2:1604456/1‑150 (MQ=255)
                                                                                                                                          aaCGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCGGGTCga       >  2:363390/1‑150 (MQ=255)
                                                                                                                                           aCGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCGGGTCgag      <  1:363390/150‑1 (MQ=255)
                                                                                                                                               cgcTGGCGGCGTTTAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCGGGTCGAGAATg  <  2:106317/150‑1 (MQ=255)
                                                                                                                                               cgcTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCGGGTCGAGAATg  >  1:669657/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTTGCGAGCGTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCGGGTCGAGAATG  >  NC_002947/1621773‑1622065

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: