Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,089,854 A→T 16.2% intergenic (‑141/‑142) yefM ← / → hisL antitoxin of the YoeB‑YefM toxin‑antitoxin system/his operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,089,8540AT16.2% 245.2 / 39.9 105intergenic (‑141/‑142)yefM/hisLantitoxin of the YoeB‑YefM toxin‑antitoxin system/his operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (47/41);  new base T (7/10);  total (54/51)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

Only 100 of 106 total aligned reads displayed.

TAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTA  >  NC_000913/2089758‑2089953
                                                                                                |                                                                                                   
taattaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTc                                                                                                  >  2:534406/1‑100 (MQ=255)
 aattaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTcc                                                                                                 <  2:1084806/100‑1 (MQ=255)
  attaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCct                                                                                                >  1:1046282/1‑100 (MQ=255)
    taaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCctct                                                                                              >  1:2663917/1‑100 (MQ=255)
        tatttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGc                                                                                          <  1:1036031/100‑1 (MQ=255)
        tatttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGc                                                                                          <  1:2772472/100‑1 (MQ=255)
        tatttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGc                                                                                          <  2:802534/100‑1 (MQ=255)
        tatttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGc                                                                                          <  2:797023/100‑1 (MQ=255)
        tatttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGc                                                                                          >  2:389456/1‑100 (MQ=255)
         atttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACTAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGct                                                                                         <  1:420312/100‑2 (MQ=255)
         atttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGct                                                                                         <  1:2442059/100‑2 (MQ=255)
         atttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAgtc                                                                                         <  2:1965910/100‑3 (MQ=255)
          tttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGccc                                                                                        <  1:331394/100‑1 (MQ=255)
          tttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGccc                                                                                        <  2:1521477/100‑1 (MQ=255)
          tttattAATTAGTTTGTAGATCAAGGTATTGTCAGCGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGccc                                                                                        <  2:888584/100‑1 (MQ=255)
           ttattaATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGcccc                                                                                       >  2:1262577/1‑100 (MQ=255)
            tattaATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGccccc                                                                                      >  1:1309862/1‑100 (MQ=255)
            tattaATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGccccc                                                                                      <  2:1091232/100‑1 (MQ=255)
              ttaattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCcctc                                                                                    >  2:2915499/1‑100 (MQ=255)
               taattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGctgatggc                                                                                   <  2:1767100/100‑8 (MQ=38)
                 attaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGctgatggcac                                                                                 <  2:2454462/100‑10 (MQ=38)
                     gTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTccc                                                                             >  1:1050754/1‑100 (MQ=255)
                     gTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCCTCGCTATTTTTGTTTTTTTTTTCTAAGAGGACAGTCCTCTTAGCCCCCTCTTTTccc                                                                             >  1:2716961/1‑100 (MQ=255)
                      tttGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTcccc                                                                            >  2:668896/1‑100 (MQ=255)
                      tttGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCACCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTcccc                                                                            <  2:2052800/100‑1 (MQ=255)
                      ttccgatctCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTcccc                                                                            <  2:877708/92‑1 (MQ=255)
                       ttGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTGTTTTTGTTGTTATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTccctg                                                                           >  1:1617744/1‑98 (MQ=255)
                        tGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGc                                                                          >  2:2762106/1‑100 (MQ=255)
                        ggTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGc                                                                          >  1:859875/2‑100 (MQ=255)
                         gTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCt                                                                         <  1:1530890/100‑1 (MQ=255)
                           aGATCAAGGTATTGTCAGTGAGACGAAACTCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGCAGACAGTAAGCTTCGCCCACTCCTTTCCCCGCtaa                                                                       >  1:712426/1‑98 (MQ=255)
                           aGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGCCAGTCCTATTAGCCCCCTCCTTTCCCCGCtca                                                                       >  1:1018635/1‑100 (MQ=255)
                           aGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtca                                                                       >  2:562208/1‑100 (MQ=255)
                           aGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtca                                                                       <  1:1446472/100‑1 (MQ=255)
                            gATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcat                                                                      <  1:305970/100‑1 (MQ=255)
                            gATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTCGTGCCTTCAGATAAGAGAAAAGTCCTCTTATCCCCCTCCTTTACCCGCtca                                                                       >  1:958669/1‑99 (MQ=255)
                              tCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcattc                                                                    <  2:138871/100‑1 (MQ=255)
                               cAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTccccagatcggaa                                                                   >  1:877705/1‑91 (MQ=255)
                               cAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcattca                                                                   >  1:625084/1‑100 (MQ=255)
                                 aGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCACTcattcatt                                                                 <  2:912551/100‑1 (MQ=255)
                                   gTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTaa                                                               <  2:204156/100‑1 (MQ=255)
                                      ttGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTCAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACa                                                            >  2:2841359/1‑100 (MQ=255)
                                       tGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTTTTTTTTTTTTCTAAGAGGACAGTCCTCTTAGCCCTCTCCTTTCCCCGCTCATTCATTaa                                                               >  1:155928/1‑96 (MQ=255)
                                        gTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTCGCCCCCTCCTTTCCCCGCTCATTCATTAACCaaa                                                          >  1:107118/1‑100 (MQ=255)
                                        gTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACaaa                                                          >  1:1375185/1‑100 (MQ=255)
                                        gTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACaaa                                                          >  1:1459279/1‑100 (MQ=255)
                                          cAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATc                                                        <  1:873997/100‑1 (MQ=255)
                                          cAGAGAGACGAAAATCCAGGATTCGCAATATTTGGTGCCATCAGCTAAGAGGAGAGGCCTCGTAGCCCCGTCCGTTTTCCGCTCATTCATTAAACAAAt                                                         >  1:240761/1‑99 (MQ=255)
                                           aGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATcc                                                       >  1:550668/1‑100 (MQ=255)
                                               agaCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTg                                                   <  2:1884030/100‑1 (MQ=255)
                                               agaCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTg                                                   <  2:2322271/100‑1 (MQ=255)
                                                gaCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGc                                                  >  1:483273/1‑100 (MQ=255)
                                                gaCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGc                                                  >  2:347326/1‑100 (MQ=255)
                                                  cGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCa                                                <  1:1453767/100‑1 (MQ=255)
                                                  cGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCa                                                >  1:597454/1‑100 (MQ=255)
                                                    aaaaTCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATa                                              <  2:517550/100‑1 (MQ=255)
                                                      aaTCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATaaa                                            >  1:362720/1‑100 (MQ=255)
                                                        tCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAat                                          <  1:617998/100‑1 (MQ=255)
                                                         ccAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAata                                         <  1:2881463/100‑1 (MQ=255)
                                                         ccAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAata                                         <  2:1459285/100‑1 (MQ=255)
                                                           aGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAatata                                       <  2:545926/100‑1 (MQ=255)
                                                              cTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATaa                                    <  1:2005137/100‑1 (MQ=255)
                                                                  gCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATaaaaaa                                >  1:2726608/1‑100 (MQ=255)
                                                                  gCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATaaaaaa                                >  1:654733/1‑100 (MQ=255)
                                                                    tATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGc                              <  2:2278489/100‑1 (MQ=255)
                                                                    tATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCACCTCCTTTCACCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGc                              <  1:820963/100‑1 (MQ=255)
                                                                     aTTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGcc                             >  2:237738/1‑100 (MQ=255)
                                                                      tttttGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAg                               >  2:134658/1‑97 (MQ=255)
                                                                      tttttGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGccc                            >  1:2605522/1‑100 (MQ=255)
                                                                       ttttGGTGCCATCAGCTAAGAGGACAGTCTTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACGAATCCATTGCCATAAAATATATAAAAAAGCCCt                           <  1:1756052/100‑1 (MQ=255)
                                                                        tttGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCtt                          >  1:1486881/1‑100 (MQ=255)
                                                                           ggTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCt                       >  1:471673/1‑100 (MQ=255)
                                                                           ggTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCt                       >  1:1175462/1‑100 (MQ=255)
                                                                           ggTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCt                       <  2:360123/100‑1 (MQ=255)
                                                                           ggTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCt                       <  2:1327868/100‑1 (MQ=255)
                                                                              gCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTc                    >  2:2786568/1‑100 (MQ=255)
                                                                              gCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTc                    >  2:795358/1‑100 (MQ=255)
                                                                              gCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTc                    >  1:228189/1‑100 (MQ=255)
                                                                              gCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTc                    <  2:2589471/100‑1 (MQ=255)
                                                                               ccATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCt                   <  1:2674854/100‑1 (MQ=255)
                                                                                 gggggCTAAGAGGACTGTCCTCTTAGCTCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTaa                 >  2:1186961/5‑100 (MQ=255)
                                                                                   catcTAAGAGGACAGCCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACg               <  1:1189288/97‑1 (MQ=255)
                                                                                   cAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACg               <  2:1309867/100‑1 (MQ=255)
                                                                                   cAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACg               <  2:2696788/100‑1 (MQ=255)
                                                                                    atcTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAAGCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGt              <  1:1559771/98‑1 (MQ=255)
                                                                                    aGCTAAGAGGACAGTCCTCTTAGCCCCCTCCAATCCCCAATCATTCATTAAACAAATCCATTGACATAAAATATAAAAAAAAGCCCTTGCTTTCTAACGt              <  1:2786550/100‑1 (MQ=255)
                                                                                     gctacaagGACAGTCCTCATAGACCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAAc                <  1:2174876/91‑1 (MQ=255)
                                                                                     gCTAAGAGGACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGt              >  1:2388327/1‑99 (MQ=255)
                                                                                     gCTAAGAGGACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTg             >  1:1598954/1‑100 (MQ=255)
                                                                                     gCTAAGAGGACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTg             >  2:2318312/1‑100 (MQ=255)
                                                                                     gCTAAGAGGACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTg             >  2:1755548/1‑100 (MQ=255)
                                                                                      cTAAGAGGACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGa            >  2:1487191/1‑100 (MQ=255)
                                                                                       tAAGAGGACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGaa           >  2:1715337/1‑100 (MQ=255)
                                                                                       tAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGaa           >  1:1459847/1‑100 (MQ=255)
                                                                                       tAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGaa           >  1:381342/1‑100 (MQ=255)
                                                                                            ggACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTgg      >  2:2132844/1‑100 (MQ=255)
                                                                                            ggACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTgg      >  2:1495741/1‑100 (MQ=255)
                                                                                            ggACAATCCGCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTgg      <  1:2230134/100‑1 (MQ=255)
                                                                                               cAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGttt   <  1:2750240/100‑1 (MQ=255)
                                                                                                aGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTa  >  1:1589644/1‑100 (MQ=255)
                                                                                                |                                                                                                   
TAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTA  >  NC_000913/2089758‑2089953

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