Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,158 C→T 10.4% intergenic (+109/‑204) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,1580CT10.4% 238.5 / 16.3 96intergenic (+109/‑204)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (42/41);  new base T (3/6);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.97e-02

CCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCCGCCA  >  NC_002947/955010‑955307
                                                                                                                                                        |                                                                                                                                                       
ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                                        <  1:2918796/150‑1 (MQ=255)
 cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTg                                                                                                                                                       <  1:2852808/150‑1 (MQ=255)
 cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTg                                                                                                                                                       <  1:2561845/150‑1 (MQ=255)
 cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCAACAAGGTTccgg                                                                                                                                                       >  2:3937322/1‑148 (MQ=255)
  aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAAGTGCTGCCTGTTTCACGAAAACCCCGGCCTCTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCGGGGCTCTTCGCGGGTAAACCCGCTCCTACCCGGTTCgtgt                                                                                                                                                      <  2:494172/150‑1 (MQ=255)
  aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                                                                                                      <  1:719996/150‑1 (MQ=255)
  aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                                                                                                      <  2:2839473/150‑1 (MQ=255)
    aCCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCGACAAGGTTGCTGTAg                                                                                                                                                    >  2:4261689/1‑150 (MQ=255)
     ccGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAgg                                                                                                                                                   <  2:1485955/150‑1 (MQ=255)
      cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGa                                                                                                                                                  <  2:4312527/150‑1 (MQ=255)
      cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGa                                                                                                                                                  <  2:3010037/150‑1 (MQ=255)
       gCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                               <  1:1886367/150‑1 (MQ=25)
        cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                               <  2:246154/150‑1 (MQ=38)
        cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                               <  2:2903532/150‑1 (MQ=38)
        cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                                                                                                               <  1:475235/150‑1 (MQ=255)
        cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                                                                                                               <  1:3093212/150‑1 (MQ=255)
        cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                                                                                                               <  1:475127/150‑1 (MQ=255)
        cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                                                                                                               <  1:1158589/150‑1 (MQ=255)
         tGCGTATGAGCTGATGACCTACGCACGATGAAACTGCTGCCTGTTTCACTAAAACCCCGCCCACTGTGCCGGGGTATTTTATTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCg                                                                                                                                              <  2:2569575/150‑1 (MQ=38)
         tGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                             <  1:3794698/150‑1 (MQ=25)
          gCGGTTGTGCTGATGCCCTACGCACGCTGAAACTGCGGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGTTTTTTTTTTGGGGGTATTTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                             <  1:4326311/150‑1 (MQ=25)
          gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                             <  2:3590650/150‑1 (MQ=255)
          gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                             <  1:3296377/150‑1 (MQ=255)
           cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                                                                                                            <  1:2237250/150‑1 (MQ=255)
           cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                                                                                                            <  1:4274752/150‑1 (MQ=255)
            ggATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGt                                                                                                                                           <  2:613526/150‑1 (MQ=255)
             gATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGtt                                                                                                                                          <  1:3050452/150‑1 (MQ=255)
                  gCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGACAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAccc                                                                                                                                     <  2:2231241/150‑1 (MQ=255)
                    tgatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCcgc                                                                                                                                   <  2:272561/150‑1 (MQ=37)
                    tgatgaCCGACGCACGCTGAAACTGCTGCCTGTATCACGAAAACCCCGGCCACTGTGCCGGGGTATTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCcgc                                                                                                                                   <  2:3358492/150‑1 (MQ=255)
                     gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCGCCGACAAGGTTCCGGTAGGAGCGGGGTTTACcgcg                                                                                                                                  >  2:3832093/1‑150 (MQ=255)
                          ccTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAATCCCGCTCCTACAAGGTTCCTGGAGGTGCGGGTTTTCCCGCGAagag                                                                                                                             >  2:3363419/1‑150 (MQ=255)
                          ccTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCGCCAACAAGGGTGCGGTAGGAGCGGGGTTTCCCCCGCagag                                                                                                                             >  1:2253098/1‑150 (MQ=255)
                          ccTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTGTGGGGGTATTTTTCCAGTGCCGGCCTCTTCGCGGGGAAAACCCCAACAACGAAGTTACGGGAGGGGCGGGTGTTCCcgcgtaaag                                                                                                                             >  1:2252846/1‑145 (MQ=255)
                           cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTGAACCCCCGCCGACAAGGTTCCGGTAGGTGCGGGGTTTCCCGCGGAGAgg                                                                                                                            >  1:714139/1‑150 (MQ=255)
                           cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCGGTAGGGGCGGGTTTACCCGCGAAGAgg                                                                                                                            >  2:713829/1‑150 (MQ=255)
                                      gAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCGCCTACAAGGTTCCTGGAGGAGCGGGTTTACCCGCGCAGAGGCCAGCGCCGaa                                                                                                                >  2:1746550/1‑150 (MQ=38)
                                         aCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaa                                                                                                             <  2:2560407/150‑1 (MQ=255)
                                         aCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTGTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaa                                                                                                             <  2:2635085/150‑1 (MQ=255)
                                         aCTGCTGCCTGTCGCACGAAACCCCCGGCCAGTGTGCGGGGTTTTTTTTTGGGGGGTATTTTTCGGGGGCCGGCCTCTGGGCGGGTAAACCCGCTCCTCCCCGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaa                                                                                                             <  2:3810503/150‑1 (MQ=25)
                                              tgcCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGTAACCCCCTCCGACAAGGTTACTGTAGGAGCGGGTTTTCCCGCGGAGAGGCCAGCCCCGacaaataaaa                                                                                                        >  2:245371/1‑150 (MQ=255)
                                                ccTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTAGCGGGTTAACACCCGCCGACAAGGTTGCAGGTGGAGGGGGGTTTACCCCGAAGGGGGCCGCACACgcaaataaaaat                                                                                                      >  1:1793927/1‑150 (MQ=255)
                                                ccTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTATTCGCGGGTAAACCCGCTCCTACCAGGTTCCGGTAGGAGCGGGTTTACCCGCGTAGTGGCCAGCAAAGaaaaataaaaat                                                                                                      >  2:364766/1‑150 (MQ=38)
                                                  tGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTAGCGGGGTATACCCCGCCCACCACGTTGGGGGTGGGGGGGGGGTTCCCGCGGAGAGGCCCGCCCAGcaaaaaaaaaaaaa                                                                                                    >  1:3129213/1‑147 (MQ=255)
                                                  tGTTTCACGAAAACCCCGGCCACGGTGCCGGGGTTTTTTTTGGGGGGTATTTTTCCGGTGCCGGCCTCTCCGCGGGTAAACCCGCCCCACCCCGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATAc                                                                                                    <  2:856846/150‑1 (MQ=37)
                                                    tttCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCACCTACAAGGTTGCTGTATGAGCGGGTTTTCCCGCGAAGAGGCCAGCAAAGAAAAATAAAAATACCt                                                                                                  >  2:4223411/1‑150 (MQ=38)
                                                    ccaaaaaaaaaCCCCGGCACAGTGGCCGGGGTTTTTTTTTTGGTGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGCCTTCAAGGTTCCGGTAGGGGCGGGGTTACCCGCGAAGAGGCCAGCACAGGAAAATAAAAATAcc                                                                                                   >  2:2233506/8‑150 (MQ=21)
                                                      tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCACCTACAAAGTTCCGGTAGGAGCGGGTGTACCCCCCCACAGGCCCGCAACGAAAAATAAAAATAccagg                                                                                                >  2:789220/1‑147 (MQ=37)
                                                        aCGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCTCCGACCAGGTTGCGGGTGGGGCGGGGTTAACCGCGCAGCGGCCAGCACCGAAAAATAAAAAAACCTGCtt                                                                                              >  2:4159108/1‑150 (MQ=255)
                                                             aaCCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCGTCGCGGGGAAACCCGCACCCACAAGGGTCCGGGAGGCGCGGGGGTACCCGCGCAGAGGCCCCCCCCGCAAAAAAAAAAAACCGGATTTTCaa                                                                                         >  2:1956709/1‑150 (MQ=255)
                                                               ccccGGCCCCTGTGCGGGTGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACGGTTCCTTTTGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaa                                                                                       <  1:734189/150‑1 (MQ=37)
                                                                       aCTGTGCCGGGGATTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCTCCTACAAGGGTGCTGTAGGGGCGGGTTTTCCCGCGACGAGGCCAGCACAGAAAAAAAAAAATACCTGATTTTGAAAAAATAGAta                                                                               >  2:2791069/1‑149 (MQ=37)
                                                                        cTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCt                                                                              <  2:3381908/150‑1 (MQ=255)
                                                                                   gttttttttttGGGGTTATTTTTCCTGTGCCGGCCTCTTCGCGGTTAAACCCGCTCCTACACGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc                                                                     <  2:429023/148‑1 (MQ=255)
                                                                                    ttttttttttGGGGTTATTTTTCCGGTGCCGGCCTCTTCGCGGTTAAACCCGCTCCTACAAGGTTCCTGTCGGAGCGGGTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCt                                                                    <  2:1877266/149‑1 (MQ=38)
                                                                                        ttttttGGGGGTTTTTTTCCCGGTGCCGCCCTCTTCGCGGGAACCCCCGCTCCTACACGGTACGTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaat                                                                 <  1:2702163/150‑1 (MQ=37)
                                                                                             tGGGGGTATTTTCCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                           <  2:1028763/150‑1 (MQ=255)
                                                                                                 ggTATTTTTCCTGTGCCGGCCTCTTCGCGGGTTAACCCGCTCCTACAAGGTTGCTGGAGGAGCGGGTTTAACCGCGGAGGGGCCAGCACCGGAAAATAAAAATACCTGTTTTTAAAAAAATAGGTATTAAAAGTTGTAATAAAACAGGaa                                                       >  2:3510506/1‑150 (MQ=25)
                                                                                                 ggTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAGCCGCTCCTACAATGTTGCTGTAGGTGCGGGTTTAGCCGCGTAGAGGCCAGCACAGAAAAATTAAAATACCTGATTTTGTAAAAATAGATCTGAAAAGTTGTAATAAAAGAGGaa                                                       >  2:1441133/1‑150 (MQ=25)
                                                                                                    aTTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAAtt                                                    <  2:3075426/150‑1 (MQ=255)
                                                                                                    aTTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAAtt                                                    <  2:3075433/150‑1 (MQ=255)
                                                                                                           cgtgtgCCGCCCTCTGGGGGGGTAACCCCGCTCCTCCACGGTTCCTGTTGGCGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCtt                                              <  1:3216641/148‑1 (MQ=255)
                                                                                                                         ttCGCGGGTAAACCCGCGCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGTTCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGt                                >  2:1183987/1‑150 (MQ=255)
                                                                                                                           cgcgGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGttt                              >  2:3702639/1‑150 (MQ=255)
                                                                                                                            gcgGGTAAACCCGCTCCTACACGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTc                             <  2:77860/150‑1 (MQ=255)
                                                                                                                            gcgGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAGTACAGGAAATTATCCTTTACAGGCATCCTGTTTc                             >  2:3148245/1‑150 (MQ=255)
                                                                                                                            gcgGGAAAACCCGCTCCTACAAGGTTCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTc                             <  1:713829/150‑1 (MQ=255)
                                                                                                                             ggggttcaaCCCGATCCTACACGGTTCCTGTTGGTGGGTGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTc                             <  1:2790275/143‑1 (MQ=255)
                                                                                                                                 gtaatgccGCTCCTACAAGGTTCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAATTAGATCTGAATAGTTCTAATAATACAGGTAATTATCGTTTACAGGCATGCTGTTTCGGct                         >  2:2341787/7‑149 (MQ=38)
                                                                                                                                   aaaCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACt                       >  2:3899452/1‑150 (MQ=255)
                                                                                                                                    aaCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTa                      >  1:1241004/1‑150 (MQ=255)
                                                                                                                                      cccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc                    >  1:1961007/1‑150 (MQ=38)
                                                                                                                                       ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct                   >  1:3496614/1‑150 (MQ=38)
                                                                                                                                       ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct                   >  2:1238897/1‑150 (MQ=38)
                                                                                                                                       ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc                    >  1:1363848/1‑150 (MQ=255)
                                                                                                                                       ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct                   >  1:1156226/1‑150 (MQ=255)
                                                                                                                                       ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct                   >  1:4346146/1‑150 (MQ=255)
                                                                                                                                       ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct                   >  2:2525492/1‑150 (MQ=255)
                                                                                                                                       ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct                   >  2:378955/1‑150 (MQ=255)
                                                                                                                                        cGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtctt                  >  2:2354765/1‑150 (MQ=255)
                                                                                                                                         gcacacacacGGTCCCGGTCGGCGGGGGTTTCCCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc                 <  1:3832093/144‑1 (MQ=255)
                                                                                                                                         gCTCCTACACGGTTCCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc                 <  2:1331475/150‑1 (MQ=255)
                                                                                                                                         gCTCCTACACGGTACGTGTAGTCGGGGGTTTCCCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc                 <  1:3937322/150‑1 (MQ=255)
                                                                                                                                         gCTCCTACAAGGTTCCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc                 <  2:2814860/150‑1 (MQ=255)
                                                                                                                                         gCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc                 >  1:1556563/1‑150 (MQ=255)
                                                                                                                                             cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt             >  2:643587/1‑150 (MQ=255)
                                                                                                                                              tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt            >  1:1209208/1‑150 (MQ=255)
                                                                                                                                              tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt            >  1:266595/1‑150 (MQ=255)
                                                                                                                                              tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAta            >  2:2952856/1‑149 (MQ=255)
                                                                                                                                               aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTGCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTATGATTc           >  1:3136085/1‑150 (MQ=255)
                                                                                                                                                cAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCa          >  1:1119233/1‑150 (MQ=255)
                                                                                                                                                cAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCa          >  2:2285089/1‑150 (MQ=255)
                                                                                                                                                 aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg         >  2:3241215/1‑150 (MQ=255)
                                                                                                                                                   ggTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgc       >  2:3757116/1‑150 (MQ=255)
                                                                                                                                                   ggTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgc       >  2:2462501/1‑150 (MQ=255)
                                                                                                                                                   ggTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATCATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgc       >  2:2858249/1‑150 (MQ=255)
                                                                                                                                                      tCCTTTGGGAGGGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGCAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgccgc    <  1:3998760/150‑1 (MQ=255)
                                                                                                                                                      tCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgccgc    <  2:714139/150‑1 (MQ=255)
                                                                                                                                                        cTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCCGCca  >  1:893569/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                       
CCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCCGCCA  >  NC_002947/955010‑955307

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: