Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,970,399 T→G 8.8% T51P (ACG→CCG)  PP_5723 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,970,3990TG 8.8% 538.3 / 13.0 171T51P (ACG→CCG) PP_5723hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/120);  new base G (14/4);  total (47/124)
Fisher's exact test for biased strand distribution p-value = 3.57e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

AAACAGCGCGTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCA  >  NC_002947/5970250‑5970548
                                                                                                                                                     |                                                                                                                                                      
 acaggcgcgaggcctggtcgctcagggcatcgagaagggctttgggcgcgagcatggggctgtccactgaaacggagaaatggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTc                                                                                                                                                       >  2:1371249‑M1/138‑150 (MQ=255)
 aaCAGCGCGTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTc                                                                                                                                                       <  1:848313/150‑1 (MQ=255)
 aaCAGCGCGTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTc                                                                                                                                                       <  1:8535595/150‑1 (MQ=255)
  caggcgcgaggcctggtcgctcagggcatcgagaagggctttgggcgcgagcatggggctgtccactgaaacggagaaatggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTcg                                                                                                                                                      >  2:5965310‑M1/137‑150 (MQ=255)
  caggcgcgaggcctggtcgctcagggcatcgagaagggctttgggcgcgagcatggggctgtccactgaaacggagaaatggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTcg                                                                                                                                                      >  2:3328077‑M1/137‑150 (MQ=255)
  aCAGCGCGTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTcg                                                                                                                                                      <  1:8412483/150‑1 (MQ=255)
  aCAGCGCGTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTcg                                                                                                                                                      <  2:2282100/150‑1 (MQ=255)
   cAGCGCGTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTcgc                                                                                                                                                     <  1:2031184/150‑1 (MQ=255)
   cAGCGCGTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTcgc                                                                                                                                                     <  1:52248/150‑1 (MQ=255)
    aGCGCGTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGTGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTcgcg                                                                                                                                                    >  1:304914/1‑150 (MQ=255)
     gcgcgaggcctggtcgctcagggcatcgagaagggctttgggcgcgagcatggggctgtccactgaaacggagaaatggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTCGCGa                                                                                                                                                   >  2:6220231‑M1/134‑150 (MQ=255)
     gcgcgaggcctggtcgctcagggcatcgagaagggctttgggcgcgagcatggggctgtccactgaaacggagaaatggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTCGCGa                                                                                                                                                   >  2:1526114‑M1/134‑150 (MQ=255)
      cgcgaggcctggtcgctcagggcatcgagaagggctttgggcgcgagcatggggctgtccactgaaacggagaaatggccgccagtgtagcacgcacacacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTCGCGAt                                                                                                                                                  >  2:75275‑M1/133‑150 (MQ=255)
       gcgaggcctggtcgctcagggcatcgagaagggctttgggcgcgagcatggggctgtccactgaaacggagaaatggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTCGCGAtt                                                                                                                                                 >  1:1039666‑M1/132‑150 (MQ=255)
        cgtgaggtgAGTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTg                                                                                                                                                <  1:6810985/144‑1 (MQ=255)
        cgTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTg                                                                                                                                                <  2:1678556/150‑1 (MQ=255)
        cgTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTg                                                                                                                                                <  1:8622592/150‑1 (MQ=255)
        cgTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGGcgcgggtg                                                                                                                                                >  2:2882914/1‑146 (MQ=255)
         gaggcctggtcgctcagggcatcgagaagggctttgggcgcgagcatggggctgtccactgaaacggagaaatggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTCGCGATTgg                                                                                                                                               >  1:6652592‑M1/130‑150 (MQ=255)
          ttATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTggg                                                                                                                                              <  2:7188352/150‑1 (MQ=255)
          ttATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTggg                                                                                                                                              <  2:8872083/150‑1 (MQ=255)
            aTGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGt                                                                                                                                            <  2:4583384/150‑1 (MQ=255)
             tgtgATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTgggggg                                                                                                                                           <  1:3677718/150‑3 (MQ=255)
              gtgATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGTGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGc                                                                                                                                          <  2:304914/150‑1 (MQ=255)
              gtgATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGc                                                                                                                                          <  1:7479683/150‑1 (MQ=255)
                  tttGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAg                                                                                                                                      >  2:5907925/1‑150 (MQ=255)
                  tcgctcagggcatcgagaagggctttgggcgcgagcatggggctgtccactgaaacggagaaatggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGGGTgggggttgggggggacgg                                                                                                                                      <  1:5900466‑M1/30‑8 (MQ=255)
                    gctcagggcatcgagaagggctttgggcgcgagcatggggctgtccactgaaacggagaaatggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTgcggcgcgggtggggtgggccggg                                                                                                                                    <  1:7477772‑M1/32‑8 (MQ=255)
                     gCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAgcgc                                                                                                                                   <  2:218329/150‑1 (MQ=255)
                      cTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGTGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGcc                                                                                                                                  <  1:871432/150‑1 (MQ=255)
                      cTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGcc                                                                                                                                  <  1:871309/150‑1 (MQ=255)
                         gggATTGGCCAGTCGCCAGCTCGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCATCAGCTTGTAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGATGCGcccac                                                                                                                               <  2:1770150/150‑3 (MQ=255)
                           gATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGGGGCGCGGGTGGGGGGGGGAGCGccccccc                                                                                                                             >  1:2241783/1‑149 (MQ=255)
                            aTTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCaa                                                                                                                            <  2:319385/150‑1 (MQ=255)
                             ttGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAg                                                                                                                           <  1:6697335/150‑1 (MQ=255)
                             ttGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAg                                                                                                                           <  1:6697280/150‑1 (MQ=255)
                                gCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTTGGGGGCGGGTTGGGGGGGAAGCGCCCCCCCAGaat                                                                                                                        >  2:4159963/1‑148 (MQ=255)
                                 ccAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGAttt                                                                                                                       <  1:4006419/150‑1 (MQ=255)
                                 ccAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGAttt                                                                                                                       <  2:1741253/150‑1 (MQ=255)
                                 ccAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCGTGCGTGACGATTTGGGTGCGCAGCGCCCCCCCAGaagt                                                                                                                       >  1:7754742/1‑147 (MQ=255)
                                  cAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGGCGGGGGTGGGGGGCGCGGCGCCCCCCCAgggttg                                                                                                                      >  2:6027490/1‑145 (MQ=255)
                                   aGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTtgt                                                                                                                     <  2:7679370/150‑1 (MQ=255)
                                     tCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGc                                                                                                                   <  1:5299278/150‑1 (MQ=255)
                                     tCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGc                                                                                                                   <  2:3835398/150‑1 (MQ=255)
                                       gCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCt                                                                                                                 <  1:1645256/150‑1 (MQ=255)
                                       gCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGACAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCt                                                                                                                 <  2:7777186/150‑1 (MQ=255)
                                        ccAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCGTGGGGGGGGGGTTGGGGGGGGGGCGCCCCCCAACAATtgtggggg                                                                                                                >  2:2502323/1‑146 (MQ=255)
                                           gCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGGGGCGCGGGAGGGGGGCGGAGCGCCCCCCCAGAGTtttggcgtttt                                                                                                             >  2:5751070/1‑143 (MQ=255)
                                            cACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATg                                                                                                            <  2:9043242/150‑1 (MQ=255)
                                             aCGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGc                                                                                                           <  1:820430/150‑1 (MQ=255)
                                               gggAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTg                                                                                                         <  2:5208693/150‑1 (MQ=255)
                                               gggAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTg                                                                                                         <  2:440458/150‑1 (MQ=255)
                                               gggAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTg                                                                                                         <  2:4878119/150‑1 (MQ=255)
                                                ggAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGa                                                                                                        <  2:3952319/150‑1 (MQ=255)
                                                 gAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGGGGCGGGGTTTTGGTGGGGAGCGCCCCCCACCAATTTTGGCATGATGCTggt                                                                                                       >  2:3462704/1‑148 (MQ=255)
                                                 gAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGGCGCGTTTGGGGGGCGAAGCGCCCCCCCACCTGTGTGCCGTGATGCTGat                                                                                                       >  2:7477772/1‑150 (MQ=255)
                                                       tGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGcc                                                                                                <  1:3234955/150‑1 (MQ=255)
                                                         cacaATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCgg                                                                                              <  1:8730044/150‑1 (MQ=255)
                                                         cacaATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCgg                                                                                              <  1:1395543/150‑1 (MQ=255)
                                                          acaATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCggg                                                                                             <  1:5641411/150‑1 (MQ=255)
                                                           caATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCgggg                                                                                            <  2:9115211/150‑1 (MQ=255)
                                                            aaTTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGc                                                                                           <  2:5439568/150‑1 (MQ=255)
                                                              tttCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCg                                                                                         <  1:1044604/150‑1 (MQ=255)
                                                                tCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCt                                                                                       <  2:420454/150‑1 (MQ=255)
                                                                tCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCt                                                                                       <  1:9019623/150‑1 (MQ=255)
                                                                 cGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCtg                                                                                      <  2:7527880/150‑1 (MQ=255)
                                                                 cGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCtg                                                                                      <  1:2372537/150‑1 (MQ=255)
                                                                  ggTGCGCATCAGCGTGGCAAGCTGCGCGAAAAATGTTCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCtgt                                                                                     <  1:5807611/150‑1 (MQ=255)
                                                                  ggTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCtgt                                                                                     <  1:8194290/150‑1 (MQ=255)
                                                                    tGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGGCGCGCGTGGGGGGCGCCGCGCCCCCCCACAGTTGTGCCGTGTTGGTGAGAGTTGCCGGGGcgcggggg                                                                                    >  2:5900466/1‑143 (MQ=255)
                                                                     gcgcAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCa                                                                                  <  1:8110614/150‑1 (MQ=255)
                                                                       gcagcaGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGc                                                                                <  1:9119125/150‑1 (MQ=255)
                                                                        cagcagCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGcc                                                                               <  1:5470298/150‑1 (MQ=255)
                                                                           gagaaatggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTCGCGATTGCGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCaa                                                                            >  2:1067361‑M1/64‑150 (MQ=255)
                                                                            agcGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAAt                                                                           <  2:5802639/150‑1 (MQ=255)
                                                                            agaaatggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAAt                                                                           >  2:1713339‑M1/63‑150 (MQ=255)
                                                                            agaaatggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAAt                                                                           >  2:6490043‑M1/63‑150 (MQ=255)
                                                                             gaaatggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATc                                                                          >  1:6102398‑M1/62‑150 (MQ=255)
                                                                               ggggCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATcgc                                                                        <  1:5324215/150‑1 (MQ=255)
                                                                                 ggCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGa                                                                      <  1:8996746/150‑1 (MQ=255)
                                                                                 ggCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGa                                                                      <  1:7851767/150‑1 (MQ=255)
                                                                                 ggCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGa                                                                      <  2:7185001/150‑1 (MQ=255)
                                                                                 ggCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGa                                                                      <  2:6273573/150‑1 (MQ=255)
                                                                                  ggccgccagtgtagcacgcacaaacctggcgcctgagtgcacattggcctttgtgtGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGAc                                                                     >  2:8440210‑M1/57‑150 (MQ=255)
                                                                                  gCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGAc                                                                     <  1:2003823/150‑1 (MQ=255)
                                                                                    aaGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGc                                                                   <  2:6066240/150‑1 (MQ=255)
                                                                                      gCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCaa                                                                 <  2:1187713/150‑1 (MQ=255)
                                                                                      gCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCaa                                                                 <  1:2940961/150‑1 (MQ=255)
                                                                                           gcgAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGGCGCGAGTGGGGTGCGCAGCCCCCCCCCACAATTTTTCCGTGATGCTGGTATTGCCCGGGCCCCGGGGCAGCCCAAAAGCGCCGCAACGCCg                                                            >  2:8957847/1‑150 (MQ=255)
                                                                                             gAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCGTTCGTGGCGGGTTGGGTGCGGAGCGCCCCCCCAGATTTGTTGCGTGTTGGTGATATTTGCGGGGGCGCGGTGGGGCCCCAAAGAGACGCCACGCCGca                                                          >  2:1422905/1‑150 (MQ=255)
                                                                                                        aaCAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGGGGGGCGGTTTGGGGGGGGAGCGCCCCCCAAGATTTTTGGCGTTATGGTGTTAATGCCGGGGGCGCGGGGCAGCCCACTAGCGACGCAACGCCGCACACACAGGccc                                               >  2:5199059/1‑149 (MQ=255)
                                                                                                              ttGAAAACTGCACGCGATCTTTGTGGGAGCGGCCGTGCGGCGCCCTTGGGGGGGGGATCCCCCCCCCCCCTTTGTGGCGGGATTGTGGTCTTTCCCGGGCCCCGGTGCCGCCCAACCCCCGCGCCACGCCCCCCCCAAAGGcccaagggg                                         >  2:6636783/1‑143 (MQ=255)
                                                                                                                                 tttGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCg                                                                                         <  2:6842107/83‑1 (MQ=255)
                                                                                                                                 tttGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCg                                                                                         >  1:6842107/1‑83 (MQ=255)
                                                                                                                                   tgtgGGCGGGGCCTTGGGGCGCGCTGGGGGTGCGACGCGCCCCCCACGCTTTGTGCCTGTTTGATTATCTGGCGGGGGCCGCTGTGCCGCCCAAGCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTgtgc                    <  2:8209894/150‑1 (MQ=255)
                                                                                                                                      gggAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCt                                                                                       >  1:6478969/1‑80 (MQ=255)
                                                                                                                                                TTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTggc       <  2:4683800‑M1/150‑1 (MQ=255)
                                                                                                                                                    GTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGcc   <  1:377231‑M1/150‑1 (MQ=255)
                                                                                                                                                    GTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGcc   <  1:8450004‑M1/150‑1 (MQ=255)
                                                                                                                                                     TCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAATGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCa  <  2:3569599‑M1/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AAACAGCGCGTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTCGGTGCGCAGCAGCGGGGCAAGCTGCGCGAAAAATGTGCAACAGCTTGAAAACTGCACGCGATCTTTGTGGGAGCGGCCTTGCGTCGCGATTGGGGTGCGAAGCGCCCCCCAAGATTTGTGCCTTGATGCTGATATTGCCGGGGCCGCTGTGCAGCCCAATCGCGACGCAAGGCCGCACACAAAGGCCAATGTGCACTCAGGCGCCAGGTTTGTGCGTGCTACACTGGCGGCCA  >  NC_002947/5970250‑5970548

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