Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,156 T→C 19.4% intergenic (+107/‑206) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,1560TC19.4% 360.9 / 70.4 168intergenic (+107/‑206)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (57/75);  new base C (10/22);  total (68/100)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

CAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCC  >  NC_002947/955007‑955303
                                                                                                                                                       |                                                                                                                                                   
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)
 aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTATTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                    <  2:2831444/150‑1 (MQ=255)
  gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTcc                                                                                                                                                   <  1:4695193/150‑1 (MQ=255)
  gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTcc                                                                                                                                                   <  1:4215473/150‑1 (MQ=255)
  gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTcc                                                                                                                                                   <  1:4553698/150‑1 (MQ=255)
    cAGACCGCTGCGTATTTGCTGATGACCTACGCACGCTTAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTg                                                                                                                                                 <  1:4375269/150‑1 (MQ=255)
    cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                                  <  1:1178731/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)
       aCCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAg                                                                                                                                              <  1:8727706/150‑1 (MQ=255)
       aCCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAg                                                                                                                                              <  1:2347217/150‑1 (MQ=255)
         cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAg                                                                                                                                           <  2:5345053/150‑1 (MQ=25)
          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                                                                                                                                        <  2:1688131/150‑1 (MQ=38)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                        <  2:5193842/150‑1 (MQ=38)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                        <  1:7550296/150‑1 (MQ=255)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                        <  2:1834684/150‑1 (MQ=255)
             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                        <  1:1078669/150‑1 (MQ=255)
             gCGGATGGGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGATTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                        <  2:6948588/150‑1 (MQ=37)
              cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                                                                                                       <  2:3926332/150‑1 (MQ=255)
              cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                                                                                                       <  1:7977925/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)
                                              tgctgcCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                        >  1:949169/1‑117 (MQ=25)
                                              tgctgcCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                        <  2:949169/117‑1 (MQ=25)
                                                ctgcCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                          <  2:6337009/113‑1 (MQ=25)
                                                 tgcCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGCCTAAAAGGTTCCTGTAGGAGCGGGTGTACCCGCGAAGAGGCCCGCACAGaaaaaaaaaa                                                                                                    >  2:8534553/1‑150 (MQ=37)
                                                   ccTGTTTCACGAAATCCCCGGCCACTGTGCCGGGGTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGAAACCCGCTCCGAAAAGGTTGCGGTAGGGGCGGGGTTACCCGCGGAGAGGCCAGCACAGAaaaataaaaata                                                                                                 >  1:2190216/1‑150 (MQ=37)
                                                         tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGGTCCTGGAGGAGCGGGTTTAACCGCGAAGAGGCCCGCACAGAAAAATAAAAATACCTGa                                                                                            >  2:2172061/1‑150 (MQ=255)
                                                            aaaaaaCCCCGGCACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCAACAAGGTTCCTGTAGGGGCGGGGTTAACCGCGAAGAGGCCCGCACAGAAAAATAAAAATACCTGAttt                                                                                         >  2:1143497/3‑150 (MQ=37)
                                                            aaaaaaCCCCGGCACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGAAAACCCCGCCGACAAGGGTGCGGGGGGTGCGGGTTTTACCGCGAAGAGGCCAGCACCCAAAAAAAAAAATAACGGAttt                                                                                         >  1:6711199/3‑150 (MQ=255)
                                                             gAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGCGGCCAGCACAGAAAAATAAAAATACCTGAtttt                                                                                        >  2:638988/1‑150 (MQ=255)
                                                              aaaaCCCCGGCACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTGCTGGAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGAttttg                                                                                       >  2:2078567/1‑149 (MQ=37)
                                                                aaCCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGTTGCGGTAGGGGCGGGTTTACCCGCGTAGGGGCCCGCACAGAAAAATAAAAATACCAGATTTTCaa                                                                                     >  2:7222595/1‑150 (MQ=37)
                                                                  ccccGGCCAGTGTGCGGGGTTTTTTTTTTGGGGGTTTTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGTTCCTGTCGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaata                                                                                   <  2:154707/150‑3 (MQ=37)
                                                                  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)
                                                                                         ttttttttGGGGTTATTTTTCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaa                                                              <  2:3168885/150‑1 (MQ=255)
                                                                                                tgcgggTATTTTTCCTGAGCCGCCCACTTCGCGGGTAAACCCGCACAAACAACGTTCCTCAAGCACCGGCTTTACCCGCGAAGAGGCCACCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                       <  2:4653713/147‑1 (MQ=255)
                                                                                                tGGGGTTTTTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCATACACGGTACCTGTAGTTGGGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                       <  1:7128212/150‑1 (MQ=37)
                                                                                                tGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCAAATAAAACa                                                       <  2:8758848/150‑1 (MQ=255)
                                                                                                   ggttATTTTTCCTGTGCCGGCCTCTTCGCGGTAAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGa                                                    <  1:7464070/147‑1 (MQ=255)
                                                                                                   gggTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCTCCAAACAGGTTACGGTAGGAGCGGGTTTTACCCCGAAGGGGCCAGCCCAGAAAAAAAAAAATAACTGATTTTAAAAAAATAGATATGTAAATTTATAATAAAACAGGa                                                    >  2:1224270/1‑150 (MQ=255)
                                                                                                       aTTTTCCCTGTGCCGGCATCTTCGCGTTAAACCCCGCTCCAAAATGGTTCCTGTAGGAGGGGTTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAAtt                                                <  1:4836678/150‑1 (MQ=255)
                                                                                                               tgtgCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAg                                                                                                                                           <  1:5531021/51‑1 (MQ=255)
                                                                                                               tgtgCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAg                                                                                                                                           >  2:5531021/1‑51 (MQ=255)
                                                                                                                 tgCCGGCCTGTTCGGGGGTAAACCCGCTCCTACATGGTTCCTGTTGGAGGGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACa                                      <  2:4076384/150‑1 (MQ=38)
                                                                                                                 tgCCGGCCTCGTCGCGGGTAAACCCGCTCCTCCACGTTTCCTCTCGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACa                                      <  1:5174033/150‑1 (MQ=38)
                                                                                                                      gCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCAt                                 <  2:4177388/150‑1 (MQ=255)
                                                                                                                             ggcgGTTAAACCCGCTCCACCATGGTTCCTGTTGGCGGGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGttt                          <  1:1971523/149‑1 (MQ=255)
                                                                                                                                gggAAAACCCGCACCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCgg                       <  2:6535491/150‑1 (MQ=255)
                                                                                                                                 ggTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGc                      >  1:2828881/1‑150 (MQ=255)
                                                                                                                                  gTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCa                     >  1:5116044/1‑150 (MQ=37)
                                                                                                                                  gTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCa                     >  1:1602752/1‑150 (MQ=37)
                                                                                                                                  gTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCa                     >  1:717846/1‑150 (MQ=255)
                                                                                                                                    aaaCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATAc                                                                                                >  1:40037/1‑73 (MQ=255)
                                                                                                                                       cccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc                >  2:4058343/1‑150 (MQ=38)
                                                                                                                                       cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc                >  2:1744399/1‑150 (MQ=255)
                                                                                                                                       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)
                                                                                                                                        ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct               >  2:6584388/1‑150 (MQ=255)
                                                                                                                                         cGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtctt              >  2:2333714/1‑150 (MQ=38)
                                                                                                                                         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)
                                                                                                                                             ccTCCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGa          <  2:543628/150‑1 (MQ=255)
                                                                                                                                             caaacacGGTTCCTTTTGGAGCGTGTTTACCCGCGAAGTGGACAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGa          <  1:6677036/147‑1 (MQ=255)
                                                                                                                                              cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt         >  2:2569365/1‑150 (MQ=255)
                                                                                                                                              cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt         >  2:2316808/1‑150 (MQ=255)
                                                                                                                                               tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt        >  1:4038707/1‑150 (MQ=255)
                                                                                                                                               tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGAATTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTAAACAGGCATCCTGTTACGGCACTATCTTCTGAta        >  2:5639099/1‑149 (MQ=255)
                                                                                                                                                  aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg     >  2:4580513/1‑150 (MQ=255)
                                                                                                                                                  aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg     >  1:1132713/1‑150 (MQ=255)
                                                                                                                                                    ggTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgc   >  2:763464/1‑150 (MQ=255)
                                                                                                                                                     ggtCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAACACAGGAAAGTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgcc  >  1:1332583/3‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
CAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCC  >  NC_002947/955007‑955303

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