Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,967,284 A→G 8.0% intergenic (+111/+8) PP_1762 → / ← yciK hypothetical protein/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,967,2840AG 8.0% 256.2 / 8.4 88intergenic (+111/+8)PP_1762/yciKhypothetical protein/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (46/34);  new base G (4/3);  total (51/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTTCTCGCCGGGGTAGGCCTGGGCCCGCATCGCCGTGCG  >  NC_002947/1967193‑1967432
                                                                                           |                                                                                                                                                    
tctTCGCGGGCTTGCCCGCTCCCAAAGGATCACCGCAACCCTCAAGGATTGCGCCCCGCCTGCGGGAGCGGGCAAGCCCGCGAAGAGGCCGCCGCAGGTTTATTGGGCGTTGAACGCCATCCCGTTCACCCCTGTGCTTTTCGGCCCCac                                                                                            >  2:500936/1‑149 (MQ=37)
       GGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtac                                                                                     <  1:1472401‑M1/150‑1 (MQ=255)
           ttgcccgctcccacaggcagggcgcaatccttgagggttgcggtgatcctgtgggagcgggcaagcccgcgaagaggccgGTGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtaca                                                                                    <  1:1631677‑M2/67‑1 (MQ=255)
           ttgcccgctcccacaggcagggcgcaatccttgagggttgcggtgatcctgtgggagcgggcaagcccgcgaagaggccgGTGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtaca                                                                                    >  2:1631677‑M2/81‑147 (MQ=255)
           TTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacaggt                                                                                 <  1:587847‑M1/150‑1 (MQ=255)
           TTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacag                                                                                   <  2:1283162‑M1/148‑1 (MQ=255)
           TTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacag                                                                                   >  1:1283162‑M1/1‑148 (MQ=255)
              cccgctcccacaggcagggcgcaatccttgagggttgcggtgatcctgtgggagcgggcaagcccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAggtacaggtaca                                                                              >  1:1265743‑M2/78‑150 (MQ=255)
              CCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtac                                                                                     <  1:2218445‑M1/143‑1 (MQ=255)
              CCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtac                                                                                     >  2:2218445‑M1/1‑143 (MQ=255)
               ccgctcccacaggcagggcgcaatccttgagggttgcggtgatcctgtgggagcgggcaagcccgcgaagaggccgGTACAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGTTACAGGTacac                                                                             >  1:1099939‑M2/77‑150 (MQ=255)
               ccgctcccacaggcagggcgcaatccttgagggttgcggtgatcctgtgggagcgggcaagcccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCt                                                                                                                           <  2:1781667‑M2/28‑1 (MQ=255)
               ccgctcccacaggcagggcgcaatccttgagggttgcggtgatcctgtgggagcgggcaagcccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCt                                                                                                                           >  1:1781667‑M2/77‑104 (MQ=255)
               CCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac                                                                             <  1:254951‑M1/150‑1 (MQ=255)
                cgctcccacaggcagggcgcaatccttgagggttgcggtgatcctgtgggagcgggcaagcccgcgaagaggccgGTGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACAcc                                                                            >  1:2225097‑M2/76‑150 (MQ=255)
                  ctcccacaggcagggcgcaatccttgagggttgcggtgatcctgtgggagcgggcaagcccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCgg                                                                          >  1:2335103‑M2/74‑150 (MQ=255)
                   tcccacaggcaggacgcaatccttgagggttgcggtgatcctgtgggagcgggcaagcccgcgaagaggccgGTACAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGc                                                                         >  2:759060‑M2/73‑150 (MQ=255)
                   TCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTTCCCCCTCCGGATAACCGCCGCATCAGGTACAGGTACACCGGc                                                                         <  2:1230681‑M1/150‑1 (MQ=255)
                   TCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGc                                                                         >  2:587847‑M1/1‑150 (MQ=255)
                   TCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGc                                                                         <  1:1157527‑M1/150‑1 (MQ=255)
                     ccacaggcagggcgcaatccttgagggttgcggtgatcctgtgggagcgggcaagcccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGcag                                                                       >  1:1612037‑M2/71‑149 (MQ=255)
                     ccacaggcagggcgcaatccttgagggttgcggtgatcctgtgggagcgggcaagcccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCAt                                                                       >  2:2336351‑M2/71‑150 (MQ=255)
                     CCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCAt                                                                       <  2:1239867‑M1/150‑1 (MQ=255)
                     CCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCAt                                                                       <  2:9837‑M1/150‑1 (MQ=255)
                      cacaggcagggcgcaatccttgagggttgcgccctgcctgtgggagcgggcaagcccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATg                                                                      >  1:381356‑M2/70‑150 (MQ=255)
                      CACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTAtt                                                                                                                                          <  1:1029667‑M1/82‑1 (MQ=255)
                      CACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTAtt                                                                                                                                          >  2:1029667‑M1/1‑82 (MQ=255)
                      CACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCGTGGACTGTCCGGCCCCATCAGGTACAGGTACACCGGCATg                                                                      <  1:22261‑M1/150‑1 (MQ=255)
                       acaggcagggcgcaatccttgagggttgcggtgatcctgtgggagcgggcaagcccgcgaagaggccgGTGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa                                                                     <  2:2225097‑M2/82‑1 (MQ=255)
                       ACAGGATCACCGCCACCCTCAAGGCTTGCGCCGTCCCTGTGTGATCGGGCATGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa                                                                     >  2:22261‑M1/1‑150 (MQ=255)
                       ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa                                                                     >  1:1239867‑M1/1‑150 (MQ=255)
                       ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTAGGATCGTGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa                                                                     <  2:559496‑M1/150‑1 (MQ=255)
                         aggcagggcgcaatccttgagggttgcggtgatcctgtgggagcgggcaagcccgcgaagaggccgGTACAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGAtc                                                                   <  2:1099939‑M2/84‑1 (MQ=255)
                           gcagggcgcaatccttgagggttgcggtgatcctgtgggagcgggcaagcccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGAtctc                                                                 >  2:2121616‑M2/65‑150 (MQ=255)
                             tCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCtt                                                               >  1:2326790/1‑150 (MQ=255)
                              cACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCAccc                                                                                                               >  2:2029671/1‑101 (MQ=255)
                              cACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCAccc                                                                                                               <  1:2029671/101‑1 (MQ=255)
                                ccGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCAGCATGATCTCTTCgg                                                            >  2:1677230/1‑150 (MQ=255)
                                 ctcAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCggg                                                           >  2:1553175/3‑150 (MQ=255)
                                 cGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCggg                                                           >  2:1424700/1‑150 (MQ=255)
                                  gCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGaa                                                                                                                                >  1:959592/1‑80 (MQ=255)
                                  gCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGaa                                                                                                                                <  2:959592/80‑1 (MQ=255)
                                        cttgagggttgcggtgatcctgtgggagcgggcaagcccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGcggc                                                    <  2:1612037‑M2/99‑2 (MQ=255)
                                        cTCAAGGATGGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCggg                                                    <  1:500936/150‑1 (MQ=255)
                                         tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtaca                                                                                    >  1:72682/1‑117 (MQ=255)
                                         tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtaca                                                                                    <  2:72682/117‑1 (MQ=255)
                                         tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCgg                                                            >  2:306699/1‑141 (MQ=255)
                                         tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCgg                                                            <  1:306699/141‑1 (MQ=255)
                                         tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGgt                                                   >  1:772203/1‑150 (MQ=255)
                                         tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACAGGCATGATCTCTTCGGGCAGCGGgt                                                   >  1:1697808/1‑150 (MQ=255)
                                         tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCATGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGgt                                                   >  1:227864/1‑150 (MQ=255)
                                             ggATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGgttgtt                                               >  1:1200213/1‑150 (MQ=255)
                                              gATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTc                                              >  2:2060564/1‑150 (MQ=255)
                                                 tGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGc                                           <  1:9837/150‑1 (MQ=255)
                                                   cggtgatcctgtgggagcgggcaagcccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCgg                                         <  1:2121616‑M2/110‑1 (MQ=255)
                                                    gcCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCggg                                        >  2:1245207/1‑150 (MQ=255)
                                                     cccTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGt                                       >  2:189656/1‑150 (MQ=255)
                                                      cctgcctgTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCgg                                         >  2:559185/1‑147 (MQ=255)
                                                      cctgcctgTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCgg                                         <  1:559185/147‑1 (MQ=255)
                                                      cctgcctgTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCAGGGTCTTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGtt                                      >  1:419204/1‑150 (MQ=255)
                                                       ctgcctgTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTtc                                     >  1:1895197/1‑150 (MQ=255)
                                                       ctgcctgTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTAGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTTTTCTGCGGGTtc                                     >  1:966687/1‑150 (MQ=255)
                                                        tgcctgTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTtct                                    >  2:2387211/1‑150 (MQ=255)
                                                         gcctgTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTtctc                                   <  2:2326790/150‑1 (MQ=255)
                                                          cctgtgggagcgggcaagcccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTTCTCg                                  >  1:1220914‑M2/34‑150 (MQ=255)
                                                              tgggagcgggcaagcccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTTCTCGCCgg                              <  2:254951‑M2/121‑1 (MQ=255)
                                                                  agcgggcaagcccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTTCTCGCCGGGGTa                          >  2:1695109‑M2/26‑150 (MQ=255)
                                                                   gcgggcaagcccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTTCTCGCCGGGGTAg                         <  1:1424700‑M2/126‑1 (MQ=255)
                                                                       gcaagcccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCAcc                                                                                                                >  1:1758616‑M2/21‑59 (MQ=255)
                                                                       gcaagcccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCAcc                                                                                                                <  2:1758616‑M2/39‑1 (MQ=255)
                                                                            cccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTTCTCGCCGGGGTAGGCCTGGGcc                >  2:1296920‑M2/16‑150 (MQ=255)
                                                                            cccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTTCTCGCCGGGGTAGGCCTGGGcc                >  2:419152‑M2/16‑150 (MQ=255)
                                                                             ccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggt                                                                                       <  1:1359997‑M2/64‑1 (MQ=255)
                                                                             ccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggt                                                                                       >  2:1359997‑M2/15‑78 (MQ=255)
                                                                             ccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCgg                                                     >  1:1512869‑M2/15‑112 (MQ=255)
                                                                             ccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCgg                                                     <  2:1512869‑M2/98‑1 (MQ=255)
                                                                             ccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTTCTCGCCGGGGTAGGcctgggagt               <  2:1157527‑M2/136‑4 (MQ=255)
                                                                             ccgcgaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTTCTCGCCGGGGTAGGCCTGGGccc               >  1:956616‑M2/15‑150 (MQ=255)
                                                                              cgcgaagaggccgATGCAGGTTTATTGGGCGTTGCACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTTCTCGCCGGGGTAGGCCTGGGCCCg              >  2:435425‑M2/14‑150 (MQ=255)
                                                                                 gaagaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTTCTCGCCGGGGTAGGCCTGGGCCCGCAt           <  2:956616‑M2/140‑1 (MQ=255)
                                                                                    gaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCtgt                                                                                                           >  1:2222452‑M2/8‑51 (MQ=255)
                                                                                    gaggccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCtgt                                                                                                           <  2:2222452‑M2/44‑1 (MQ=255)
                                                                                       gccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCgg                                                            <  1:1506567‑M2/91‑1 (MQ=255)
                                                                                       gccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCgg                                                            >  2:1506567‑M2/5‑95 (MQ=255)
                                                                                       gccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTTCTCGCCGGGGTAGGCCTGGGCCCGCATCGCCGt     <  1:1553175‑M2/146‑1 (MQ=255)
                                                                                       gccgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCATGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTTCTCGCCGGGGTAGGCCTGGGCCCGCATCGCCGt     <  2:227864‑M2/146‑1 (MQ=255)
                                                                                         cgATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTAACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTTCTCGCCGGGGTAGGCCTGGGCCCGCATCGCCGTgc   >  2:273696‑M2/3‑150 (MQ=255)
                                                                                          gATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTTCTCGCCGGGGTAGGCCTGGGCCCGCATCGCCGTgcg  >  1:2388811‑M2/2‑150 (MQ=255)
                                                                                           |                                                                                                                                                    
TCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTTGTTCTGCGGGTTCTCGCCGGGGTAGGCCTGGGCCCGCATCGCCGTGCG  >  NC_002947/1967193‑1967432

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: