Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,330,402 Δ1 bp 7.1% intergenic (+21/+15) sohB → / ← yciN inner membrane protein, S49 peptidase family protein/DUF2498 protein YciN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,330,4020A. 7.1% 344.3 / 18.2 113intergenic (+21/+15)sohB/yciNinner membrane protein, S49 peptidase family protein/DUF2498 protein YciN
Reads supporting (aligned to +/- strand):  ref base A (47/57);  new base . (3/5);  total (51/62)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 113 total aligned reads displayed.

GTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCACG  >  NC_000913/1330254‑1330546
                                                                                                                                                     |                                                                                                                                                
gTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGAc                                                                                                                                                 <  1:620061/150‑1 (MQ=255)
    aGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTcg                                                                                                                                             <  2:1599940/150‑1 (MQ=255)
    aGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTcg                                                                                                                                             <  1:737949/150‑1 (MQ=255)
    aGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTcg                                                                                                                                             <  1:1369029/150‑1 (MQ=255)
     gTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTcgc                                                                                                                                            >  2:620761/1‑150 (MQ=255)
     gTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTcgc                                                                                                                                            <  2:1248345/150‑1 (MQ=255)
      tGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTcgcg                                                                                                                                           <  1:715555/150‑1 (MQ=255)
      gggTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGTTCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTcgcg                                                                                                                                           <  2:1048232/149‑1 (MQ=255)
       ggTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGt                                                                                                                                          >  1:856459/1‑150 (MQ=255)
       ggTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGt                                                                                                                                          >  1:599822/1‑150 (MQ=255)
        ggcAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGTCTTAGCCTCGCGtt                                                                                                                                         <  1:955923/148‑1 (MQ=255)
        gTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGTCTTAGCCTCGCGtt                                                                                                                                         <  1:1287702/150‑1 (MQ=255)
          cAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGTCTTAGCCTCGCGtttg                                                                                                                                       <  1:904962/150‑2 (MQ=255)
          cAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGtttt                                                                                                                                       >  2:959694/1‑150 (MQ=255)
                aCGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTa                                                                                                                                 <  2:673029/150‑1 (MQ=255)
                   cTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATc                                                                                                                              <  2:359134/150‑1 (MQ=255)
                   cTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGGCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATc                                                                                                                              <  1:837658/150‑1 (MQ=255)
                    tataTGCAGCGTAAACGACTCATTGACCGTTTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTAGGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGCAAAACGCGAGGTCTTAGCCTCGCGTTTGTCttttacat                                                                                                                             <  1:469387/150‑6 (MQ=37)
                    tataTGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGTGGCTAAGACCTCGCGTTTTGCTTTAATCa                                                                                                                             >  1:996833/1‑150 (MQ=255)
                    tataTGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGTGGCTAAGACCTCGCGTTTTGCTTTAATCa                                                                                                                             >  1:1003320/1‑150 (MQ=255)
                     atatGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGTCTTAGCCTCGCGTTTGTCttttacatc                                                                                                                            <  1:1406626/150‑7 (MQ=38)
                       atGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAAcc                                                                                                                          <  1:1242854/150‑1 (MQ=255)
                             cGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAgatga                                                                                                                    >  2:381401/1‑150 (MQ=255)
                               tAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGata                                                                                                                  <  2:833843/150‑1 (MQ=255)
                                aaaCGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGatat                                                                                                                 >  2:55614/1‑150 (MQ=255)
                                aaaCGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGatat                                                                                                                 <  1:1108553/150‑1 (MQ=255)
                                  aCGACTCATTGACCGAATCACCGGCAGCGCGGCAGCGAGCGCCGATCGAATGTTGCTACGCTGGGGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGCGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATAttt                                                                                                               <  2:307565/150‑1 (MQ=255)
                                     aCTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTc                                                                                                            >  2:97308/1‑150 (MQ=255)
                                          ttGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGaaa                                                                                                       <  2:977407/150‑1 (MQ=255)
                                           tGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAg                                                                                                      >  2:1191017/1‑150 (MQ=255)
                                            gACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGc                                                                                                     >  1:1421890/1‑150 (MQ=255)
                                             cccGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGca                                                                                                    <  1:1559250/149‑1 (MQ=255)
                                              ccGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTATGACCTCTCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGcac                                                                                                   >  2:1072428/1‑150 (MQ=255)
                                                    cACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGc                                                                                             >  2:1239537/1‑150 (MQ=255)
                                                       cGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAg                                                                                          >  2:155233/1‑150 (MQ=255)
                                                          cAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGgtg                                                                                       >  2:1100195/1‑150 (MQ=255)
                                                          cAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGgtg                                                                                       <  1:1601941/150‑1 (MQ=255)
                                                           aGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGgtgt                                                                                      >  2:859888/1‑150 (MQ=255)
                                                            gcgcGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGtt                                                                                     <  2:599822/150‑1 (MQ=255)
                                                              gcgGCAGTGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTg                                                                                   <  2:1552623/150‑1 (MQ=255)
                                                              gcgGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTg                                                                                   <  2:856459/150‑1 (MQ=255)
                                                               cgGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAGCCATTGATGTAAACGACAAACGCGAGGCTAAGACCTCGCGTGTTGCTTTAATCAACCAGATGATAGTTTTCTGAAAGCACAAGGGCCAGGTGTTTgca                                                                                 >  1:320283/1‑148 (MQ=255)
                                                                 gCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAAc                                                                                <  2:779857/150‑1 (MQ=255)
                                                                 gCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAAc                                                                                <  2:1639612/150‑1 (MQ=255)
                                                                 gCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAAc                                                                                <  1:417634/150‑1 (MQ=255)
                                                                  cAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACa                                                                               >  1:1299595/1‑150 (MQ=255)
                                                                     agagCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACatat                                                                            <  1:1430605/150‑1 (MQ=255)
                                                                     agagCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACatat                                                                            >  2:1100059/1‑150 (MQ=255)
                                                                          gcCGATCGATCGCGGCTACGCTGGTGGCAGCGGGGGCCAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAac                                                                       <  2:320210/150‑1 (MQ=255)
                                                                           ccGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAaca                                                                      >  1:42327/1‑150 (MQ=255)
                                                                              atcgatTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACcg                                                                   >  1:225281/1‑150 (MQ=255)
                                                                                      ttgCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCtt                                                           <  1:1119767/150‑1 (MQ=255)
                                                                                       tgCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCttt                                                          >  2:1099077/1‑150 (MQ=255)
                                                                                        gCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCtttt                                                         >  1:1064930/1‑150 (MQ=255)
                                                                                              cTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGt                                                   <  1:1024025/150‑1 (MQ=255)
                                                                                                ggtggCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTg                                                 <  2:1230112/150‑1 (MQ=255)
                                                                                                  tggCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGc                                               <  2:658410/150‑1 (MQ=255)
                                                                                                  tgccAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACata                                                                             <  2:773862/117‑1 (MQ=255)
                                                                                                  tgccAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACata                                                                             >  1:773862/4‑120 (MQ=255)
                                                                                                   ggCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCa                                              >  2:395941/1‑150 (MQ=255)
                                                                                                       gCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATcc                                          <  1:644838/150‑1 (MQ=255)
                                                                                                        cGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCt                                         <  2:554580/150‑1 (MQ=255)
                                                                                                        cGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCt                                         <  1:1296659/150‑1 (MQ=255)
                                                                                                          gggTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTg                                       >  1:899207/1‑150 (MQ=255)
                                                                                                           ggTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGt                                      <  1:1679524/150‑1 (MQ=255)
                                                                                                            gTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGtt                                     >  2:1141953/1‑150 (MQ=255)
                                                                                                             tCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTc                                    <  2:360364/150‑1 (MQ=255)
                                                                                                              cAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCa                                   <  1:1500000/150‑1 (MQ=255)
                                                                                                              cAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTATTGCTTTAATCAACCAGATGATAATATTCTGAAAGCACATGGGACAGCTGTTTGAAAATAATAAACACCGAGGTTCTTTAGGCTGTTGGCAATCgttgttaa                                   >  1:719407/1‑148 (MQ=255)
                                                                                                                aaaGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATc                                 <  1:48628/150‑1 (MQ=255)
                                                                                                                aaaGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATc                                 <  1:381401/150‑1 (MQ=255)
                                                                                                                aaaGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATc                                 <  1:713603/150‑1 (MQ=255)
                                                                                                                aaaGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATc                                 <  1:837982/150‑1 (MQ=255)
                                                                                                                  aGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                               >  1:776993/1‑150 (MQ=255)
                                                                                                                   gCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTaa                              >  1:739546/1‑150 (MQ=255)
                                                                                                                   gCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTaa                              >  1:829817/1‑150 (MQ=255)
                                                                                                                     cATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAAACAACCAGATGATATTTTTCTGAAAGCACATGTTACAGGTGATTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGTCAATCGGTGTTCATATaaaa                            >  2:93717/1‑150 (MQ=255)
                                                                                                                       ttGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAag                          <  1:1513750/150‑1 (MQ=255)
                                                                                                                        tGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagt                         <  1:685628/150‑1 (MQ=255)
                                                                                                                         gATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagta                        >  2:705373/1‑150 (MQ=255)
                                                                                                                          aTGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagtag                       <  1:97308/150‑1 (MQ=255)
                                                                                                                          aTGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagtag                       >  1:1596661/1‑150 (MQ=255)
                                                                                                                            gTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGTTATTTTTCTGAAAGCACATGGGCCAGGTTTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTc                     <  2:1066537/150‑1 (MQ=255)
                                                                                                                            gTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTc                     >  1:391600/1‑150 (MQ=255)
                                                                                                                            gTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTCTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagtagta                     >  2:1307745/1‑149 (MQ=255)
                                                                                                                               aaaGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGcc                  >  1:185289/1‑150 (MQ=255)
                                                                                                                               aaaGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGcc                  <  1:1601532/150‑1 (MQ=255)
                                                                                                                                aaGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAGCCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACCCCGCGGTGCTTTTGGCTGTTGGCAATCCTAGTTCATCTAAAAAGTAGTCGCCg                 >  1:330377/1‑150 (MQ=255)
                                                                                                                                  gcaaaaCGCGAGGTCTTAGCCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGt               >  2:594128/4‑150 (MQ=255)
                                                                                                                                  gcaaaaCGCGAGGTCTTAGCCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGt               >  2:1325395/4‑150 (MQ=255)
                                                                                                                                  gACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGt               <  1:155233/150‑1 (MQ=255)
                                                                                                                                   aCAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTa              <  1:1202290/150‑1 (MQ=255)
                                                                                                                                     aaaCGCGAGGTCTTAGCCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTaaa            >  2:990491/1‑150 (MQ=255)
                                                                                                                                        cgcgAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATAc         >  2:323070/1‑150 (MQ=255)
                                                                                                                                        cgcgAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATAc         <  1:1239537/150‑1 (MQ=255)
                                                                                                                                         gcgAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATAcc        <  2:517436/150‑1 (MQ=255)
                                                                                                                                         gcgAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATAcc        <  1:543719/150‑1 (MQ=255)
                                                                                                                                             ggCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTATAAAGTAGTCGCCGGTAAATACCAGCa    >  2:1339227/1‑150 (MQ=255)
                                                                                                                                             ggCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCa    >  1:1410766/1‑150 (MQ=255)
                                                                                                                                                cTAAGACCTCTCGTTTTGCTTTTATCAACCTGATGATAATTTTGTTAAAGCTCATTGGTCAGGTGTTTGTACATTTTAAACTCCGCGGTGCTTTTGGCTGTTGGCATTCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCACg  <  1:25365/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCACG  >  NC_000913/1330254‑1330546

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: