Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,177,172 T→A 5.9% intergenic (+44/‑186) tufB → / → secE translation elongation factor Tu 2/Sec translocon subunit SecE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,177,1720TA 5.9% 603.1 / 18.3 187intergenic (+44/‑186)tufB/secEtranslation elongation factor Tu 2/Sec translocon subunit SecE
Reads supporting (aligned to +/- strand):  ref base T (91/85);  new base A (6/5);  total (97/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 193 total aligned reads displayed.

TGTTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGT  >  NC_000913/4177023‑4177321
                                                                                                                                                     |                                                                                                                                                     
tgttACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAtt                                                                                                                                                       <  2:1011095/150‑1 (MQ=255)
   tACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGa                                                                                                                                                    >  2:799146/1‑150 (MQ=255)
    aCCCTGATCCACCCGAGCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGAt                                                                                                                                                   <  1:641131/150‑1 (MQ=255)
     cccTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                                                                                  >  2:1418904/1‑150 (MQ=255)
     cccTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                                                                                  >  1:1061791/1‑150 (MQ=255)
     cccTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                                                                                  >  1:875077/1‑150 (MQ=255)
       cTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGcc                                                                                                                                                >  1:480866/1‑150 (MQ=255)
           tCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttt                                                                                                                                            <  2:127763/150‑1 (MQ=255)
             cACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttttt                                                                                                                                          <  1:1156686/150‑1 (MQ=255)
             cACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGCGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttttt                                                                                                                                          >  1:1254879/1‑150 (MQ=255)
              aCCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAATTTCTGAGCTAATTTCCGATATCATTTGAAGCAATTCGCACTAAAAGTGTATAATTTTATTCGCTTTTTg                                                                                                                                         <  2:966604/150‑1 (MQ=18)
               cccGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCCTTTTGc                                                                                                                                        <  2:543545/150‑1 (MQ=255)
               cccGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagt                                                                                                                                        <  1:175170/150‑4 (MQ=255)
                ccGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                                                                       <  1:115227/150‑1 (MQ=255)
                ccGATCGCGATGGACGACGGTCTGCGGTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                                                                       <  1:1234036/150‑1 (MQ=255)
                      gcgATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCttt                                                                                                                                 <  1:25458/150‑1 (MQ=255)
                        gATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                                                                               <  1:1192765/150‑1 (MQ=255)
                        gATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                                                                               <  2:1237024/150‑1 (MQ=255)
                        gATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                                                                               <  2:1099470/150‑1 (MQ=255)
                           ggACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAc                                                                                                                            <  2:733420/150‑1 (MQ=255)
                            gacgacGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAcc                                                                                                                           >  1:1003900/1‑150 (MQ=255)
                             acgacgGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCa                                                                                                                          <  1:379133/150‑1 (MQ=255)
                                  ggTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAAcc                                                                                                                     <  1:881092/150‑1 (MQ=255)
                                    tCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTg                                                                                                                   >  1:1417264/1‑150 (MQ=255)
                                           ttCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtca                                                                                                              >  2:576414/1‑148 (MQ=255)
                                            tCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcatca                                                                                                           >  2:1153512/1‑150 (MQ=255)
                                             cGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAg                                                                                                          <  2:1417766/150‑1 (MQ=255)
                                             cGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAg                                                                                                          >  1:46683/1‑150 (MQ=255)
                                             cGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAg                                                                                                          <  2:647450/150‑1 (MQ=255)
                                             cGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAg                                                                                                          <  1:1239912/150‑1 (MQ=255)
                                             cGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCGCATCAg                                                                                                          <  2:1061791/150‑1 (MQ=255)
                                              gCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGt                                                                                                         >  2:568030/1‑150 (MQ=255)
                                               cAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTg                                                                                                        <  2:235174/150‑1 (MQ=255)
                                               cAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGAGCACTAAAAGGGCATCATTTGATGAACTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTg                                                                                                        >  1:1178786/1‑150 (MQ=255)
                                                  tCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAtt                                                                                                     <  1:919275/150‑1 (MQ=255)
                                                    cGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGAAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAtttt                                                                                                   >  1:1199189/1‑150 (MQ=255)
                                                       gAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCtt                                                                                                <  2:875077/150‑1 (MQ=255)
                                                       gAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCtt                                                                                                <  1:854761/150‑1 (MQ=255)
                                                        aaGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCttt                                                                                               <  2:452884/150‑1 (MQ=255)
                                                         aGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                                                                              <  2:1402452/150‑1 (MQ=255)
                                                         aGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                                                                              <  2:262629/150‑1 (MQ=255)
                                                           gcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTc                                                                                            <  1:984081/150‑1 (MQ=255)
                                                            cggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCa                                                                                           <  1:1310774/150‑1 (MQ=255)
                                                             ggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCAt                                                                                          <  2:702794/150‑1 (MQ=255)
                                                             ggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCAt                                                                                          >  1:361205/1‑150 (MQ=255)
                                                                      gTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                                                                                <  2:1331372/119‑1 (MQ=255)
                                                                      gTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                                                                                >  1:1331372/1‑119 (MQ=255)
                                                                       ttGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCt                                                                                <  1:721291/150‑1 (MQ=255)
                                                                        tGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTg                                                                               <  1:695799/150‑1 (MQ=255)
                                                                        tGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTg                                                                               <  1:460579/150‑1 (MQ=255)
                                                                          gcgcGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTgag                                                                             <  1:194539/150‑1 (MQ=255)
                                                                                    tGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTg                                                                   <  1:799146/150‑1 (MQ=255)
                                                                                    tGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTg                                                                   >  1:31433/1‑150 (MQ=255)
                                                                                     gTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTTAGACAGGCTCTGt                                                                  >  2:1030371/1‑150 (MQ=255)
                                                                                      tAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGtt                                                                 >  1:592493/1‑150 (MQ=255)
                                                                                       aGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTg                                                                >  2:1148782/1‑150 (MQ=255)
                                                                                        gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa                                                               >  1:1012974/1‑150 (MQ=255)
                                                                                        gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa                                                               >  1:704595/1‑150 (MQ=255)
                                                                                          aaaaGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAgg                                                             <  2:468846/150‑1 (MQ=255)
                                                                                          aaaaGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAgg                                                             <  2:1375319/150‑1 (MQ=255)
                                                                                               ttCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGta                                                        >  1:1133800/1‑150 (MQ=255)
                                                                                                   gAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAAt                                                    <  2:185503/150‑1 (MQ=255)
                                                                                                         aTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGaaa                                              >  1:1293633/1‑150 (MQ=255)
                                                                                                          atGCCGATAACATTTCACGCCATGCGCACTAAAAGGGCATCATTTGATGCACTTTTTGCACGCTATCGTAACAGCACCTGGCTCATCAGAGATTTTCTTTGTCATAATCATTGCTGAGAAAGGCTCACTTGAGGGCGGATTATCCCaaaa                                             >  2:486208/2‑150 (MQ=255)
                                                                                                           tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAg                                            >  1:745384/1‑150 (MQ=255)
                                                                                                           tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAg                                            >  1:76268/1‑150 (MQ=255)
                                                                                                            gCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGc                                           >  1:1356798/1‑150 (MQ=255)
                                                                                                            gCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGc                                           >  1:1043564/1‑150 (MQ=255)
                                                                                                              cGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTa                                         >  1:814871/1‑150 (MQ=255)
                                                                                                                aTAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAAt                                       <  2:1291500/150‑1 (MQ=255)
                                                                                                                 tAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATa                                      >  1:706329/1‑150 (MQ=255)
                                                                                                                   aCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATAcg                                    >  1:1201303/1‑150 (MQ=255)
                                                                                                                    cATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACa                                                                          >  1:159288/1‑111 (MQ=255)
                                                                                                                     aTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATAcgcg                                  >  1:610363/1‑150 (MQ=255)
                                                                                                                      tttGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCa                                                                                           >  1:904767/1‑92 (MQ=255)
                                                                                                                         gACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTc                              <  2:1193489/150‑1 (MQ=255)
                                                                                                                                 gcgtgcaaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCGCATCAGGGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATACGAAAAGCTAATACGCGTTTCGATTTGGt                      >  1:1268516/8‑150 (MQ=255)
                                                                                                                                 gcgcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGt                      <  2:745384/150‑1 (MQ=255)
                                                                                                                                  cgcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGtt                     >  2:840643/1‑150 (MQ=255)
                                                                                                                                  cgcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTACAATCCGAAAAGCTAATACGCGTTTCGATTTGGtt                     <  2:1113891/150‑1 (MQ=255)
                                                                                                                                    tgcaaaaaGGGCATCAAATGATGCCCTTTTAGTGCGCATTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTg                   >  2:787462/5‑150 (MQ=38)
                                                                                                                                    cACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTg                   >  2:918119/1‑150 (MQ=255)
                                                                                                                                    cACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTg                   <  2:162007/150‑1 (MQ=255)
                                                                                                                                       tAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCt                <  1:66465/150‑1 (MQ=255)
                                                                                                                                       tAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCt                <  1:945458/150‑1 (MQ=255)
                                                                                                                                        aaaaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTc               >  2:1133504/1‑150 (MQ=255)
                                                                                                                                        aaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTTGTTTGCCTc               >  2:1384898/1‑150 (MQ=255)
                                                                                                                                         aaaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTcg              >  1:980539/1‑150 (MQ=255)
                                                                                                                                          aaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTcgc             >  2:945452/1‑150 (MQ=255)
                                                                                                                                            gggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGa           >  1:1227439/1‑150 (MQ=255)
                                                                                                                                             ggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGAt          <  1:1418904/150‑1 (MQ=255)
                                                                                                                                               catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGt                                                         <  2:388991/101‑1 (MQ=255)
                                                                                                                                               catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGt                                                         >  1:388991/1‑101 (MQ=255)
                                                                                                                                               catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATcg        >  1:32733/1‑150 (MQ=255)
                                                                                                                                               catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATcg        <  1:28766/150‑1 (MQ=255)
                                                                                                                                                 tcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATcgcg      >  1:1385838/1‑150 (MQ=255)
                                                                                                                                                   atTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCggg    <  2:592493/150‑1 (MQ=255)
                                                                                                                                                   atTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCggg    <  2:1238997/150‑1 (MQ=255)
                                                                                                                                                    tttGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGATCATCAGTGATTTTCTTTGTCATAATCATTTCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCTATTTGGTTTGCATGGCGATcgcgttt   <  2:900879/150‑4 (MQ=255)
                                                                                                                                                     ttGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGt  <  1:1111654/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGT  >  NC_000913/4177023‑4177321

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