Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,898,366 T→C 12.4% intergenic (+32/‑176) PP_4304 → / → sbp‑I VIC family cation transporter/sulfate ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,898,3660TC12.4% 559.9 / 48.6 214intergenic (+32/‑176)PP_4304/sbp‑IVIC family cation transporter/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (72/111);  new base C (19/7);  total (94/120)
Fisher's exact test for biased strand distribution p-value = 1.42e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 228 total aligned reads displayed.

CAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCA  >  NC_002947/4898217‑4898515
                                                                                                                                                     |                                                                                                                                                      
    gCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGc                                                                                                                                                    <  1:6053623/150‑1 (MQ=255)
    gCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGc                                                                                                                                                    <  1:6430261/150‑1 (MQ=255)
      cATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTc                                                                                                                                                  <  2:439109/150‑1 (MQ=255)
       aTGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAATTCCGCTcc                                                                                                                                                 >  2:5378305/1‑150 (MQ=255)
         gCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTa                                                                                                                                               <  2:4329814/150‑1 (MQ=255)
         gCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTa                                                                                                                                               <  1:2006504/150‑1 (MQ=255)
            tGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCATGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAg                                                                                                                                            <  1:6955220/150‑1 (MQ=255)
              gggAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGgc                                                                                                                                          <  1:152962/150‑1 (MQ=255)
                gAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGgcgc                                                                                                                                        <  2:2052078/150‑1 (MQ=255)
                gAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGgcgc                                                                                                                                        <  2:1165948/150‑1 (MQ=255)
                   cagcTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGa                                                                                                                                     <  1:2737094/150‑1 (MQ=255)
                   cagcTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGCTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGa                                                                                                                                     <  2:7088160/150‑1 (MQ=255)
                            cttgttcttgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGcgcg                                                                                                                            >  1:7051954‑M1/101‑150 (MQ=255)
                            caTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGcgcg                                                                                                                            <  1:2514590/150‑1 (MQ=255)
                            caTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGcgcg                                                                                                                            <  1:4039552/150‑1 (MQ=255)
                              tGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGt                                                                                                                          >  1:5269473/1‑150 (MQ=255)
                               gATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTc                                                                                                                         <  1:7599403/150‑1 (MQ=255)
                               gATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTc                                                                                                                         <  1:8501538/150‑1 (MQ=255)
                                 ttGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTccc                                                                                                                       <  1:6412083/150‑1 (MQ=255)
                                  tGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAGATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcccc                                                                                                                      <  1:8552856/150‑1 (MQ=255)
                                  tGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcccc                                                                                                                      <  1:8080592/150‑1 (MQ=255)
                                   gtgccgggtgttggggcagagtgccagcgactatcccgggggggttgatacgggttacagaacaaaaaaggctttggttatgcgtctgcctgtGCCGGCCTCGTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTccccc                                                                                                                     >  2:5415896‑M1/94‑150 (MQ=255)
                                   gCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTccccc                                                                                                                     <  1:6070490/150‑1 (MQ=255)
                                   gCCCGACCTCCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTccccc                                                                                                                     <  2:8471860/150‑1 (MQ=255)
                                    cccGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCg                                                                                                                    <  1:7647006/150‑1 (MQ=255)
                                    agaagggtgttggggcagattgccagccactataacgggtgggattatacgggataaagaacaaaaaatgcattggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCg                                                                                                                    >  2:6785243‑M1/93‑150 (MQ=255)
                                      aagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTa                                                                                                                  >  2:5122737‑M1/91‑150 (MQ=255)
                                       gACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAg                                                                                                                 <  1:6710939/150‑1 (MQ=255)
                                       gACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAg                                                                                                                 <  2:4572919/150‑1 (MQ=255)
                                       agggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAg                                                                                                                 >  1:3306552‑M1/90‑150 (MQ=255)
                                         ccTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGa                                                                                                               <  1:7386653/150‑1 (MQ=255)
                                           tgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGc                                                                                                             >  1:4769733‑M1/86‑150 (MQ=255)
                                           tgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGc                                                                                                             >  1:3244322‑M1/86‑150 (MQ=255)
                                           tGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGc                                                                                                             <  1:5575881/150‑1 (MQ=255)
                                             ttggggcagatttccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCgg                                                                                                           >  2:7524766‑M1/84‑150 (MQ=255)
                                             ttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCgg                                                                                                           >  1:1672292‑M1/84‑150 (MQ=255)
                                             ttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCgg                                                                                                           >  1:7319431‑M1/84‑150 (MQ=255)
                                             ttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCgg                                                                                                           >  1:8848469‑M1/84‑150 (MQ=255)
                                             ttcgggcagattgccagcgactataacgggtgggattatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCgg                                                                                                           >  1:8669289‑M1/84‑150 (MQ=255)
                                             ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCgg                                                                                                           <  2:1431021/150‑1 (MQ=255)
                                             ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCgg                                                                                                           <  2:4620478/150‑1 (MQ=255)
                                               ggggcagattgccagcgactataacgggtgggtttatatgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGt                                                                                                         >  2:6595880‑M1/82‑150 (MQ=255)
                                               ggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGt                                                                                                         >  1:4481539‑M1/82‑150 (MQ=255)
                                               ggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGt                                                                                                         >  2:5173019‑M1/82‑150 (MQ=255)
                                               acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGt                                                                                                         >  2:2872562‑M2/1‑149 (MQ=255)
                                                caAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtt                                                                                                        >  2:5863248‑M2/1‑148 (MQ=255)
                                                caAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCGTAGCGGGTTAAAACCCCCCGACAGGCGCGGCCCCGGCGCGGGCCCCCCGAGGGGCggggt                                                                                                        >  1:1955200/1‑148 (MQ=255)
                                                 aagaagACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAACGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGttt                                                                                                       >  2:8774886‑M2/1‑147 (MQ=255)
                                                 aagaagACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGttt                                                                                                       >  2:4169388‑M2/1‑147 (MQ=255)
                                                       ttgccagcgactagaacgggtgggtttatacgggttaaagaacaaaaaatgctgtggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCcgc                                                                                                >  2:8368523‑M1/74‑150 (MQ=255)
                                                         ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGGCCCCCGTAGGAGCGGGtttacccgcgc                                                                                               <  2:4918805‑M2/150‑12 (MQ=255)
                                                          ccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcccgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGaa                                                                                             >  2:3673762‑M1/71‑150 (MQ=255)
                                                           aCGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtttacccgcgaag                                                                                             <  1:8297402‑M2/150‑14 (MQ=255)
                                                                cAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCAAATGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGTTAAATCCGCTCCTACAGGCGCGGACGGTGCGCGGTCCCCCGTTGGAGCGGGtttacccgcgaagaggcc                                                                                        >  1:6643334‑M2/1‑132 (MQ=255)
                                                                 aGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtttacccgcgaagaggccg                                                                                       >  2:8684876‑M2/1‑131 (MQ=255)
                                                                 aGGGGGCGGCGTTCGGTTCTCGCTGTGGGATTGGGTTGTTTCCCAAAGGGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACCGGCGCGGACCGTGCGCGGTCCCCCGTTGGAGCGGAtttacccgcgaagaggccg                                                                                       >  2:5461533‑M2/1‑131 (MQ=255)
                                                                       cggcgTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAg                                                                                >  2:4412458/1‑150 (MQ=255)
                                                                           gttctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAAACGCTACAACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAg                                                                            >  2:3406486/1‑150 (MQ=255)
                                                                               tgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTATATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGc                                                                        <  2:1323421/150‑1 (MQ=255)
                                                                                        tgtgGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGTTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGTAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCaaa                                                               <  2:1220062/150‑1 (MQ=255)
                                                                                         ttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTa                                                                                                      <  2:1113629‑M1/72‑1 (MQ=255)
                                                                                             gAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCACTTCGCGGGTAAATCCGCTCCTACAGGCCCGGACCCTGCGCGGTCCCCCCTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCACAATCAAAACCAAAACAtt                                                          >  1:7567004/1‑150 (MQ=38)
                                                                                               aTGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACACCCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAtttt                                                        <  2:3582023/150‑1 (MQ=255)
                                                                                                tGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAttttt                                                       <  2:5983201/150‑1 (MQ=255)
                                                                                                 gggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAtttttt                                                      <  2:2208148/150‑1 (MQ=255)
                                                                                                    ttgttTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCCCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACACGCAGAAGCAAAACCAAAGCATTTTTTGtt                                                   >  2:7227869/1‑150 (MQ=255)
                                                                                                          ccaaaaGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTCCAGGCGCGGACCGTGCGCGGTCCCCCTTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTaa                                             <  2:3113048/147‑1 (MQ=255)
                                                                                                           ccAAAGCGTGTTAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAc                                            <  2:1365482/150‑1 (MQ=255)
                                                                                                           ccAAAGCGTGATAGCGCCTGCGCCGGCCTCGTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCCTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCCGAAGCATAACCCCAGCATTTTTTGGTCTTTAAc                                            >  2:6737339/1‑150 (MQ=38)
                                                                                                           ccAAAGCGTGATAGCGCCTGCCCCGGCCTCTACGCGGGAAAAACCGCTCCTACAGGCGCGCACCGTGCGCGGACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAACCATTTTTTGTTCTTTAAc                                            <  2:5731068/150‑1 (MQ=38)
                                                                                                             aaaGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCACCTACCGGCGCGGGCCGTGCGCGGTCCCCCGGCGGAGCGGGTTTACCCGCGCAGAGGCCGGGCCAGGCAGGAGCAAAACCAAAGCAATTTTTTTTGTTTAAccc                                          >  2:1642434/1‑150 (MQ=255)
                                                                                                                    tatcgcgcGTGCGCGGGCCTCTTGGGGTTTAATCCCGCTCATACCGGGGGGGACGGGGCGGGGCCCCCGTTCGGCGGGGTTTTACCCGCTATGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATaaa                                   <  1:2433095/146‑1 (MQ=255)
                                                                                                                    gctcgcgcCTGCGCCGGCCTCTGCGCGGGCAAAGCCGCTCCGACCGGCGCGGCCCGTGGGCGGTCCCCCGTTGGAGCGGGTTTACCCGGGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATaaa                                   <  1:1316332/146‑1 (MQ=255)
                                                                                                                    gATCGCGCCTGCGCCGGCCTCTTCGCGGGTAACTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATaaa                                   <  2:2692729/150‑1 (MQ=255)
                                                                                                                      tAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAacc                                 <  2:2810737/150‑1 (MQ=255)
                                                                                                                          tgcTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGt                                                                                                                                 >  2:6317724/3‑51 (MQ=255)
                                                                                                                               cgcCGGCCTCTTCGCGGTTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTta                        <  2:3142900/150‑1 (MQ=255)
                                                                                                                                 cgGGCCTCTGCGCGGGTAAACCCGCCCCTCCGGGGGGGGCCCGTGCGCGGCCCCCCGTAGGCGGGGTTTACCCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTtata                      >  1:5378305‑M1/3‑150 (MQ=255)
                                                                                                                                 ccGCCCTCTGCGGGGGTAACCCCGCACCCCCCGGCGCGGCCGGGGCGGGGGCCCCCGTCGGTGGGGGTTAACCGGCGAAGCGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTtata                      <  1:1148078/150‑1 (MQ=255)
                                                                                                                                 CCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAATCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTtata                      <  2:7021173‑M1/150‑1 (MQ=255)
                                                                                                                                 CCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGGGCGCGGTCCCCCGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTtata                      >  1:8186882‑M1/1‑150 (MQ=255)
                                                                                                                                    GCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTc                   <  2:7278532‑M1/150‑1 (MQ=255)
                                                                                                                                      CTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGc                 <  2:8548043‑M1/150‑1 (MQ=255)
                                                                                                                                      CTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGc                 <  2:3759718‑M1/150‑1 (MQ=255)
                                                                                                                                             CGGGTAAATCCGCTCCTACCGGCGCGGACCGTGCGCGGTCCCCCGTCGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAACGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAAt          >  2:6252224‑M1/1‑150 (MQ=255)
                                                                                                                                              GGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAac                                  <  1:3819301‑M1/125‑1 (MQ=255)
                                                                                                                                               GGTAAATCCGCTCCTACAGGCGCGGACCGGGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCt        >  1:6738226‑M1/1‑150 (MQ=255)
                                                                                                                                                gTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccaca                                             >  1:845972‑M2/1‑52 (MQ=255)
                                                                                                                                                   actCCGCTCCTCCCGGCGGGGACGGTGCGCGGTCCCCCGTTGGAGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGccc    <  1:5749061/148‑1 (MQ=255)
                                                                                                                                                   aaTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttg     <  1:1045947‑M2/150‑102 (MQ=255)
                                                                                                                                                    aTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt    <  1:7477609‑M2/150‑103 (MQ=255)
                                                                                                                                                    aTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt    <  1:6130354‑M2/150‑103 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtg   <  2:7667089‑M2/150‑104 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtg   <  1:3677159‑M2/150‑104 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgggggtcttctgggg   <  2:7755072‑M2/150‑104 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccccccctgctcgtgggtcttcttgtgg   <  1:6623760‑M2/150‑104 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  1:1115818/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  2:823673/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa  >  2:398423/1‑150 (MQ=255)
                                                                                                                                                     tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGGCGCTGGCAATCTGCCCCa  >  2:8877750/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCA  >  NC_002947/4898217‑4898515

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: