Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,155 T→A 18.0% intergenic (+106/‑207) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,1550TA18.0% 363.0 / 78.1 168intergenic (+106/‑207)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (54/80);  new base A (10/20);  total (68/100)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-02

Only 100 of 170 total aligned reads displayed.

CAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGC  >  NC_002947/955007‑955302
                                                                                                                                                     |                                                                                                                                                    
 cgCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                    <  1:4646991/149‑1 (MQ=255)
 aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                    <  2:2831330/150‑1 (MQ=255)
 aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                    <  2:7777013/150‑1 (MQ=255)
 aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                    <  2:7597133/150‑1 (MQ=255)
   ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                                  <  2:1221622/150‑1 (MQ=255)
   ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                                  <  1:4706483/150‑1 (MQ=255)
    cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                                  <  1:1178731/150‑1 (MQ=255)
    cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTg                                                                                                                                                 <  2:8401238/150‑1 (MQ=255)
     aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                                                                                                <  1:6944581/150‑1 (MQ=255)
     aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                                                                                                <  1:6944607/150‑1 (MQ=255)
      gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTa                                                                                                                                               <  1:5885239/150‑1 (MQ=255)
       aCCGCTGCTGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCAGGCCACTGTGCCGGGGTTTTTTTTTTGGGGTATTTTTCCTGTGCCGGCCTATTCTCGTGTAAACCCTCTCCTACAAGGTTCCTGTAg                                                                                                                                              <  1:2347541/150‑1 (MQ=255)
       aCCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTTCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAg                                                                                                                                              <  1:8727079/150‑1 (MQ=255)
       aCCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAg                                                                                                                                              <  1:2347217/150‑1 (MQ=255)
         cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAg                                                                                                                                           <  2:5345053/150‑1 (MQ=25)
          gCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAg                                                                                                                                           <  2:2469325/150‑1 (MQ=38)
           cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                          <  2:2219955/150‑1 (MQ=38)
           cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                                                                                                          <  1:8978601/150‑1 (MQ=255)
           cTGCGGATGAGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                                                                                                          <  1:4267245/150‑1 (MQ=255)
           cTGCGGAGGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                                                                                                          <  2:1143475/150‑1 (MQ=255)
            tGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCg                                                                                                                                         <  2:3944963/150‑1 (MQ=38)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                        <  1:7485472/150‑1 (MQ=38)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                        <  1:5179078/150‑1 (MQ=38)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                        <  2:6674668/150‑1 (MQ=38)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                        <  2:5193842/150‑1 (MQ=38)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                        <  2:1834684/150‑1 (MQ=255)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                        <  1:1078669/150‑1 (MQ=255)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                        <  2:1834040/150‑1 (MQ=255)
             gCGGATGGGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGATTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                        <  2:6948588/150‑1 (MQ=37)
              cGGATGTGCTGCTGACCTACGCACGCTGAAAGTGCTGCCTGTTACAGGAAAACCCCGGCAACTGTGCCGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                       <  1:984683/150‑1 (MQ=25)
              cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                       <  2:2186034/150‑1 (MQ=38)
              cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                       <  2:3422056/150‑1 (MQ=38)
              cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                       <  2:1128998/150‑1 (MQ=38)
              cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                                                                                                       <  1:441274/150‑1 (MQ=255)
              cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                                                                                                       <  1:7977925/150‑1 (MQ=255)
              cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTCCAAGGTTCCTGTAGGAGCggg                                                                                                                                       <  2:5663907/150‑1 (MQ=255)
                     gCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGACAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAccc                                                                                                                                <  1:7160211/150‑1 (MQ=255)
                        gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACAAGGTTGCTGTAGGAGCGGGTTTACCcgcg                                                                                                                             >  2:5157591/1‑150 (MQ=255)
                        gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCGCCGACAAGGGTCCGGTAGGCGCGGGGTTTCCcgcg                                                                                                                             >  1:543628/1‑150 (MQ=255)
                              cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCGGTAGGAGCGGGTTTACCCGCGAAGAgg                                                                                                                       >  2:6677036/1‑150 (MQ=255)
                              cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCGCCGACAAGGTGGCTGGGGGGGCGGGGTTTGCCGCGGCGAgg                                                                                                                       >  1:8064237/1‑150 (MQ=255)
                                    acgcTGAAAGTGCTGCCTTTTTCACTAAACCCCCGGCCACGGGGCGGGGTTTTTTTTTTGGGGTTTTTTTTCGTGTGCCGGCCTGTTCGGGGGTAAACCCGCTCCTACATGGTTCCTGTCGGAGGGGGTTTACCCGCGAAGAGGCCAGca                                                                                                                 <  2:866589/150‑1 (MQ=25)
                                        tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTTAACCCGCTCCGACAAGGTTGCGGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGa                                                                                                            >  2:6310587/1‑150 (MQ=38)
                                        tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAAACGCGACAACAAGGTTACGGTTGGTGCGGGTTTAACCGCGAAGAGGCCCGCACAGa                                                                                                            >  2:4098925/1‑150 (MQ=255)
                                              tgctgcCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                        <  2:949169/117‑1 (MQ=25)
                                                ctgcCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                          >  1:6337009/1‑113 (MQ=25)
                                                 tgcCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGCCTAAAAGGTTCCTGTAGGAGCGGGTGTACCCGCGAAGAGGCCCGCACAGaaaaaaaaaa                                                                                                   >  2:8534553/1‑150 (MQ=37)
                                                           aCGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                    >  1:1703450/1‑92 (MQ=255)
                                                           aCGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                    <  2:1703450/92‑1 (MQ=255)
                                                             gAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGCGGCCAGCACAGAAAAATAAAAATACCTGAtttt                                                                                       >  2:638988/1‑150 (MQ=255)
                                                                aaCCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGTTGCGGTAGGGGCGGGTTTACCCGCGTAGGGGCCCGCACAGAAAAATAAAAATACCAGATTTTCaa                                                                                    >  2:7222595/1‑150 (MQ=37)
                                                                 aCCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCAATTCGCGGGTAAACCCCCAACTACAAGGTTGCTGTGGGAGCGGGGTTAACCGGGAAGGGGCCCGCGCAGAAAAAAAAAAATAGCTTATTTTaaaa                                                                                   >  2:3130475/1‑150 (MQ=255)
                                                                  ccccGGCCAGTGTGCGGGGTTTTTTTTTTGGGGGTTTTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGTTCCTGTCGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaata                                                                                  <  2:154707/150‑3 (MQ=37)
                                                                  ccccGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                        <  2:2098340/97‑1 (MQ=18)
                                                                  ccccGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                        >  1:2098340/1‑97 (MQ=18)
                                                                  ccccGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaa                                                                                  >  2:6088633/1‑150 (MQ=255)
                                                                              ggCCGGGGTTTTTTTTTGGGGGTATTTTGCCTGTGCCGGCCTCTTCGCGGGTAACCCCGCTCCTCCACGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGaa                                                                      <  2:8514670/149‑1 (MQ=255)
                                                                                cgggggTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGGGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGTTTACCCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGaaaa                                                                    <  2:812473/148‑1 (MQ=255)
                                                                                ccGGGGTTTTTTTTTGGGGGTATTTGCCCTGTGCCGGCCGCTTCGCGGGAAAACCCGCTCCTACAAGGTTCCTGTCGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGaaaa                                                                    <  2:7317912/150‑1 (MQ=38)
                                                                                       tttttttttGGGGTTATTTTTCGTGGGCCGGCCTCTTCGCGGGAAACCCCGCTCCTACACGTTGCCTGTCGGAGGGGGTTACCCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCta                                                              <  1:7222804/150‑1 (MQ=37)
                                                                                       tttttttttGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTTATACCGGTGCGTGAAGGTTGCTGTTTGTGGGGGGTTTGCCGCGGAGAGGGCCGCACAGAAAAATAAAAATAACGGATTTTGAAAAAATATATATGAGAATTTTta                                                              >  2:5598263/1‑150 (MQ=255)
                                                                                       tttttttttGGGGGTATTTTCCCTGGGCCGGCCTCTGCGCGGGAAAACCCGCGCCACCACGGTCCCTGTCGGCGCGGGTTTCCCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCta                                                              <  1:7222826/150‑1 (MQ=37)
                                                                                         tttttttGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTCCAACAAGGTTGCAGTAGGAGCGGGGTTACCCCCGAACAGGCCAGCACAGAAAAATAAAAATACCAGATTTTCAAAAAAAATTTATGAAAAGTTCtaat                                                            >  1:7620582/1‑150 (MQ=37)
                                                                                               tGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCAAATAAAACa                                                      <  2:8758848/150‑1 (MQ=255)
                                                                                               tGGGGGTATTTTTCCGGTGCCGGCCACTTCGCGGGTAACCCCGCTCCTACAAGGTTCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                      <  2:6382751/150‑1 (MQ=255)
                                                                                                         tttCGTGGGCCGGCCTCTCCGCGGTTAAACCCGCCCCACCCCGGTGCCTGTCGGCGGGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATc                                            <  1:4130035/150‑1 (MQ=255)
                                                                                                              tgtgCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAg                                                                                                                                           >  2:5531021/1‑51 (MQ=255)
                                                                                                              tgtgCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAg                                                                                                                                           <  1:5531021/51‑1 (MQ=255)
                                                                                                                     gCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCAt                                <  2:4177388/150‑1 (MQ=255)
                                                                                                                             gcgGTTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGTATTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTAGAGGCATCCTGTTTc                        <  2:5521216/150‑1 (MQ=255)
                                                                                                                               gggTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGAAAAATACAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCg                       >  1:673114/1‑150 (MQ=21)
                                                                                                                               gggAAAACCCGCACCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCgg                      <  2:6535491/150‑1 (MQ=255)
                                                                                                                                 gTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCa                    >  1:5116044/1‑150 (MQ=37)
                                                                                                                                 gTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCa                    >  1:1602752/1‑150 (MQ=37)
                                                                                                                                 gTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCa                    >  2:1271137/1‑150 (MQ=255)
                                                                                                                                 gTAAACCCGCTCATACCCGGTCCCTGTCGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCa                    <  1:2808654/150‑1 (MQ=39)
                                                                                                                                   aaaCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATAc                                                                                               >  1:40037/1‑73 (MQ=255)
                                                                                                                                      cccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc               >  2:4058343/1‑150 (MQ=38)
                                                                                                                                      cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc               >  2:3642728/1‑150 (MQ=255)
                                                                                                                                       ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAG‑‑‑‑‑GAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGa         >  2:7264246/1‑150 (MQ=21)
                                                                                                                                       ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct              >  2:3854019/1‑150 (MQ=38)
                                                                                                                                        cGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtctt             >  2:5675711/1‑150 (MQ=255)
                                                                                                                                         gccccTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc            <  1:9053392/147‑1 (MQ=255)
                                                                                                                                         gCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTAAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc            >  1:3265710/1‑150 (MQ=255)
                                                                                                                                           tCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCACGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAAGTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTg          >  2:129133/1‑150 (MQ=255)
                                                                                                                                            caaacacGGTTCCTTTTGGAGCGTGTTTACCCGCGAAGTGGACAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGa         <  1:6677036/147‑1 (MQ=255)
                                                                                                                                             cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTCTCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt        >  1:4359230/1‑150 (MQ=255)
                                                                                                                                             cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt        >  1:4625860/1‑150 (MQ=255)
                                                                                                                                             cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt        >  2:2569365/1‑150 (MQ=255)
                                                                                                                                             cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt        >  2:2316808/1‑150 (MQ=255)
                                                                                                                                              tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt       >  2:5090274/1‑150 (MQ=255)
                                                                                                                                               aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTc      >  1:4013424/1‑150 (MQ=255)
                                                                                                                                               aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTc      >  1:4900102/1‑150 (MQ=255)
                                                                                                                                               aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTc      >  1:6856156/1‑150 (MQ=255)
                                                                                                                                                cccggTTCCTTGAGGTGGGGGTTTACCCGAGAAGAGTCCAGAACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTATTATACAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCGGATTct     <  2:3935046/147‑2 (MQ=255)
                                                                                                                                                cacggTTCCTTTAGGTGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCa     <  1:5157591/147‑1 (MQ=255)
                                                                                                                                                 aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg    >  2:4580513/1‑150 (MQ=255)
                                                                                                                                                 aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg    >  1:1132713/1‑150 (MQ=255)
                                                                                                                                                   ggTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAATATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgc  >  2:763407/1‑150 (MQ=255)
                                                                                                                                                   ggTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgc  >  2:763464/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGC  >  NC_002947/955007‑955302

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: