Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,089,854 A→T 21.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,8540AT21.2% 283.8 / 74.8 136intergenic (‑141/‑142)yefM/hisLantitoxin of the YoeB‑YefM toxin‑antitoxin system/his operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (60/47);  new base T (14/15);  total (74/62)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-02

Only 100 of 137 total aligned reads displayed.

TTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTA  >  NC_000913/2089757‑2089953
                                                                                                 |                                                                                                   
ttaattaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGt                                                                                                   >  1:469778/1‑100 (MQ=255)
ttaattaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGt                                                                                                   <  1:2469585/100‑1 (MQ=255)
  aattaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTcc                                                                                                 >  1:93085/1‑100 (MQ=255)
  aattaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTcc                                                                                                 <  1:78781/100‑1 (MQ=255)
      aaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCtt                                                                                             <  1:1006969/100‑1 (MQ=255)
       aCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTa                                                                                            <  1:2623189/100‑1 (MQ=255)
       aCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTa                                                                                            <  2:2603526/100‑1 (MQ=255)
        cTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAgatcggaaga                                                                                           >  1:1493971/1‑91 (MQ=255)
         tatttattAATTAGTTTGTAGATCAAGTTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGc                                                                                          <  1:673584/100‑1 (MQ=255)
         tatttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGATACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGc                                                                                          >  2:3248761/1‑100 (MQ=255)
         tatttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGc                                                                                          <  1:523422/100‑1 (MQ=255)
         tatttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGc                                                                                          <  1:3047850/100‑1 (MQ=255)
         tatttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGc                                                                                          <  1:208394/100‑1 (MQ=255)
          atttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGct                                                                                         <  2:803299/100‑2 (MQ=255)
          atttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGct                                                                                         <  1:2763985/100‑2 (MQ=255)
          atttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGcc                                                                                         <  1:1845127/100‑1 (MQ=255)
           ttttttAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATTCAGTCTTCGCTATTTTTTTTGTCTTTATCTTAGAGGTCTGTCCTTTTAGctg                                                                                        <  2:1661801/100‑3 (MQ=255)
           tttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGccc                                                                                        <  1:2633927/100‑1 (MQ=255)
           tttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGccc                                                                                        >  2:1465204/1‑100 (MQ=255)
            ttagtaATTAGTTTGTAGATAAAGGTATTGCTAGTGAGACGAAAATCCAGACAAGGCTATTTTTGGAGCCATCAGAGAAGAGGACAGTcaaatt                                                                                             >  2:3527608/1‑89 (MQ=255)
             tattaATTAGTTTGTAGATCAAGGTATTGTCAGTGAGATGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGccccc                                                                                      <  1:3385330/100‑1 (MQ=255)
             tattaATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGccccc                                                                                      >  1:1341966/1‑100 (MQ=255)
               ttaattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCctcctc                                                                                    >  2:300912/1‑100 (MQ=255)
               ttaattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCcctc                                                                                    <  2:68358/100‑1 (MQ=255)
               ttaattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCcct                                                                                     <  1:2471970/99‑1 (MQ=255)
                taattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGctgatggc                                                                                   <  1:731577/100‑8 (MQ=38)
                taattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCcctcc                                                                                   <  1:886951/100‑1 (MQ=255)
                taattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCcctcc                                                                                   >  1:3150737/1‑100 (MQ=255)
                 aattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCcctcct                                                                                  >  2:728574/1‑100 (MQ=255)
                 aattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCcctcct                                                                                  <  2:376651/100‑1 (MQ=255)
                  attaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCtt                                                                                 <  1:1775877/100‑1 (MQ=255)
                    taGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTc                                                                               <  1:1817916/100‑1 (MQ=255)
                      gTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTccc                                                                             >  2:1320045/1‑100 (MQ=255)
                      gTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTccc                                                                             >  1:529173/1‑100 (MQ=255)
                         tGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGc                                                                          <  1:443386/100‑1 (MQ=255)
                              aTCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGTTTTCATTATCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcatt                                                                     >  2:1224808/1‑100 (MQ=255)
                              aTCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcatt                                                                     <  2:2685136/100‑1 (MQ=255)
                              aTCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcatt                                                                     >  1:2770554/1‑100 (MQ=255)
                               tCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtcattc                                                                    <  1:2770303/100‑1 (MQ=255)
                                  aGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTcattcatt                                                                 >  2:2056423/1‑100 (MQ=255)
                                   ggTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTa                                                                <  1:2079734/100‑1 (MQ=255)
                                      aTTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAAc                                                             <  1:2238158/100‑1 (MQ=255)
                                         gTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACaaa                                                          >  1:1632930/1‑100 (MQ=255)
                                           cAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATc                                                        <  1:1166626/100‑1 (MQ=255)
                                             gTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCa                                                      >  1:3224334/1‑100 (MQ=255)
                                               gagaCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCACCCTCCTTTCCCCGCTCATTCATTAAACAAATCAAtt                                                    >  1:225785/1‑100 (MQ=255)
                                                agaCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTg                                                   >  2:2794159/1‑100 (MQ=255)
                                                  aCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGcc                                                 <  2:1632915/100‑1 (MQ=255)
                                                   cGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGc                                                  >  1:1672029/1‑98 (MQ=255)
                                                   cGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCa                                                <  1:1851737/100‑1 (MQ=255)
                                                   cGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCa                                                >  1:3356041/1‑100 (MQ=255)
                                                   cGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCa                                                >  2:1808425/1‑100 (MQ=255)
                                                     aaaaTCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATa                                              >  1:3335532/1‑100 (MQ=255)
                                                     aaaaTCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATa                                              >  2:1747326/1‑100 (MQ=255)
                                                        aTCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATaaaa                                           >  2:2450763/1‑100 (MQ=255)
                                                        aTCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATaaaa                                           >  1:2356186/1‑100 (MQ=255)
                                                        aTCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATaaaa                                           >  1:3829678/1‑100 (MQ=255)
                                                         tCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAat                                          <  2:3804411/100‑1 (MQ=255)
                                                           cAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAatat                                        >  2:3443676/1‑100 (MQ=255)
                                                           cAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAatat                                        >  1:3896254/1‑100 (MQ=255)
                                                                 tCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTATTAGCCAACTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATaaaa                                  >  2:3692182/1‑100 (MQ=255)
                                                                    cTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAg                               <  2:2160426/100‑1 (MQ=255)
                                                                    cTATTTTTGGTGCCATCAGCAAAGATGACAGTCCTCTTAGCTACCTGCTTTCCCCGCTAATTCATTAAACAAAACCATTACCATAAAATATATAAAAAAg                               >  2:1727076/1‑100 (MQ=255)
                                                                      aTTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCTCTTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGcc                             <  1:726987/100‑1 (MQ=255)
                                                                      aTTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGcc                             <  1:3244932/100‑1 (MQ=255)
                                                                      aTTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGcc                             <  1:3777342/100‑1 (MQ=255)
                                                                       tttttGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGcc                             >  1:1084001/1‑99 (MQ=255)
                                                                             gTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCtt                      >  2:685453/1‑100 (MQ=255)
                                                                             gTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCtt                      <  1:3440075/100‑1 (MQ=255)
                                                                              tGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCttt                     <  2:3451054/100‑1 (MQ=255)
                                                                               gCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTc                    <  1:1673847/100‑1 (MQ=255)
                                                                                ccATCAGCTAAGAGGACAGTACTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTACCATCACATATATAAAAAAGACCTTGCTTTCt                   >  2:1565290/1‑100 (MQ=255)
                                                                                   tCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAAc                >  1:3582137/1‑100 (MQ=255)
                                                                                    gggCTAAGAGGACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACg               >  2:3483468/3‑100 (MQ=255)
                                                                                    gggCTAAGAGGACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACg               >  2:3791678/3‑100 (MQ=255)
                                                                                    cAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACg               <  2:2030281/100‑1 (MQ=255)
                                                                                    cAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACg               <  1:566912/100‑1 (MQ=255)
                                                                                     ggCTAAGAGGACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGt              >  2:3890595/2‑100 (MQ=255)
                                                                                     ggCTAAGAGGACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGt              >  2:3612689/2‑100 (MQ=255)
                                                                                     ggCTAAGAGGACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAAACCATTACAATAAAATAAAAAAAAAAGCCATTGCTTTCTAACg               >  1:204125/2‑99 (MQ=255)
                                                                                     agaaaaGAGAACAGGCCTAATAACGCCCTCCTTCCCCCACTCAATAATTTAACCGATCTTTTGTGAAAAAAAAAAAAAAAAAGCCCTTGCTTTCTAACGt              <  2:1900357/96‑1 (MQ=255)
                                                                                      gatAAGAGGACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGt              >  1:2658934/3‑99 (MQ=255)
                                                                                      gCTAAGAGGACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAAAAAATCCAATACCAAAAAATATATAAAAAAGCACTTGATTTCTAACGt              >  1:5984/1‑99 (MQ=255)
                                                                                        tAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGaa           >  2:2199853/1‑100 (MQ=255)
                                                                                        tAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGaa           >  1:2990199/1‑100 (MQ=255)
                                                                                        tAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGaa           <  1:3569881/100‑1 (MQ=255)
                                                                                         aaGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGaaa          <  2:1812418/100‑1 (MQ=255)
                                                                                          agagGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAg         <  2:2827723/100‑1 (MQ=255)
                                                                                           gagGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGt        >  2:1038991/1‑100 (MQ=255)
                                                                                            agGATAGTCCTCTTAGCTCCCTCATTTCCCCTCTAATTCATTAAACAAATCCATTGCAATAAAATATATAAAAAATTTCTTGCTTTCTAACGTGAAAg         <  2:571589/98‑1 (MQ=255)
                                                                                            agGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTg       >  2:3460755/1‑100 (MQ=255)
                                                                                             ggACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTgg      >  1:1954240/1‑100 (MQ=255)
                                                                                             ggACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTgg      >  2:1301416/1‑100 (MQ=255)
                                                                                                cAGTCCTCTTAGCTACCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGTCCTTGCTTTCTAACTTGAAAGTGGtt    >  2:1639745/1‑99 (MQ=255)
                                                                                                cAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGttt   <  1:1381314/100‑1 (MQ=255)
                                                                                                cAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGttt   >  1:972200/1‑100 (MQ=255)
                                                                                                cAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGttt   <  1:3364339/100‑1 (MQ=255)
                                                                                                cAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGttt   >  2:2229295/1‑100 (MQ=255)
                                                                                                cAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCATGCTTTCTAACGTGAAAGTGGtt    >  2:2419471/1‑99 (MQ=255)
                                                                                                 aGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTa  <  1:3492398/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
TTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTA  >  NC_000913/2089757‑2089953

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