Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,152 A→G 14.0% 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,1520AG14.0% 200.7 / 22.8 87intergenic (+103/‑210)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (40/31);  new base G (6/6);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

TGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAG  >  NC_002947/955003‑955300
                                                                                                                                                      |                                                                                                                                                    
tgGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAACACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                      <  2:3896373/150‑1 (MQ=255)
tgGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATATTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                      <  1:1895368/150‑1 (MQ=255)
    cAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGtt                                                                                                                                                  <  1:2070262/150‑1 (MQ=255)
     aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                 <  2:4143842/150‑1 (MQ=255)
     aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                 <  2:133775/150‑1 (MQ=255)
       ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCACCGAAAAgggtgct                                                                                                                                               >  1:277785/1‑145 (MQ=255)
        cAGACCGCGGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTg                                                                                                                                              <  1:3739314/150‑1 (MQ=255)
           aCCGCTGCGGATGTGCTGATGACCGACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAg                                                                                                                                           <  1:1658699/150‑1 (MQ=255)
            ccGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAgg                                                                                                                                          <  1:5030897/150‑1 (MQ=255)
              gCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                       <  1:2042945/150‑1 (MQ=25)
               cTGCGGATGTGCTGATGACCTACGCACGCTGAACCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                                                                                                       <  2:1370335/150‑1 (MQ=255)
               cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                       <  1:3727574/150‑1 (MQ=38)
               cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                       <  2:4927226/150‑1 (MQ=38)
               cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGATTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                       <  2:3350669/150‑1 (MQ=38)
               cTGCGGACGCGCCGATGACCTCCGCACGCTGAACCTGCTGCCTGTGTCACGCACACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTCCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGGAGGAGc                                                                                                                                       <  2:3772557/150‑1 (MQ=38)
                tGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCg                                                                                                                                      <  2:2514490/150‑1 (MQ=255)
                tGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCGGTAGGAGCg                                                                                                                                      >  2:1253158/1‑150 (MQ=255)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                     <  2:1575820/150‑1 (MQ=38)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                     <  2:1905101/150‑1 (MQ=255)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                     <  1:4985418/150‑1 (MQ=255)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                     <  1:4985368/150‑1 (MQ=255)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                     <  2:733155/150‑1 (MQ=255)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                     <  1:4945882/150‑1 (MQ=255)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                     <  1:4904020/150‑1 (MQ=255)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                     <  1:4903550/150‑1 (MQ=255)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                     <  1:3791844/150‑1 (MQ=255)
                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                     <  1:379254/150‑1 (MQ=255)
                  cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                    <  2:4774467/150‑1 (MQ=38)
                    gATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGtt                                                                                                                                  <  2:142069/150‑1 (MQ=255)
                        tgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAcc                                                                                                                              <  1:4175642/150‑1 (MQ=255)
                          cTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGCCTACAACGGTCCGGGAGGCGCGGGTTTACCcg                                                                                                                            >  1:3361374/1‑150 (MQ=39)
                           tgatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCcgc                                                                                                                           <  2:2987876/150‑1 (MQ=255)
                                aCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACGGTTCCTGTAGGAGCGGGTTTACCCGCGAaga                                                                                                                      <  2:1681960/150‑1 (MQ=255)
                                   tccgcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCAGGGTGGCGGGTTTTTTTTTTGGGGGTATTTTTCCTGGGCGGGCCTCTGGGGGGGTAACCCCGCGCCTACACGGTGCCTGTTGGTGGGGGTTTACCCGCGAAGAGGc                                                                                                                   <  2:1296461/148‑1 (MQ=255)
                                       cacgcTGAAACTGCGGCCTGTTTCACGAAAACCCCGGCCAGTGTGCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGc                                                                                                               <  2:3938417/150‑1 (MQ=37)
                                            tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCGGGAGGAGCGGGGTTACCCGCGGAGAGGGCAGCACAGa                                                                                                          >  1:5076925/1‑150 (MQ=38)
                                            tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGCCGAAAAGGGTGCGGTAGGAGCGGGGTTTCCCGCGCAGAGGCCCGcacaca                                                                                                          >  2:3514616/1‑148 (MQ=38)
                                            tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCATCGCGGGTTAAACCCCGACAAAAAGGGTGCAGGAGGTGCCGGGTTAGCCCCGGAGAGGGCCGcacacg                                                                                                          >  2:1694098/1‑148 (MQ=255)
                                                            ttCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCGCCTAAAAGGTTCCGGTAGGAGCGGGTTTTACCGCGAAGAGGCCAGCACCGAAAAAAAAAAATACCTg                                                                                          >  2:1585161/1‑150 (MQ=38)
                                                             tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGGAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCGGa                                                                                         >  2:2911537/1‑150 (MQ=255)
                                                               aCGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTGGCGGGTTATACCGGTGCTTGAAGGTTGGTGTTGGTGCGGGTTTTTGCGGGTATTGGGGAGCGCAGTAAAATAAATATAGCTGttt                                                                                       >  2:4900372/1‑150 (MQ=255)
                                                                       cccGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCGACAAGGTTCCTGGAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGTTTTTaaaaaa                                                                               >  2:4466938/1‑150 (MQ=255)
                                                                       cccGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTACTACAAGGTTCCGGGAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGCAAAATAAAAATAACAGATTTTCaaaaa                                                                               >  2:3158818/1‑150 (MQ=37)
                                                                        ccGGCCACTGGGCCGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAACCCCGATCCTACAAGGTTCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa                                                                              <  2:1857297/150‑1 (MQ=38)
                                                                                     cGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCACTTCGCGGGTAAACCCCCACCGAAAAGGGTCCGGGAGGAGCCGGTTTAACCGCGAAGAGGCCAGCCCCCAAAAAAAAAAATAACAGCATTTCAAAAAATAGATCTGaaaat                                                                 >  2:232241/1‑149 (MQ=255)
                                                                                        tttttttttttGGGGGTATTTTTCCTGTGCCGCCCTCTCCGCGGTAAACCCGCCCCCAACCAGGTACCTGCCGGAGGGGTTTCACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc                                                              <  1:493559/148‑1 (MQ=255)
                                                                                                tgtGGGGGTATTTTTCGTGTGCCGGCCTCTTCGGGGTAAAACCCGCTCCAACACGGTTCCTGTCGGCGGGGTTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATaaaa                                                      <  1:3514616/148‑1 (MQ=255)
                                                                                                  tGGGGTTATTTTACCTGGGCCGGCCTATACGCGGTTAAACCCGAACCTCCTAGGTTCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                    <  1:1682387/150‑1 (MQ=38)
                                                                                                  tGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAACCCCGCTCCTACACGGTTCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                    <  1:1682563/150‑1 (MQ=255)
                                                                                                          tttttCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGCCTACAAGGTTCCGGTAGGAGCGGGGTTACCCGCGCAGAGGCCCGCCCCGAAAAATAACAATACCTGATTTTCCAAAAATAGATCTGGAAAGGTCTAATAAAACAGGAAATTa                                            >  2:4509919/1‑150 (MQ=37)
                                                                                                            tttCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCCACAACAAGGTTCCGGTAGGGGCGGGTTTAACCCCCAACAGGCCAGCACACAAAAATAAAAAAAACAGATTTTAAAAAAATATATCAGAAAAGTTATAATAAAAAAAGAAAatata                                          >  2:2306681/1‑149 (MQ=255)
                                                                                                                        gCCTCTTCGCGGGTAAAACCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCACCACACAAAAATAAAAATACCAGATTTTCAAAAAATAGAACTGAAAAGTTCTAATAAAAAAGGAAAATATACTTTAAACGCaa                              >  2:3481814/1‑149 (MQ=25)
                                                                                                                         ccTCTTCGCGGGGAAACCCGCTCCTACAAGGTTCCGGTAGGAGCGGGTTTACCCGCCAAGAGGCCACCACCGAAAAATAAAAATACCTGATTTTCCAAAAAAAGATCAGAAAAGTTCTAATAAAACAGGAAAGTATCCATTACAGGCagc                             >  2:3570135/1‑148 (MQ=37)
                                                                                                                           tctTCGCGGGTTAAACCCCTCCTACCAGGTTCCGGGAGGAGCGGGTTTACCCGCGCAGAGGCCACCACAGAAAAATAAAAATACCTGATTTTGAAAAAATAGATCTGTAAAGTTCTAATAAAAAAGGAAATTATCCTTTAAAGGCCTccg                           >  2:2149008/1‑149 (MQ=255)
                                                                                                                                  gggTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAATACAGGAAATTATCCTTTACAGGCATCCTGTTTCgg                    >  2:4259579/1‑150 (MQ=255)
                                                                                                                                    gTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCa                  >  2:1792334/1‑150 (MQ=255)
                                                                                                                                        aCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAt              >  2:2455967/1‑150 (MQ=255)
                                                                                                                                         cccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc             >  1:2538117/1‑150 (MQ=38)
                                                                                                                                         cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc             >  1:4220629/1‑150 (MQ=255)
                                                                                                                                         cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc             >  1:4220623/1‑150 (MQ=255)
                                                                                                                                         cccGCTCCTACAAGGTTCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGTAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTGTCGGCACtctc             >  2:2970779/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct            >  2:2676055/1‑150 (MQ=38)
                                                                                                                                          ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct            >  2:2825586/1‑150 (MQ=38)
                                                                                                                                          ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct            >  1:2554150/1‑150 (MQ=38)
                                                                                                                                          ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct            >  1:1716120/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct            >  1:693247/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct            >  1:2059357/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct            >  1:4633831/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct            >  1:3809211/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct            >  1:2413023/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct            >  1:3809204/1‑150 (MQ=255)
                                                                                                                                           cGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGTCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGTAAATTATCCGTTACAGGCATCCTGTTTCGGCACTAtctt           >  2:2620164/1‑150 (MQ=37)
                                                                                                                                           cGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtctt           <  2:1716120/150‑1 (MQ=255)
                                                                                                                                            gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCTCTAtcttc          >  1:4690376/1‑150 (MQ=38)
                                                                                                                                               tCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTg        >  2:663722/1‑150 (MQ=255)
                                                                                                                                                  tACAAGGTTCCTGTAGGAGGGGTTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt     <  1:1505671/150‑1 (MQ=255)
                                                                                                                                                  tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGaaaa                                                                  >  2:259481/1‑89 (MQ=255)
                                                                                                                                                  tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGaaaa                                                                  <  1:259481/89‑1 (MQ=255)
                                                                                                                                                  tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt     >  2:1909286/1‑150 (MQ=255)
                                                                                                                                                  tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt     >  1:2019522/1‑150 (MQ=255)
                                                                                                                                                  tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt     >  2:919458/1‑150 (MQ=255)
                                                                                                                                                  tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTACAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt     >  1:4337191/1‑150 (MQ=255)
                                                                                                                                                   aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTc    <  2:3361374/150‑1 (MQ=255)
                                                                                                                                                   aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTc    <  2:2113430/150‑1 (MQ=255)
                                                                                                                                                   aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTc    >  2:1818146/1‑150 (MQ=255)
                                                                                                                                                   aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTc    >  2:5023318/1‑150 (MQ=255)
                                                                                                                                                    cAAGGTTCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCa   <  1:1253158/150‑1 (MQ=255)
                                                                                                                                                    cAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCa   >  1:604120/1‑150 (MQ=255)
                                                                                                                                                    cAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCa   >  1:3643930/1‑150 (MQ=255)
                                                                                                                                                     aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg  >  1:699870/1‑150 (MQ=255)
                                                                                                                                                     aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg  >  1:3489271/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: