Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,855,731 C→T 7.2% intergenic (+48/+34) PP_3401 → / ← PP_3402 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,855,7310CT 7.2% 589.0 / 23.7 195intergenic (+48/+34)PP_3401/PP_3402hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (88/93);  new base T (5/9);  total (93/102)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 196 total aligned reads displayed.

AGTGTGGCCTTCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGTTGAGG  >  NC_002947/3855582‑3855879
                                                                                                                                                     |                                                                                                                                                    
  tgtgGCCTTCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCa                                                                                                                                                    <  1:5345163/150‑1 (MQ=255)
    tgGCCTTCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGc                                                                                                                                                  <  2:4874249/150‑1 (MQ=255)
     ggCCTTCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGcc                                                                                                                                                 <  1:5214021/150‑1 (MQ=255)
     ggCCTTCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGcc                                                                                                                                                 >  1:531382/1‑150 (MQ=255)
      gCCTTCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCt                                                                                                                                                <  2:3495971/150‑1 (MQ=255)
      gCCTTCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCt                                                                                                                                                <  2:2990629/150‑1 (MQ=255)
      gCCTTCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCt                                                                                                                                                <  1:885814/150‑1 (MQ=255)
       ccTTCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTg                                                                                                                                               <  2:4117968/150‑1 (MQ=255)
        cTTCGTCTGGTATGGGGGGTGTTATCACGGCGTTTTTCGCGAGATGGTAATGTGGGGCGAGGGTGGCCTTAACGTGCTGCTGGAACTGCAATAAGCCGATCAGGGGGGGCCGTGGCCACGCTCGTGGGTGTTTATACCGGGCAAGCCTgg                                                                                                                                              <  2:1438921/150‑1 (MQ=255)
        cTTCGTCTGGGATGGGGGGTGTTATCACGGCGTGTTTCGCGAGATGGTAATGTGCGGCGAGGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGGGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTgg                                                                                                                                              <  2:1438827/150‑1 (MQ=255)
        cTTCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTgg                                                                                                                                              >  2:637451/1‑150 (MQ=255)
        cTTCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTgg                                                                                                                                              >  2:177771/1‑150 (MQ=255)
         ttCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTggg                                                                                                                                             <  2:1004638/150‑1 (MQ=255)
            ggcggggATGGGCGGAGTTAGCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCAGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGGGGCTGCGTATACCCGGCCAGCCTGGGCTg                                                                                                                                          >  2:2675308/5‑150 (MQ=255)
             tCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGa                                                                                                                                         >  2:1627595/1‑150 (MQ=255)
                 ggATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGTAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGTCAGCCCAGGCTGGCCgg                                                                                                                                     <  2:3879592/150‑1 (MQ=255)
                 ggATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGTCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCgg                                                                                                                                     <  1:664270/150‑1 (MQ=255)
                 ggATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCgg                                                                                                                                     >  2:5232505/1‑150 (MQ=255)
                 ggATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCgg                                                                                                                                     >  1:911628/1‑150 (MQ=255)
                 ggATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCgg                                                                                                                                     <  2:3085302/150‑1 (MQ=255)
                  gATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCggg                                                                                                                                    >  2:3613674/1‑150 (MQ=255)
                   aTGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGTCAGCCCAGGCTGGCCgggt                                                                                                                                   <  2:1218749/150‑2 (MQ=255)
                    tGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGc                                                                                                                                  <  1:2239129/150‑1 (MQ=255)
                       gCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGTCAGCCCAGGCTGGCCgggtataa                                                                                                                               <  2:211145/150‑6 (MQ=255)
                       gCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCgggtataa                                                                                                                               <  2:4316826/150‑6 (MQ=255)
                        cGGAGTTATCACGGCGGGGTTCGCGAGAAGGGACAGTGCGGCGACGGTGGCCTGAACCGGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGa                                                                                                                              <  2:3047412/150‑1 (MQ=255)
                         ggAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGAc                                                                                                                             <  2:4142280/150‑1 (MQ=255)
                           aGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAg                                                                                                                           <  1:763656/150‑1 (MQ=255)
                           aGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAg                                                                                                                           <  2:3349699/150‑1 (MQ=255)
                             ttATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGgc                                                                                                                         >  2:2012999/1‑150 (MQ=255)
                             ttATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGgc                                                                                                                         >  1:4977042/1‑150 (MQ=255)
                              tATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGgcg                                                                                                                        >  2:143082/1‑150 (MQ=255)
                               aTCACGGCGTGGTTCGCGAGAAGGGCAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGgcgc                                                                                                                       <  1:5016774/150‑1 (MQ=255)
                               aTCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGgcgc                                                                                                                       >  2:5442192/1‑150 (MQ=255)
                                tCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGgcgcg                                                                                                                      >  2:3690713/1‑150 (MQ=255)
                                  aCGGCGTGGTTCGCGATAAGGGAAAGTGCGGCGACGGAGGCCTGAACCTGCTGCTGCAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCAACGCTCGAGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCt                                                                                                                    >  2:853144/1‑150 (MQ=255)
                                      cGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAgg                                                                                                                <  2:2809398/150‑1 (MQ=255)
                                       gTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGc                                                                                                               >  1:1032642/1‑150 (MQ=255)
                                           ttCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATc                                                                                                           <  2:2192623/150‑1 (MQ=255)
                                              gcgAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGcc                                                                                                        <  2:2081499/150‑1 (MQ=255)
                                               cgAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCt                                                                                                       >  2:4330761/1‑150 (MQ=255)
                                                gagaAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCtt                                                                                                      <  2:1508551/150‑1 (MQ=255)
                                                gagaAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCtt                                                                                                      <  2:4424890/150‑1 (MQ=255)
                                                  gaAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGa                                                                                                    <  1:5205937/150‑1 (MQ=255)
                                                    aGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGt                                                                                                  >  2:1892126/1‑150 (MQ=255)
                                                    aGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGt                                                                                                  >  2:388036/1‑150 (MQ=255)
                                                      ggAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTgg                                                                                                <  1:148274/150‑1 (MQ=255)
                                                       gAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGa                                                                                                    <  2:3161191/145‑1 (MQ=255)
                                                       gAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGa                                                                                                    >  1:3161191/1‑145 (MQ=255)
                                                            tGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCgg                                                                                          >  2:980085/1‑150 (MQ=255)
                                                              cggcgACGGTTGCATGTACCTGCTGCTGTAACTGCTATAAGCCGATCAGCTGCGTCCGTGGCCACGCTCGTGGCTGCTTATACCCTGCCAGGCTTGGCTGACCGGGTCACGTCAGGCGAGATGATTGAATGGCATTGTGTTGCGCCGGGt                                                                                        <  2:4746355/150‑1 (MQ=255)
                                                               ggcgACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTg                                                                                       >  2:2414293/1‑150 (MQ=255)
                                                                  gACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCGTAGACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCt                                                                                    <  2:4590693/150‑1 (MQ=255)
                                                                   aCGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTa                                                                                   <  1:2013111/150‑1 (MQ=255)
                                                                            tGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCacgcaac                                                                          >  1:2596974/1‑150 (MQ=255)
                                                                                ccTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCAt                                                                      <  1:4330718/150‑1 (MQ=255)
                                                                                 ctgctgctgGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATg                                                                     >  1:2097303/1‑150 (MQ=255)
                                                                                 ctgctgctgGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATg                                                                     >  1:4069222/1‑150 (MQ=255)
                                                                                  tgctgctgGAACTGCAATACGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGTCGCGCTCAGGCAATCGCCTTGAGTGGAGCCGGGTGGCTATCACGCAACGCATGc                                                                    <  1:3685346/150‑1 (MQ=255)
                                                                                  tgctgctgGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGc                                                                    <  2:4000400/150‑1 (MQ=255)
                                                                                       ctgGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAgctg                                                               <  1:1631817/150‑1 (MQ=255)
                                                                                        tgGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAgctgc                                                              <  1:3043633/150‑1 (MQ=255)
                                                                                        tgGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAgctgc                                                              >  2:2040940/1‑150 (MQ=255)
                                                                                        tgGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAgctgc                                                              >  1:3585395/1‑150 (MQ=255)
                                                                                          gAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTc                                                            >  1:3591425/1‑150 (MQ=255)
                                                                                              tGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTccacc                                                        >  1:4206689/1‑150 (MQ=255)
                                                                                               gCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGTCTCGCTCAGGCAATCGCCTTGAGTGGCGCAGGTTGGCTATCTCGCTACGCATGCAgctgctgaagca                                                       <  2:5195564/150‑7 (MQ=255)
                                                                                               gCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTccacca                                                       >  1:4817999/1‑150 (MQ=255)
                                                                                                cAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAg                                                      >  2:2971756/1‑150 (MQ=255)
                                                                                                   tAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGtt                                                   <  1:3039104/150‑1 (MQ=255)
                                                                                                   tAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGtt                                                   <  1:4663232/150‑1 (MQ=255)
                                                                                                       ccGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGAcc                                               >  1:528065/1‑150 (MQ=255)
                                                                                                         gATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAg                                             <  1:1994274/150‑1 (MQ=255)
                                                                                                         gATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAg                                             >  1:5369019/1‑150 (MQ=255)
                                                                                                                ggcggcCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGa                                      >  1:2338735/1‑150 (MQ=255)
                                                                                                                 gcggcCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCTGGTCAACCTGGtggagca                                     <  2:911628/150‑7 (MQ=255)
                                                                                                                   cgcCGTGGCCACGCTCGTGGCTGCATATACCCCGCCACCCTGGGCTCACCGGCGCACCACACGCGCGCTCACGCAATCGCCTAGAGTGCCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGATGACCAGGTCACGACaa                                   <  2:154498/149‑1 (MQ=255)
                                                                                                                    gcCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAAc                                  <  2:5369019/150‑1 (MQ=255)
                                                                                                                       gTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGt                               >  1:2695770/1‑150 (MQ=255)
                                                                                                                       gTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGt                               >  1:3091522/1‑150 (MQ=255)
                                                                                                                            cACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAgg                          >  1:3891731/1‑150 (MQ=255)
                                                                                                                            cACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAgg                          >  1:4524389/1‑150 (MQ=255)
                                                                                                                                cTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCACCGGGTGGCTATCACGCAACGCATGCTGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCtt                      <  1:932818/150‑1 (MQ=255)
                                                                                                                                 tCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAgg                          <  1:2905037/145‑1 (MQ=255)
                                                                                                                                 tCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTg                     >  2:3477920/1‑150 (MQ=255)
                                                                                                                                 tCGTGGCTGCTTAGACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACGGTTGAGGTCTTg                     <  2:3591425/150‑1 (MQ=255)
                                                                                                                                     ggCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGa                            >  1:1606235/1‑139 (MQ=255)
                                                                                                                                           ttATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGAcc           >  2:4995165/1‑150 (MQ=255)
                                                                                                                                             ataCCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGAcccc         <  1:5004618/150‑1 (MQ=255)
                                                                                                                                                cccGGTCAGCCCAGGCTGGCCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGtt      >  1:4192044/1‑150 (MQ=255)
                                                                                                                                                cccGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGtt      >  2:3734701/1‑150 (MQ=255)
                                                                                                                                                 ccGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGAAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGTTg     >  1:2179739/1‑150 (MQ=255)
                                                                                                                                                 ccGGCCAGCCTGGGCTGACCGGGGCACGAAAGGCGCGCTCAGGCAATAGCCTTTAGTGGCGCCGGGTTGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTTTCAAGGTCACGACAAATGTTGATGTGTTGAAGCGCGACCCCGTTg     >  2:877154/1‑150 (MQ=255)
                                                                                                                                                  cGGTCAGCCCAGGCTGGCCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCtt                      >  2:4393718/1‑132 (MQ=255)
                                                                                                                                                  cGGTCAGCCCAGGCTGGCCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCtt                      <  1:4393678/132‑1 (MQ=255)
                                                                                                                                                  cGGTCAGCCCAGGCTGGCCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGTTGa    >  1:1411559/1‑150 (MQ=255)
                                                                                                                                                  cGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGTTGa    <  1:331520/150‑1 (MQ=255)
                                                                                                                                                  cGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGTTGa    <  1:3802351/150‑1 (MQ=255)
                                                                                                                                                   ggCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGTTGAg   <  2:4102627/150‑1 (MQ=255)
                                                                                                                                                    gCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGTTGAgg  <  1:5442192/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGTGTGGCCTTCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGTTGAGG  >  NC_002947/3855582‑3855879

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