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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,177,1720TA 8.8% 356.9 / 19.5 125intergenic (+44/‑186)tufB/secEtranslation elongation factor EF‑Tu 2/preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (52/62);  new base A (4/7);  total (56/69)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

Only 100 of 128 total aligned reads displayed.

GTTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGT  >  NC_000913/4177024‑4177321
                                                                                                                                                    |                                                                                                                                                     
gttACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATGCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAatt                                                                                                                                                      <  2:1160363/150‑1 (MQ=255)
  tACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGCCGCAATGCGCACTAAAAGGGCATCATTTGa                                                                                                                                                    >  1:349131/1‑150 (MQ=255)
  tACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGa                                                                                                                                                    >  2:191894/1‑150 (MQ=255)
   aCCCTGATACACCCGATCGCGATGGACGACGGTCTGCGTTTCGGAATCCTTGAAGGCGGCCGTAACGTTTGCGCGGGCGTTGTATCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGAt                                                                                                                                                   <  1:1304268/150‑1 (MQ=21)
    cccTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                                                                                  <  1:61930/150‑1 (MQ=255)
    cccTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                                                                                  >  2:1955415/1‑150 (MQ=255)
    cccTGATCCACACGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGGTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                                                                                  >  1:1212541/1‑150 (MQ=255)
      cTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGcc                                                                                                                                                >  1:1472276/1‑150 (MQ=255)
          tCCACCCGATCGCGATGGAGGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttt                                                                                                                                            <  1:576564/150‑1 (MQ=255)
           ccACCCGATCGCGATTGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCTTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttt                                                                                                                                           <  2:1726474/150‑1 (MQ=255)
           ccACCCGATCGCGATGTACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttt                                                                                                                                           <  1:1366780/150‑1 (MQ=255)
           ccACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGGGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttt                                                                                                                                           <  1:272487/150‑1 (MQ=255)
            cACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttttt                                                                                                                                          >  2:1758591/1‑150 (MQ=255)
            cACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttttt                                                                                                                                          >  2:611733/1‑150 (MQ=255)
            cACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttttt                                                                                                                                          <  2:962972/150‑1 (MQ=255)
             tcccGATCGCGAAGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGTCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttag                                                                                                                                         <  1:725982/149‑3 (MQ=25)
             ccccGATCTCGATGGACGACGGTCTGCGTTTCGCAAACCGTGAAGGCTGTCGTACCGTTTGCGCGGGCGTTGTAGCAAAATTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttag                                                                                                                                         <  1:335261/149‑3 (MQ=25)
             ccccGATCGCGATGGACGACGGGCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTg                                                                                                                                         <  2:1648927/149‑1 (MQ=255)
             aCCCGATCGCGATTGACGACGGTCTGAGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTg                                                                                                                                         <  2:1084676/150‑1 (MQ=255)
              cccGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGc                                                                                                                                        <  2:1836089/150‑1 (MQ=255)
               ccGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                                                                       >  1:1497731/1‑150 (MQ=255)
               ccGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTGCTGAGCTAATTGCCGATAACATTTGAAGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                                                                       <  2:879595/150‑1 (MQ=255)
               ccGAGCGCGATGGACGACGGGCTGCGTTGCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                                                                       <  1:1998509/150‑1 (MQ=255)
                 gATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGGTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACg                                                                                                                                     <  1:833514/150‑1 (MQ=255)
                  aTCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAGGTTCTGAGCTAAGTGCCGATAACATTTGACGCAATGCGCCCTAAAAGGGCATCATTTGATGCCCTTGTTGCACGc                                                                                                                                    >  1:1005034/1‑150 (MQ=255)
                    cgcgATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTAGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCtt                                                                                                                                  <  1:987969/150‑1 (MQ=255)
                      cgATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTc                                                                                                                                >  2:828661/1‑150 (MQ=255)
                       gtttgacgacGGGCTGCGTTTCGCAATGCGTGAAGGCGGCCGGACCGTTGGCGCGGGCGTTGTAGCAACAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                                                                               <  2:763737/146‑1 (MQ=37)
                       gATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGTGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                                                                               <  1:1485864/150‑1 (MQ=255)
                       gATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                                                                               >  2:1872044/1‑150 (MQ=255)
                        aTGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGt                                                                                                                              <  1:1990190/150‑1 (MQ=255)
                           gacgacGGTCTGCGTATCGCAATCCGTGAAGGCGGACGTACAGTTGGCGCGGCCGTTGTAGCAAAAGTTCTGAGCTATTTGCCGATTACATATGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCATTTTGCACGCTTTCGTAcc                                                                                                                           >  1:1927627/1‑150 (MQ=37)
                            ccgacgGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGAGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCa                                                                                                                          <  2:95179/149‑1 (MQ=255)
                            acgacgGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGGTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCAGCATTTGATGCCCTTTTTGCACGCTTTCGTACCa                                                                                                                          <  1:1216044/150‑1 (MQ=255)
                            acgacgGTCTGCGATTCGCAATCCGTGAAGGCGGCCGGACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAACGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCa                                                                                                                          <  1:1215938/150‑1 (MQ=255)
                             cgacgGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAg                                                                                                                         <  2:111106/150‑1 (MQ=255)
                                 ggTCTGCGTTTCGCGATCCTTGAATGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAATTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAAcc                                                                                                                     <  1:972615/150‑1 (MQ=255)
                                     tGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                                                                                 <  2:751508/150‑1 (MQ=255)
                                     tGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                                                                                 <  2:1448559/150‑1 (MQ=255)
                                      gCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                                                                                >  1:1867938/1‑150 (MQ=255)
                                       cGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGGTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtc                                                                                                               >  1:1767216/1‑150 (MQ=255)
                                          ttCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTGAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcatc                                                                                                            <  1:406245/150‑1 (MQ=255)
                                           tCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcatca                                                                                                           >  1:8699/1‑150 (MQ=255)
                                            cGCAATCCGTGAAGTCGGCGGTACCGTTGGAGCGGGCGTTGTAGCAAAAGTGCTGAGATATTTGCCGATACCATTTGACGCAATGCGCACTAACAGGGCATCATTTGATGCCGTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAg                                                                                                          <  1:334573/150‑1 (MQ=255)
                                            cGCAATCCGTGAAGGCGTCCGTACCGTTGGAGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCTATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGCCCCTGGCtcagcat                                                                                                          >  2:1421255/1‑149 (MQ=255)
                                              cAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTg                                                                                                        >  1:405431/1‑150 (MQ=255)
                                                     tGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCt                                                                                                 <  1:711897/150‑1 (MQ=255)
                                                     tGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCACCAGTGATTTTCt                                                                                                 <  1:557829/150‑1 (MQ=255)
                                                       aaGGCGGCCGTACCGTTGGCGCGGGCATTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCttt                                                                                               <  2:1436986/150‑1 (MQ=255)
                                                         ggcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGACCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTATTTGt                                                                                             >  2:1239284/1‑150 (MQ=255)
                                                         ggcggcCGTACCGGTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGTGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGt                                                                                             <  2:1296249/150‑1 (MQ=255)
                                                           cggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCa                                                                                           <  2:1563929/150‑1 (MQ=255)
                                                              ccGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCGTGTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATaa                                                                                        >  2:321761/1‑150 (MQ=255)
                                                                 tACCGTTTGCGCGGGCGGTGTAGCATAAGTTCTGAGCTAATTGCCGATAACTTTTGACGGAATGCGAACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCTGAACCTGTCTCATCTGTTATTTTCTTTTTCATAATCa                                                                                     <  1:1938041/150‑1 (MQ=255)
                                                                    cGTTTGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCAGattgattg                                                                                  >  1:1687215/1‑150 (MQ=255)
                                                                    cGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTg                                                                                  >  2:1232964/1‑150 (MQ=255)
                                                                       tGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTg                                                                               <  2:722568/150‑1 (MQ=255)
                                                                         gcgcGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTgag                                                                             <  1:12099/150‑1 (MQ=255)
                                                                          cgcgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTgaga                                                                            >  2:1836857/1‑150 (MQ=255)
                                                                           gcgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGAc                                                                           <  2:211357/150‑1 (MQ=255)
                                                                                   tGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTg                                                                   <  1:1619712/150‑1 (MQ=255)
                                                                                    gTAGCAAAAGTTCTGAGCTGATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGt                                                                  <  1:528009/150‑1 (MQ=255)
                                                                                      aGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTg                                                                <  2:210808/150‑1 (MQ=255)
                                                                                      aGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACATGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTg                                                                >  1:1422669/1‑150 (MQ=255)
                                                                                       gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa                                                               >  2:1067028/1‑150 (MQ=255)
                                                                                       gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa                                                               >  1:20053/1‑150 (MQ=255)
                                                                                       gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa                                                               >  1:1494960/1‑150 (MQ=255)
                                                                                        cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAg                                                              >  2:1467439/1‑150 (MQ=255)
                                                                                         aaaaGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAgg                                                             <  2:1679753/150‑1 (MQ=255)
                                                                                           aaGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGc                                                           <  2:357859/150‑1 (MQ=255)
                                                                                                 tGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGACCCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATaa                                                     >  2:1695791/1‑150 (MQ=255)
                                                                                                 tGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa                                                               >  1:1859407/1‑140 (MQ=255)
                                                                                                 tGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa                                                               <  2:1859407/140‑1 (MQ=255)
                                                                                                  gAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAAt                                                    >  2:1818213/1‑150 (MQ=255)
                                                                                                    gCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATcc                                                  <  1:1404300/150‑1 (MQ=255)
                                                                                                    gCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATcc                                                  <  1:1290602/150‑1 (MQ=255)
                                                                                                     cTAATTGCCGATAACATTTGACGCAAAGCGCACAAAAATGGCATCATTTGATGCGCTATTTGCACGCTTTCGTACCAGCACCTGGCTCATCAGATATATTCATTGTCAAAATCATAGCTTAGAAAGCCTCTGTTGAGGGCCTATAATCCg                                                 >  2:1199199/1‑150 (MQ=255)
                                                                                                      tAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGa                                                <  1:717882/150‑1 (MQ=255)
                                                                                                        aTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGaaa                                              >  1:11422/1‑150 (MQ=255)
                                                                                                        aTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATATTCTTAGTCATAATCATTGCAGAGACAGGCTCTGAAGAGGGCGTATAATCCTaaa                                              >  1:227378/1‑150 (MQ=255)
                                                                                                           gCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGc                                           <  1:1817586/150‑1 (MQ=255)
                                                                                                            ccGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCt                                          <  1:1615973/150‑1 (MQ=255)
                                                                                                              gATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCGTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTaa                                        <  1:1759291/150‑1 (MQ=255)
                                                                                                                tAACATTTGACGCAATGAGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATa                                      <  1:1445387/150‑1 (MQ=255)
                                                                                                                  aCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATAcg                                    >  2:259649/1‑150 (MQ=255)
                                                                                                                        gACGCAATGCGCACTAACAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTc                              <  2:618967/150‑1 (MQ=255)
                                                                                                                        gACGCAATGCGCACTAAAAGTGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTTGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTc                              <  2:610718/150‑1 (MQ=255)
                                                                                                                        gACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGCCAGGCGCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTc                              >  1:358522/1‑150 (MQ=255)
                                                                                                                          cGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGCGTATAATCCGAAAAGCTAATACGCGTTTCGAt                           <  1:1967963/150‑1 (MQ=255)
                                                                                                                           gCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCCCGCTTTCGTACCAGAACCTGGCTCTTCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGAt                           >  2:223185/1‑150 (MQ=255)
                                                                                                                              aTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTg                        <  2:11890/150‑1 (MQ=255)
                                                                                                                              aTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGAGATTTTCTTTGTCATAATCATAGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTg                        <  1:108084/150‑1 (MQ=255)
                                                                                                                                      aaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCt                >  2:1355883/2‑150 (MQ=255)
                                                                                                                                       aaaaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTc               >  2:386276/1‑150 (MQ=255)
                                                                                                                                        aaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTGTCGATTTGGTTTGACTcg              >  1:1647454/1‑150 (MQ=255)
                                                                                                                                        aaaGGGCATCAAATGATGCCCTTTTTGCACGCGTTCGTACCAGAACCTGGCTCATCAGTGATTTTCGTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATACGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTcg              >  1:1763956/1‑150 (MQ=255)
                                                                                                                                           gggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGa           >  2:2007583/1‑150 (MQ=255)
                                                                                                                                                 catTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGGGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCgg     <  1:1758591/150‑1 (MQ=255)
                                                                                                                                                   tttGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCgggg   <  1:1139889/150‑1 (MQ=255)
                                                                                                                                                    ttGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGt  <  2:1557087/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGT  >  NC_000913/4177024‑4177321

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