Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,089,198 T→A 5.3% intergenic (‑70/+13) yeeZ ← / ← yoeB putative epimerase/toxin of the YoeB‑YefM toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,089,1980TA 5.3% 321.7 / 7.5 94intergenic (‑70/+13)yeeZ/yoeBputative epimerase/toxin of the YoeB‑YefM toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (46/43);  new base A (2/3);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

Only 100 of 104 total aligned reads displayed.

TTGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTAC  >  NC_000913/2089051‑2089344
                                                                                                                                                   |                                                                                                                                                   
ttGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTatg                                                                                                                                                   >  4:175981/1‑149 (MQ=255)
ttGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTatg                                                                                                                                                   >  1:508098/1‑149 (MQ=255)
  gCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTatgaa                                                                                                                                                 <  1:595143/149‑1 (MQ=255)
  gCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTatgaa                                                                                                                                                 <  2:626206/149‑1 (MQ=255)
      gCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAtgaatgaa                                                                                                                                             <  2:132789/149‑1 (MQ=255)
      gCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCATagaatgaa                                                                                                                                             <  1:164399/149‑1 (MQ=255)
       ccTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATgaatgaat                                                                                                                                            <  2:144263/149‑1 (MQ=255)
       ccTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATgaatgaat                                                                                                                                            <  4:404078/149‑1 (MQ=255)
        cTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCATAgaatgaatt                                                                                                                                           <  3:463025/149‑2 (MQ=255)
         tCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCt                                                                                                                                          >  4:434748/1‑149 (MQ=255)
         tCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCt                                                                                                                                          <  1:228401/149‑1 (MQ=255)
            cgcTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGtt                                                                                                                                       <  2:543833/149‑1 (MQ=255)
            cgcTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCATAgaatgaattgtca                                                                                                                                       <  3:113196/149‑6 (MQ=255)
             gcTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTc                                                                                                                                      <  2:508098/149‑1 (MQ=255)
                   aGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCaataat                                                                                                                                <  2:609069/149‑1 (MQ=255)
                      gACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGTGACCTTTTTCATCCGTTACCTCCTGATACACCATAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATGTGTTCAATAATGAt                                                                                                                             >  3:633712/1‑149 (MQ=255)
                      gACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGAt                                                                                                                             >  2:327999/1‑149 (MQ=255)
                       aCATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATa                                                                                                                            >  1:426391/1‑149 (MQ=255)
                         aTCGCCAGCGGCATGCCTACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACg                                                                                                                         >  2:63507/1‑149 (MQ=255)
                          tCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACg                                                                                                                         >  1:454029/1‑149 (MQ=255)
                           cGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGa                                                                                                                        <  3:295447/149‑1 (MQ=255)
                             ccAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTatgaat                                                                                                                                                <  2:38441/123‑1 (MQ=255)
                             ccAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTatgaat                                                                                                                                                >  1:38441/1‑123 (MQ=255)
                               aGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCACCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACAtg                                                                                                                    <  2:454029/149‑1 (MQ=255)
                                   gCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCATAgaatgaatt                                                                                                                                           <  1:747376/122‑2 (MQ=255)
                                   gCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCATAgaatgaatt                                                                                                                                           >  2:747376/1‑121 (MQ=255)
                                     aTGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCaa                                                                                                              >  1:364158/1‑149 (MQ=255)
                                           aaCCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGa                                                                                                                              <  3:659091/127‑1 (MQ=255)
                                           aaCCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGa                                                                                                                              >  4:659091/1‑127 (MQ=255)
                                           aaCCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCa                                                                                                        >  3:473983/1‑149 (MQ=255)
                                            accaccCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAg                                                                                                       >  3:753629/1‑149 (MQ=255)
                                            accaccCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAg                                                                                                       <  4:290052/149‑1 (MQ=255)
                                               accCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCaataat                                                                                                                                <  4:409981/121‑1 (MQ=255)
                                                 ccTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAAt                                                                                                  <  1:680961/149‑1 (MQ=255)
                                                     aCCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAAtcgtc                                                                                              <  2:85081/149‑1 (MQ=255)
                                                     aCCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAAtcgtc                                                                                              <  4:131195/149‑1 (MQ=255)
                                                      cccGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAAgcgtct                                                                                             <  4:516808/149‑2 (MQ=255)
                                                        cGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGt                                                                                           >  1:319711/1‑149 (MQ=255)
                                                         gACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTa                                                                                          >  2:35852/1‑149 (MQ=255)
                                                                 cGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGta                                                                                  <  4:472235/149‑1 (MQ=255)
                                                                  gACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGtat                                                                                 >  2:418167/1‑149 (MQ=255)
                                                                   aCCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTAAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGtata                                                                                >  3:654818/1‑149 (MQ=255)
                                                                     cTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATAcc                                                                              >  4:676598/1‑149 (MQ=255)
                                                                      tttttCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCa                                                                             >  2:434801/1‑149 (MQ=255)
                                                                          tCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACg                                                                         >  1:508814/1‑149 (MQ=255)
                                                                           cATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACgg                                                                        >  3:605124/1‑149 (MQ=255)
                                                                             tCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTg                                                                      >  2:202468/1‑149 (MQ=255)
                                                                               cGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGct                                                                    >  3:700743/1‑149 (MQ=255)
                                                                                  tATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGctcct                                                                 <  1:63507/149‑1 (MQ=255)
                                                                                      tcCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGt                                                             >  1:333797/1‑149 (MQ=255)
                                                                                      tcCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGt                                                             >  3:218814/1‑149 (MQ=255)
                                                                                       ccTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTa                                                            <  1:149038/149‑1 (MQ=255)
                                                                                        cTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTaa                                                           <  3:536281/149‑1 (MQ=255)
                                                                                        cTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTaa                                                           >  1:738730/1‑149 (MQ=255)
                                                                                              cacCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTc                                                     >  1:459087/1‑149 (MQ=255)
                                                                                               acCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCg                                                    <  4:146521/149‑1 (MQ=255)
                                                                                                ccTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCgg                                                   >  1:189838/1‑149 (MQ=255)
                                                                                                 cTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCggg                                                  >  2:103577/1‑149 (MQ=255)
                                                                                                  ttAGATCTATAAGGCTACGCTAGCGTATCATAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGa                                                 >  3:232430/1‑149 (MQ=255)
                                                                                                  ttAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCgg                                                   >  2:645726/1‑147 (MQ=255)
                                                                                                  ttAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCgg                                                   <  1:645726/147‑1 (MQ=255)
                                                                                                   tAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGAc                                                <  1:410250/149‑1 (MQ=255)
                                                                                                      aTCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAg                                             >  1:207456/1‑149 (MQ=255)
                                                                                                         tataAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGaagccaga                                          >  4:585303/1‑142 (MQ=255)
                                                                                                          ataAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAAc                                         >  2:764454/1‑149 (MQ=255)
                                                                                                           taAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAAcc                                        >  1:771634/1‑149 (MQ=255)
                                                                                                            aaGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCATTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGACACCt                                       <  2:587158/149‑1 (MQ=255)
                                                                                                            aaGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCt                                       <  1:133450/149‑1 (MQ=255)
                                                                                                                 tACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACaa                                  <  4:103151/149‑1 (MQ=255)
                                                                                                                  aCGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACaaa                                 <  1:327999/149‑1 (MQ=255)
                                                                                                                    gCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAAtt                               >  2:531471/1‑149 (MQ=255)
                                                                                                                        gCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGt                           >  1:496254/1‑149 (MQ=255)
                                                                                                                          gTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCt                                       >  2:520164/1‑135 (MQ=255)
                                                                                                                          gTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCt                                       <  1:520164/135‑1 (MQ=255)
                                                                                                                            aTCAAAACTGACAATTCATTCTATGAATTAATCTGTTAAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGATTACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATGATGTTTCa                       >  4:149847/1‑149 (MQ=255)
                                                                                                                            aTCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCa                       <  1:495932/149‑1 (MQ=255)
                                                                                                                             tCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTc                                                     <  1:479853/118‑1 (MQ=255)
                                                                                                                             tCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTc                                                     >  2:479853/1‑118 (MQ=255)
                                                                                                                             tCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAg                      <  4:605124/149‑1 (MQ=255)
                                                                                                                             tCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAg                      <  2:279458/149‑1 (MQ=255)
                                                                                                                              cAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAgg                     <  2:319711/149‑1 (MQ=255)
                                                                                                                              cAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAgg                     <  2:480078/149‑1 (MQ=255)
                                                                                                                              cAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAgg                     <  1:661096/149‑1 (MQ=255)
                                                                                                                               aaaaCTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAggg                    >  3:589956/1‑149 (MQ=255)
                                                                                                                                aaaCTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTAt                             <  1:702685/139‑1 (MQ=255)
                                                                                                                                aaaCTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTAt                             >  2:702685/1‑139 (MQ=255)
                                                                                                                                aaaCTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAgggg                   >  2:424607/1‑149 (MQ=255)
                                                                                                                                 aaCTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGt                  >  1:548706/1‑149 (MQ=255)
                                                                                                                                  aCTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGtt                 <  2:189838/149‑1 (MQ=255)
                                                                                                                                   cTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTc                >  2:127642/1‑149 (MQ=255)
                                                                                                                                    tGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCt               <  2:207456/149‑1 (MQ=255)
                                                                                                                                      aCAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTgg             <  2:281607/149‑1 (MQ=255)
                                                                                                                                       cAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGc            <  3:140814/149‑1 (MQ=255)
                                                                                                                                       agattcattcATAGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACcgc                                                                                     <  1:519513/74‑1 (MQ=38)
                                                                                                                                         attcattcATAGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGc            >  2:590741/1‑147 (MQ=255)
                                                                                                                                                 tATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTAc  >  3:164631/1‑149 (MQ=255)
                                                                                                                                                 tATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTAc  <  2:738730/149‑1 (MQ=255)
                                                                                                                                                 tATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTAc  >  1:178879/1‑149 (MQ=255)
                                                                                                                                                  atgaatgaATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGc                                                                     >  2:409000/1‑81 (MQ=255)
                                                                                                                                                  atgaatgaATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGc                                                                     <  1:409000/81‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTAC  >  NC_000913/2089051‑2089344

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