Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,154 G→A 18.5% intergenic (+105/‑208) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,1540GA18.5% 358.8 / 71.1 165intergenic (+105/‑208)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (56/76);  new base A (10/20);  total (66/99)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-02

Only 100 of 170 total aligned reads displayed.

CAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGC  >  NC_002947/955007‑955302
                                                                                                                                                    |                                                                                                                                                    
cAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTTTTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGtt                                                                                                                                                    <  1:5073726/150‑1 (MQ=255)
 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)
  gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTcc                                                                                                                                                  <  1:4695193/150‑1 (MQ=255)
  gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTcc                                                                                                                                                  <  1:4553698/150‑1 (MQ=255)
  gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTcc                                                                                                                                                  <  1:4215473/150‑1 (MQ=255)
   ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                                 <  2:1221622/150‑1 (MQ=255)
    cAGACCGCTGCGTATTTGCTGATGACCTACGCACGCTTAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTg                                                                                                                                                <  1:4375269/150‑1 (MQ=255)
    cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTg                                                                                                                                                <  2:8401238/150‑1 (MQ=255)
     aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                                                                                               <  1:6944607/150‑1 (MQ=255)
     aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                                                                                               <  1:6944138/150‑1 (MQ=255)
       aCCGCTGCTGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCAGGCCACTGTGCCGGGGTTTTTTTTTTGGGGTATTTTTCCTGTGCCGGCCTATTCTCGTGTAAACCCTCTCCTACAAGGTTCCTGTAg                                                                                                                                             <  1:2347541/150‑1 (MQ=255)
       aCCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAg                                                                                                                                             <  1:8727706/150‑1 (MQ=255)
         cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAg                                                                                                                                          <  2:5345053/150‑1 (MQ=25)
          gCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAg                                                                                                                                          <  2:2469325/150‑1 (MQ=38)
          gCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAg                                                                                                                                          <  1:5335330/150‑1 (MQ=255)
           cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                         <  2:2219955/150‑1 (MQ=38)
           cTGCGGAGGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                                                                                                         <  2:1143475/150‑1 (MQ=255)
            tGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCg                                                                                                                                        <  2:3944963/150‑1 (MQ=38)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                       <  1:8560577/150‑1 (MQ=38)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                       <  2:5193842/150‑1 (MQ=38)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                       <  2:6674668/150‑1 (MQ=38)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                       <  1:1078669/150‑1 (MQ=255)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                       <  1:1281086/150‑1 (MQ=255)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                       <  2:1834040/150‑1 (MQ=255)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                       <  1:7550296/150‑1 (MQ=255)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                       <  2:1834684/150‑1 (MQ=255)
             gCGGATGGGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGATTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                       <  2:6948588/150‑1 (MQ=37)
              cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                      <  2:2186034/150‑1 (MQ=38)
              cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                      <  2:3422056/150‑1 (MQ=38)
              cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                                                                                                      <  1:7977925/150‑1 (MQ=255)
              cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTCCAAGGTTCCTGTAGGAGCggg                                                                                                                                      <  2:5663907/150‑1 (MQ=255)
                       tgatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCcgc                                                                                                                             <  1:3483475/150‑1 (MQ=255)
                        gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACAAGGTTGCTGTAGGAGCGGGTTTACCcgcg                                                                                                                            >  2:5157591/1‑150 (MQ=255)
                              cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCGGTAGGAGCGGGTTTACCCGCGAAGAgg                                                                                                                      >  2:6677036/1‑150 (MQ=255)
                              cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCGACAAGGTTCCTGTAGGAGCGGGGTTAACCGCGCAGAgg                                                                                                                      >  1:8064219/1‑150 (MQ=39)
                                 cgcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCGGGGTTTTTTTTTTGGGGGTATTTTACCGGGGCCGGCCTCTTCGCGGGTAACCCCGCTCCTACACGTTTCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCa                                                                                                                   <  1:6790377/150‑1 (MQ=37)
                                    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)
                                              tgctgcCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                       >  1:949169/1‑117 (MQ=25)
                                                 tgcCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGCCTAAAAGGTTCCTGTAGGAGCGGGTGTACCCGCGAAGAGGCCCGCACAGaaaaaaaaaa                                                                                                   >  2:8534553/1‑150 (MQ=37)
                                                   ccTGTTTCACGAAATCCCCGGCCACTGTGCCGGGGTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGAAACCCGCTCCGAAAAGGTTGCGGTAGGGGCGGGGTTACCCGCGGAGAGGCCAGCACAGAaaaataaaaata                                                                                                >  1:2190216/1‑150 (MQ=37)
                                                           aCGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                   >  1:1703450/1‑92 (MQ=255)
                                                           aCGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                   <  2:1703450/92‑1 (MQ=255)
                                                            aaaaaaCCCCGGCACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCAACAAGGTTCCTGTAGGGGCGGGGTTAACCGCGAAGAGGCCCGCACAGAAAAATAAAAATACCTGAttt                                                                                        >  2:1143497/3‑150 (MQ=37)
                                                            aaaaaaCCCCGGCACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGAAAACCCCGCCGACAAGGGTGCGGGGGGTGCGGGTTTTACCGCGAAGAGGCCAGCACCCAAAAAAAAAAATAACGGAttt                                                                                        >  1:6711199/3‑150 (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)
                                                                  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)
                                                                                    ggggTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGAAAACCCCGCCGACAAGGTTGCGGTTGGGGCGGGGTTTCCCGCGGAGTGGCCCGCCCACAAAAATAAAAATAACGGTTGTTGGAAAAAAAGATGAGAAAAggt                                                                 >  2:1053851/1‑148 (MQ=255)
                                                                                       tttttttttGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTTATACCGGTGCGTGAAGGTTGCTGTTTGTGGGGGGTTTGCCGCGGAGAGGGCCGCACAGAAAAATAAAAATAACGGATTTTGAAAAAATATATATGAGAATTTTta                                                              >  2:5598263/1‑150 (MQ=255)
                                                                                       tttttttttGGGGGTATTTTCCCTGGGCCGGCCTCTGCGCGGGAAAACCCGCGCCACCACGGTCCCTGTCGGCGCGGGTTTCCCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCta                                                              <  1:7222826/150‑1 (MQ=37)
                                                                                               tGGGGTTTTTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCATACACGGTACCTGTAGTTGGGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                      <  1:7128212/150‑1 (MQ=37)
                                                                                               tGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCAAATAAAACa                                                      <  2:8758848/150‑1 (MQ=255)
                                                                                               tGGGGGTATTTTTCCGGTGCCGGCCACTTCGCGGGTAACCCCGCTCCTACAAGGTTCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                      <  2:6382751/150‑1 (MQ=255)
                                                                                                  ggttATTTTTCCTGTGCCGGCCTCTTCGCGGTAAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGa                                                   <  1:7464070/147‑1 (MQ=255)
                                                                                                      aTTTTCCCTGTGCCGGCATCTTCGCGTTAAACCCCGCTCCAAAATGGTTCCTGTAGGAGGGGTTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAAtt                                               <  1:4836678/150‑1 (MQ=255)
                                                                                                              tgtgCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAg                                                                                                                                          >  2:5531021/1‑51 (MQ=255)
                                                                                                              tgtgCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAg                                                                                                                                          <  1:5531021/51‑1 (MQ=255)
                                                                                                                tgCCGGCCTCGTCGCGGGTAAACCCGCTCCTCCACGTTTCCTCTCGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACa                                     <  1:5174033/150‑1 (MQ=38)
                                                                                                                     gCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCAt                                <  2:4177388/150‑1 (MQ=255)
                                                                                                                            ggcgGTTAAACCCGCTCCACCATGGTTCCTGTTGGCGGGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGttt                         <  1:1971523/149‑1 (MQ=255)
                                                                                                                             gcgGTTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGTATTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTAGAGGCATCCTGTTTc                        <  2:5521216/150‑1 (MQ=255)
                                                                                                                               gggAAAACCCGCACCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCgg                      <  2:6535491/150‑1 (MQ=255)
                                                                                                                                ggTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGc                     >  1:2828881/1‑150 (MQ=255)
                                                                                                                                 gTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCa                    >  2:1271137/1‑150 (MQ=255)
                                                                                                                                 gTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCa                    >  1:717846/1‑150 (MQ=255)
                                                                                                                                      cccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc               >  2:4058343/1‑150 (MQ=38)
                                                                                                                                      cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc               >  1:2543195/1‑150 (MQ=255)
                                                                                                                                      cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc               >  2:3642728/1‑150 (MQ=255)
                                                                                                                                       ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAG‑‑‑‑‑GAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGa         >  2:7264246/1‑150 (MQ=21)
                                                                                                                                       ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct              >  2:3854019/1‑150 (MQ=38)
                                                                                                                                       ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct              >  1:6189083/1‑150 (MQ=255)
                                                                                                                                        cGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtctt             >  2:5675711/1‑150 (MQ=255)
                                                                                                                                        cGCTCCCACACGGTTCCTGTCGGTGGGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtctt             <  1:5762143/150‑1 (MQ=255)
                                                                                                                                         gccccTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc            <  1:9053392/147‑1 (MQ=255)
                                                                                                                                           tCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCACGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAAGTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTg          >  2:129133/1‑150 (MQ=255)
                                                                                                                                            caaacacGGTTCCTTTTGGAGCGTGTTTACCCGCGAAGTGGACAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGa         <  1:6677036/147‑1 (MQ=255)
                                                                                                                                             cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTCTCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt        >  1:4359230/1‑150 (MQ=255)
                                                                                                                                             cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt        >  2:2569365/1‑150 (MQ=255)
                                                                                                                                             cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt        >  2:2316808/1‑150 (MQ=255)
                                                                                                                                              tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt       >  1:4038707/1‑150 (MQ=255)
                                                                                                                                              tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt       >  2:5090274/1‑150 (MQ=255)
                                                                                                                                               aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTc      >  1:4013424/1‑150 (MQ=255)
                                                                                                                                                cccggTTCCTTGAGGTGGGGGTTTACCCGAGAAGAGTCCAGAACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTATTATACAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCGGATTct     <  2:3935046/147‑2 (MQ=255)
                                                                                                                                                cacggTTCCTTTAGGTGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCa     <  1:5157591/147‑1 (MQ=255)
                                                                                                                                                 aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg    >  1:1132713/1‑150 (MQ=255)
                                                                                                                                                 aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg    >  2:4580513/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: