Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,104,942 G→C 6.9% intergenic (‑61/‑144) PP_1875 ← / → PP_1876 two‑component system sensor histidine kinase/response regulator/AAA family ATP/GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,104,9420GC 6.9% 523.9 / 20.3 158intergenic (‑61/‑144)PP_1875/PP_1876two‑component system sensor histidine kinase/response regulator/AAA family ATP/GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (58/89);  new base C (4/7);  total (62/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 167 total aligned reads displayed.

AACCGGGCTTGCTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGAC  >  NC_002947/2104793‑2105091
                                                                                                                                                     |                                                                                                                                                        
aaCCGGGCTTGCTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCggg                                                                                                                                                          <  1:1818549/150‑1 (MQ=255)
 aCCGGGCTTGCTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCgggg                                                                                                                                                         >  1:3006915/1‑150 (MQ=255)
 aCCGGGCTTGCTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCgggg                                                                                                                                                         <  1:1038680/150‑1 (MQ=255)
      gcttgcttGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGc                                                                                                                                                    <  2:6112165/150‑1 (MQ=255)
        ttgcttgTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCtt                                                                                                                                                  <  2:3441324/150‑1 (MQ=255)
         tgcttgTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTg                                                                                                                                                 <  2:412581/150‑1 (MQ=255)
          gCTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGG                                                                                                                                                >  1:1095923‑M1/1‑150 (MQ=255)
           cTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGC                                                                                                                                               <  1:5386101‑M1/150‑1 (MQ=255)
           cTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGCCCCCGCTTGGC                                                                                                                                               <  1:5313562‑M1/150‑1 (MQ=255)
            ttGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCG                                                                                                                                              >  2:2715658‑M1/1‑150 (MQ=255)
            ttGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCG                                                                                                                                              >  1:6829524‑M1/1‑150 (MQ=255)
                 tGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCT                                                                                                                                         <  2:5680843‑M1/150‑1 (MQ=255)
                     cggCGTTTGGTTAATTTGATATCCTTGAGTCCCGTGCGGGTTTCCCGTTCCTTGCGCTGCGCAAGCTTAGCCGATTAGCCCCGGGTGGTGGTGGTTCTGTCTGTTAATTTCTGTTGATCGCCAAGGGGGGGCCGCTTGGCGCCCCTCcgc                                                                                                                                     <  1:3203528/149‑1 (MQ=255)
                     aGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATGGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGC                                                                                                                                     >  2:5526526‑M1/1‑150 (MQ=255)
                       gCGATTGGTGAACTTGATAGCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATCGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTGGATCGCCAAGGGGGGGCCGCTTGGCGCCCCTCCGCGG                                                                                                                                   <  2:1697886‑M1/150‑1 (MQ=255)
                         gATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGC                                                                                                                                 <  2:4599977‑M1/150‑1 (MQ=255)
                         gATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGC                                                                                                                                 <  2:2604116‑M1/150‑1 (MQ=255)
                         gATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCACCGCGGGC                                                                                                                                 >  2:6663319‑M1/1‑150 (MQ=255)
                         gATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCGGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGC                                                                                                                                 <  2:4438062‑M1/150‑1 (MQ=255)
                         gATTGGTGAACTTGATATCCATGAGGCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGGGGGGGCCGCTTGGCGCCCCTCCGCGGGC                                                                                                                                 <  1:3846748‑M1/150‑1 (MQ=255)
                           ttGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTT                                                                                                                               <  1:2161737‑M1/150‑1 (MQ=255)
                             ggTGACCTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGC                                                                                                                             <  2:1163001‑M1/150‑1 (MQ=255)
                              gTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCC                                                                                                                            >  1:1263806‑M1/1‑150 (MQ=255)
                               tGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGGGTGCCC                                                                                                                           >  2:6484194‑M1/1‑150 (MQ=255)
                                   cTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTc                                                                                                                       <  2:2721867/150‑1 (MQ=255)
                                    ttGATATCCATGAGTCCCGAGCCGGGTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGACTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCAAGCCCGcgga                                                                                                                      <  2:5036644‑M1/150‑5 (MQ=255)
                                     ttatatccATGAGTCCCGAGCGGGGTTCCGGTTCCTTGCGCTGCGCACGCATCGCCGAATCGCCCGGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCACGCGGGGGCGGCTTGGCGCCCCTCCGCGGGCAAGCCCGcggag                                                                                                                     <  1:24855‑M1/148‑6 (MQ=255)
                                     tGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTccc                                                                                                                     <  2:6239756/150‑1 (MQ=255)
                                      gATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctccca                                                                                                                    <  2:4774439‑M1/150‑7 (MQ=255)
                                      gATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCca                                                                                                                    <  2:1177191/150‑1 (MQ=255)
                                        tatCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTTCCCCTTACcaca                                                                                                                  >  2:685850/1‑150 (MQ=255)
                                        tatCCAGGAGTCCCGAGCCGGTTTCCCGGTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCcaca                                                                                                                  <  1:6227833/150‑1 (MQ=255)
                                         ctCCATGAGTCCCGAGCCGGTGTCCCCCTCCTCGCGCAGCGCAAGCAGAGCCGCATAGCCCCGGCCGCTGGTGGGTCTGTCTGTTCATTACAGAGGATCGCCAAGCGGCGGCCGCTGGGCGCCCCTCCGCGGGCTTGCCCGTCCCcaccg                                                                                                                 <  2:1974659/149‑3 (MQ=255)
                                           ccATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCGTGCCCGTTCCCACAGGt                                                                                                               >  2:6364258/1‑150 (MQ=255)
                                            cATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCGTGGCGCCCCTCCGCGGGCGGGCCCCGTCCcccacgga                                                                                                              >  2:1863141/1‑146 (MQ=255)
                                             aTGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAggggc                                                                                                             >  2:1084078/1‑147 (MQ=255)
                                             aTGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGCCCCCGCTTGGCGatc                                                                                                                                           >  1:5476614‑M1/1‑117 (MQ=255)
                                             aTGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGCCCCCGCTTGGCGatc                                                                                                                                           <  2:5476614‑M1/120‑4 (MQ=255)
                                                 gTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCt                                                                                                         <  2:1511500/150‑1 (MQ=255)
                                                 gTCCCGAGCCGGTTTCCCGTTCCTTGCGCGGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCt                                                                                                         <  2:6115396/150‑1 (MQ=255)
                                                     cGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttca                                                                                                     <  2:4822646‑M1/150‑22 (MQ=255)
                                                     cGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGc                                                                                                     <  2:1555107/150‑1 (MQ=255)
                                                        gcgggTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCTAATCGCCCCGGGTGCTGGGGGGTCTGTCTGTTCTTTTCTGTTGATCGCCACGGGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGcccc                                                                                                  <  1:195213/147‑1 (MQ=38)
                                                         ccGGTTTCCCGTTCCTGGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCa                                                                                                 <  2:905133/150‑1 (MQ=255)
                                                         ccGGTTGCCCGTTCCTTGCGCTGCGCAAGCAGAGCCGAATAGCCCCGGCTGCTGGTGGTTCTGTCTGTTCTTTTCTGTTGATCGCCATGGGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCa                                                                                                 <  1:2473246/150‑1 (MQ=39)
                                                                 ccGTTCCTTGCGCTGCGCAAGCATCGCCGAATAGCCCCGGCTGCTGGTGGGTCTTTCTGTTCATTTCAGTTGTTCGCCAAGGGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAAc                                                                                        <  1:3082699/150‑1 (MQ=255)
                                                                 ccGTTCCTTGCGCGGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTATTTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAAc                                                                                        <  1:3488236/150‑1 (MQ=255)
                                                                 ccGGTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAAc                                                                                        <  1:5078031/150‑1 (MQ=255)
                                                                        ttGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCATTCCCGGTCCCACAGGTACAGCTTAGCCCCAGGGAGAACGCAAGcc                                                                                 >  2:3975254/1‑150 (MQ=255)
                                                                         tgcgctgcgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacc                                                                                 <  1:2686837‑M1/150‑42 (MQ=255)
                                                                           cgctgcgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCGTGCCCCTTCCCCCCGGTAAAGCGTAGCCCCCGGCGGAAACCAAGcccgg                                                                              >  2:890868/1‑148 (MQ=255)
                                                                             ctgcgcAAGCTTAGCCGAATCGCCCCGGCTGGTGGTGGGTCTGTCTTTTAATTTCAGTTGATCGCCACGGGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCtgtg                                                                           <  1:590780/150‑1 (MQ=39)
                                                                             ctgcgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGTCCtgtg                                                                           <  2:685077/150‑1 (MQ=255)
                                                                             ctgcgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCtgtg                                                                           <  2:4819672/150‑1 (MQ=255)
                                                                             ctgcgcAAGCAGAGCCGAATAGCCCCGGTTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCtgtg                                                                           <  2:6017735/150‑1 (MQ=255)
                                                                              tgcgcATTCATAGCCTTATAGCCCCTGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTgg                                                                          <  2:3453565/150‑1 (MQ=39)
                                                                              ggcgcaAGCATAGCCGAATCGCCCCGGGTGCTGGTGGGTCTTTCTGTTCTTTTCATTTGATCGCCAAGGGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgg                                                                            <  1:2400927‑M1/149‑47 (MQ=255)
                                                                                cgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGa                                                                        <  1:4689421/150‑1 (MQ=255)
                                                                                 gcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaa                                                                         <  2:144849‑M1/150‑50 (MQ=255)
                                                                                  cAAGCATAGCCGAATAGCCCCGGTTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaac                                                                        <  2:3376501‑M1/150‑51 (MQ=255)
                                                                                  cAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTTTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGttcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaac                                                                        <  1:6034371‑M1/150‑51 (MQ=255)
                                                                                  cAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGc                                                                      >  2:6635146‑M2/1‑149 (MQ=255)
                                                                                   aaGCATCGCCGAATCGCCCCGGCTGGTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCACGCGGGGGCCGCTGGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcg                                                                     >  1:2558036‑M2/1‑148 (MQ=255)
                                                                                   aaGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTTTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcg                                                                     >  2:2379018‑M2/1‑148 (MQ=255)
                                                                                   aaGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacg                                                                       <  2:4536919‑M1/150‑52 (MQ=255)
                                                                                    aGCATCCCCGAATCGCCCCGGCTGCTGGGGGGTGTGTCTTTTAATTTCAGTTGATCGCCAAGGGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgg                                                                    >  1:2449023‑M2/1‑147 (MQ=255)
                                                                                    aGCATAGCCGAATCGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgg                                                                      <  1:33738‑M1/150‑53 (MQ=255)
                                                                                    aGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgg                                                                      <  1:1541996‑M1/150‑53 (MQ=255)
                                                                                    aGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgg                                                                      <  2:2813326‑M1/150‑53 (MQ=255)
                                                                                    aGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgg                                                                    >  2:3793799‑M2/1‑147 (MQ=255)
                                                                                    aGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgg                                                                    >  2:3753366‑M2/1‑147 (MQ=255)
                                                                                    aGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgg                                                                      <  2:5371278‑M1/150‑53 (MQ=255)
                                                                                    aGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgg                                                                    >  1:4353967‑M2/1‑147 (MQ=255)
                                                                                    aGCATAGCCGAATAGCCCCGGCGGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgg                                                                      <  2:4060855‑M1/150‑53 (MQ=255)
                                                                                     gcctAGCCGACTCGCCCCGGCTGCTGGTGGGTCTGTCTGGTCATGTCAGTGGATCGCCACGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCggg                                                                   <  2:4986641/147‑1 (MQ=39)
                                                                                     gCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGttcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacggg                                                                     <  1:5415900‑M1/150‑54 (MQ=255)
                                                                                     gCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcggg                                                                   >  2:3173349‑M2/1‑146 (MQ=255)
                                                                                      cATAGCCGAATAGCCCGGGCTGCTGGGGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggc                                                                  >  1:1767058‑M2/1‑145 (MQ=255)
                                                                                        tAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGTGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaa                                                                >  2:6583477‑M2/1‑143 (MQ=255)
                                                                                           ccGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCC                                                                                                                           >  1:6444618‑M1/1‑90 (MQ=255)
                                                                                           ccGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCC                                                                                                                           <  2:6444618‑M1/90‑1 (MQ=255)
                                                                                             gAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaagcccg                                                           >  2:4325042‑M2/1‑138 (MQ=255)
                                                                                               atcgCCCGGGCTGCTGGTGGTTCTGTCTTTTCTTTTCTGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTGGCCCGTTCCCACAGTTACAGCTTAGCCCCAGGCTTAACGCAAGCCCTGTGGGAGCGGGCAAGCCcgcg                                                         <  1:5116532/147‑1 (MQ=255)
                                                                                                   ccccGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCGTGCCCCGTGCCCCAGGGTCAGCGGATCCCCCGGCGGGAACCCAACCCAGGGGGGGCGGGCCAGCCCCCGgggg                                                     >  1:2240073/1‑150 (MQ=255)
                                                                                                          tgctgGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCACCGCGGGCGTGCCCCGTGCCCCAGCGACACCCTAGCCCCAACCAGAACGCAACCCCCGGGGGGGCCGGGAAGCCCCCGGGGGGgcgccc                                              >  1:1669217/1‑148 (MQ=255)
                                                                                                                 gggggtgtgtgtTCTTTTCTGTTGATGGCCATGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTGAGCCCCAGGCTGAACGCACGCCCTGTGGGAGCGGGCACGCCCGCGGAGGGGCGCGCAGCGGcc                                       <  2:5044854/145‑1 (MQ=38)
                                                                                                                           gtTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCCTTACCCCACGTAAAGCATAGCCCCACGCTGAAAGCAACCCCGGGGGGAGCCGGCCACCCCCCCGAGGGGCGCGCACCCGCCCCCAAAAAAc                            >  2:4339759/1‑150 (MQ=255)
                                                                                                                           gtTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCACCGCGGGCGGGCCCCCTGCCCCAAGTACACCGTAGCCCCACGCGGCACCCCACCCCAGGGGGGGCCGGCCCCCCCCCCGGGGGGCGCCCCGCGGCCCCCCCAAAAc                            >  2:3431704/1‑150 (MQ=255)
                                                                                                                            ttCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTa                                                                                                              <  1:4380788/70‑1 (MQ=255)
                                                                                                                              cttttCTGTTGAGGGCCAGGGGGGGGCCGGTGGCCCCCCCGCGGGGGGGTGGCCGGTACCCCCCGGTCCCGCTTCCCCCCCGGCTGACCGAACGCCGTGTGGGAGGGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTga                         <  2:6544191/148‑1 (MQ=255)
                                                                                                                              cttttCCGTTGAGCCCCACGGGGGGGCCGCTTGGCGCCCCCCCGCGGGCGGGCCCGTTCCCACAGGAACCGCTTCGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCACGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTga                         <  2:5582282/148‑1 (MQ=255)
                                                                                                                                    aGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTgg                                                                          <  2:4053379/96‑1 (MQ=255)
                                                                                                                                    aGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGC                                                                                                                             <  1:2726464‑M1/47‑1 (MQ=255)
                                                                                                                                     gTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCACCGCGGGCGTTCCCCGTACCACAGGTACCCCGTATCCCCACCCAGAAACCAATCCCCGTGGGCGCGGGCAACCCCCCCGGCGGGCGCGCACCGGCCCCCAAAAAAAAgagaatttt                  >  2:6597813/1‑145 (MQ=255)
                                                                                                                                          gcGCCAAGCGGCCCCCGCTTGGCGCCCCTCCGCGGGCGTGCCCGTTCCCACCGGTACAGCTTAGCCCCAGGCGGAAAGCAAGCCCGGGGGGGGCGGGCAAGCCCGCCGAGGGGCGCGCCGCGGCCCCCAAAAAACAGAGATTGTTAGGAg             >  2:2157650/2‑150 (MQ=255)
                                                                                                                                            gCCAAGCGGCCCCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGaa                                                                                         >  1:5735871/1‑74 (MQ=37)
                                                                                                                                             ccAAGCGGCCCCCGCTTGGCGCCCCTCCGCGGGCGTGCCCCTTGCCACAGGGACAGCGGAGCCCCACGCGGGACGCCACCCCCGGGGGGGCGGGCCCCCCCCCGGCGGGGCGCCCACCGGCCCCCCAAAAACAGAGCGTGTTAGGAgtag          >  2:5232382/1‑148 (MQ=255)
                                                                                                                                              cAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCt                                                                                            >  1:1565985/1‑69 (MQ=255)
                                                                                                                                              cAAGCGGGGGCCGCTTGGCGCCCCGCCCCGGGCGGGGCCCGTCCCCCACGGAAACCAGATCCCCACGCGGGACCCCAACCCCGTGGGGGCGGGGCAGCCCCCGGACGGGGGCGCCGCCGCCCCCCAAAAACAGGCAGTGTTTGGGgtcgg         >  2:3011044/1‑149 (MQ=255)
                                                                                                                                                   gggggCCGCTTGGCGCCCCTCCGCGGGCTTGCCGGTTCCCCCAGGTCCAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAAACTGAGATTGTTAGGAGTCGTCAt      <  2:6825800/150‑1 (MQ=38)
                                                                                                                                                     |                                                                                                                                                        
AACCGGGCTTGCTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGAC  >  NC_002947/2104793‑2105091

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