Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,152 A→G 19.7% intergenic (+103/‑210) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,1520AG19.7% 304.4 / 69.1 163intergenic (+103/‑210)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (55/59);  new base G (10/20);  total (65/98)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-01

Only 100 of 169 total aligned reads displayed.

TGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAG  >  NC_002947/955003‑955300
                                                                                                                                                      |                                                                                                                                                    
tgGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                      <  1:4757085/150‑1 (MQ=255)
tgGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                      <  2:3094305/150‑1 (MQ=255)
tgGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTAATTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                      <  1:7381077/150‑1 (MQ=255)
    cAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTTTTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGtt                                                                                                                                                  <  1:5073726/150‑1 (MQ=255)
     cgCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                 <  1:4646991/149‑1 (MQ=255)
     aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                 <  2:7777013/150‑1 (MQ=255)
     aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                 <  2:7597172/150‑1 (MQ=255)
     aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                 <  2:7597133/150‑1 (MQ=255)
     aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                 <  2:2831330/150‑1 (MQ=255)
     aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCACCAACAAgggtg                                                                                                                                                 >  2:1196189/1‑147 (MQ=255)
      gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTcc                                                                                                                                                <  1:4695193/150‑1 (MQ=255)
      gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTcc                                                                                                                                                <  1:4215473/150‑1 (MQ=255)
      gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTcc                                                                                                                                                <  1:4553698/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)
        cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                                                                                             <  2:2627821/150‑1 (MQ=255)
         aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                                                                                             <  1:6944581/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)
               cTGCGGATGAGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                                                                                                       <  1:4267245/150‑1 (MQ=255)
               cTGCGGAGGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                                                                                                       <  2:1143475/150‑1 (MQ=255)
                tGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCg                                                                                                                                      <  2:3944963/150‑1 (MQ=38)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                     <  1:8560577/150‑1 (MQ=38)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                     <  1:7485472/150‑1 (MQ=38)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                     <  2:5193842/150‑1 (MQ=38)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                     <  2:1834684/150‑1 (MQ=255)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                     <  2:1834040/150‑1 (MQ=255)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                     <  1:7550296/150‑1 (MQ=255)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                     <  1:1078669/150‑1 (MQ=255)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                     <  1:1281086/150‑1 (MQ=255)
                 gCGGATGGGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGATTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                     <  2:6948588/150‑1 (MQ=37)
                  cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                    <  2:3422056/150‑1 (MQ=38)
                  cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                    <  2:2186034/150‑1 (MQ=38)
                  cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                                                                                                    <  1:7977925/150‑1 (MQ=255)
                  cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                                                                                                    <  2:3926332/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)
                            gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCGCCGACAAGGGTCCGGTAGGCGCGGGGTTTCCcgcg                                                                                                                          >  1:543628/1‑150 (MQ=255)
                                  cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCGGTAGGAGCGGGTTTACCCGCGAAGAgg                                                                                                                    >  2:6677036/1‑150 (MQ=255)
                                  cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCGCCGACAAGGTGGCTGGGGGGGCGGGGTTTGCCGCGGCGAgg                                                                                                                    >  1:8064237/1‑150 (MQ=255)
                                     cgcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCGGGGTTTTTTTTTTGGGGGTATTTTACCGGGGCCGGCCTCTTCGCGGGTAACCCCGCTCCTACACGTTTCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCa                                                                                                                 <  1:6790377/150‑1 (MQ=37)
                                        acgcTGAAAGTGCTGCCTTTTTCACTAAACCCCCGGCCACGGGGCGGGGTTTTTTTTTTGGGGTTTTTTTTCGTGTGCCGGCCTGTTCGGGGGTAAACCCGCTCCTACATGGTTCCTGTCGGAGGGGGTTTACCCGCGAAGAGGCCAGca                                                                                                              <  2:866589/150‑1 (MQ=25)
                                        acgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTACTGTCGGAGCGGGGTTACCCGCGAAGAGGCCAGca                                                                                                              >  2:1332288/1‑150 (MQ=255)
                                            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)
                                                    ctgcCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                       <  2:6337009/113‑1 (MQ=25)
                                                    ctgcCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                       >  1:6337009/1‑113 (MQ=25)
                                                     tgcCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGCCTAAAAGGTTCCTGTAGGAGCGGGTGTACCCGCGAAGAGGCCCGCACAGaaaaaaaaaa                                                                                                 >  2:8534553/1‑150 (MQ=37)
                                                       ccTGTTTCACGAAATCCCCGGCCACTGTGCCGGGGTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGAAACCCGCTCCGAAAAGGTTGCGGTAGGGGCGGGGTTACCCGCGGAGAGGCCAGCACAGAaaaataaaaata                                                                                              >  1:2190216/1‑150 (MQ=37)
                                                               aCGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                 <  2:1703450/92‑1 (MQ=255)
                                                               aCGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                 >  1:1703450/1‑92 (MQ=255)
                                                                aaaaaaCCCCGGCACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGAAAACCCCGCCGACAAGGGTGCGGGGGGTGCGGGTTTTACCGCGAAGAGGCCAGCACCCAAAAAAAAAAATAACGGAttt                                                                                      >  1:6711199/3‑150 (MQ=255)
                                                                 gAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGCGGCCAGCACAGAAAAATAAAAATACCTGAtttt                                                                                     >  2:638988/1‑150 (MQ=255)
                                                                    aaCCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGTTGCGGTAGGGGCGGGTTTACCCGCGTAGGGGCCCGCACAGAAAAATAAAAATACCAGATTTTCaa                                                                                  >  2:7222595/1‑150 (MQ=37)
                                                                      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)
                                                                                           tttttttttGGGGTTATTTTTCGTGGGCCGGCCTCTTCGCGGGAAACCCCGCTCCTACACGTTGCCTGTCGGAGGGGGTTACCCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCta                                                            <  1:7222804/150‑1 (MQ=37)
                                                                                           tttttttttGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTTATACCGGTGCGTGAAGGTTGCTGTTTGTGGGGGGTTTGCCGCGGAGAGGGCCGCACAGAAAAATAAAAATAACGGATTTTGAAAAAATATATATGAGAATTTTta                                                            >  2:5598263/1‑150 (MQ=255)
                                                                                            ttttttttGGGGTTATTTTTCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaa                                                           <  2:3168885/150‑1 (MQ=255)
                                                                                                   tGGGGTTTTTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCATACACGGTACCTGTAGTTGGGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                    <  1:7128212/150‑1 (MQ=37)
                                                                                                   tGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCAAATAAAACa                                                    <  2:8758848/150‑1 (MQ=255)
                                                                                                                  tgtgCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAg                                                                                                                                        >  2:5531021/1‑51 (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)
                                                                                                                                  cgGTTAAACCCGCTCCTACAAGGTTCCTGGAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCg                     <  2:3851653/150‑1 (MQ=255)
                                                                                                                                   gggAAAACCCGCACCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCgg                    <  2:6535491/150‑1 (MQ=255)
                                                                                                                                    ggTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGc                   >  1:2828881/1‑150 (MQ=255)
                                                                                                                                     gTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCa                  >  1:717846/1‑150 (MQ=255)
                                                                                                                                     gTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCa                  >  2:1271137/1‑150 (MQ=255)
                                                                                                                                       aaaCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATAc                                                                                             >  1:40037/1‑73 (MQ=255)
                                                                                                                                          cccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc             >  2:4058343/1‑150 (MQ=38)
                                                                                                                                          cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc             >  1:2543195/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAG‑‑‑‑‑GAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGa       >  2:7264246/1‑150 (MQ=21)
                                                                                                                                           ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct            >  2:6584388/1‑150 (MQ=255)
                                                                                                                                           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)
                                                                                                                                                 cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt      >  2:2316808/1‑150 (MQ=255)
                                                                                                                                                  tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt     >  1:4038707/1‑150 (MQ=255)
                                                                                                                                                  tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt     >  2:5090274/1‑150 (MQ=255)
                                                                                                                                                   aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTc    >  1:4900102/1‑150 (MQ=255)
                                                                                                                                                     aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg  >  2:4580513/1‑150 (MQ=255)
                                                                                                                                                     aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg  >  1:1132713/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAG  >  NC_002947/955003‑955300

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: