Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,089,196 T→A 6.1% intergenic (‑68/+15) 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,1960TA 6.1% 338.9 / 10.2 99intergenic (‑68/+15)yeeZ/yoeBputative epimerase/toxin of the YoeB‑YefM toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (47/45);  new base A (2/4);  total (49/50)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

Only 100 of 103 total aligned reads displayed.

GACTTGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTAC  >  NC_000913/2089048‑2089344
                                                                                                                                                    |                                                                                                                                                     
gACTTGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCt                                                                                                                                                       >  2:495932/1‑149 (MQ=255)
   ttGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTatg                                                                                                                                                   >  1:508098/1‑149 (MQ=255)
   ttGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTatg                                                                                                                                                   >  4:175981/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                                                                                                                                            <  1:589779/149‑1 (MQ=255)
          ccTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATgaatgaat                                                                                                                                            <  4:404078/149‑1 (MQ=255)
          ccTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATgaatgaat                                                                                                                                            <  2:144263/149‑1 (MQ=255)
           cTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCATAgaatgaatt                                                                                                                                           <  3:463025/149‑2 (MQ=255)
            tCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCt                                                                                                                                          <  1:228401/149‑1 (MQ=255)
            tCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCt                                                                                                                                          >  4:434748/1‑149 (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                                                                                                                                           >  2:747376/1‑121 (MQ=255)
                                      gCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCATAgaatgaatt                                                                                                                                           <  1:747376/122‑2 (MQ=255)
                                        aTGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCaa                                                                                                              >  1:364158/1‑149 (MQ=255)
                                              aaCCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGa                                                                                                                              >  4:659091/1‑127 (MQ=255)
                                              aaCCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGa                                                                                                                              <  3:659091/127‑1 (MQ=255)
                                              aaCCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCa                                                                                                        >  3:473983/1‑149 (MQ=255)
                                               accaccCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAg                                                                                                       <  4:290052/149‑1 (MQ=255)
                                               accaccCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAg                                                                                                       >  3:753629/1‑149 (MQ=255)
                                                  accCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCaataat                                                                                                                                <  4:409981/121‑1 (MQ=255)
                                                  accCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCaataat                                                                                                                                >  3:409981/1‑121 (MQ=255)
                                                    ccTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAAt                                                                                                  <  1:680961/149‑1 (MQ=255)
                                                        aCCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAAtcgtc                                                                                              <  4:131195/149‑1 (MQ=255)
                                                        aCCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAAtcgtc                                                                                              <  2:85081/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                                                             >  3:218814/1‑149 (MQ=255)
                                                                                         tcCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGt                                                             >  1:333797/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)
                                                                                                  acCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCg                                                    <  2:426391/149‑1 (MQ=255)
                                                                                                   ccTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCgg                                                   >  1:189838/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)
                                                                                                              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                                       <  1:520164/135‑1 (MQ=255)
                                                                                                                             gTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCt                                       >  2:520164/1‑135 (MQ=255)
                                                                                                                               aTCAAAACTGACAATTCATTCTATGAATTAATCTGTTAAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGATTACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATGATGTTTCa                       >  4:149847/1‑149 (MQ=255)
                                                                                                                               aTCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCa                       <  1:495932/149‑1 (MQ=255)
                                                                                                                                tCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTc                                                     >  2:479853/1‑118 (MQ=255)
                                                                                                                                tCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTc                                                     <  1:479853/118‑1 (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                             >  2:702685/1‑139 (MQ=255)
                                                                                                                                   aaaCTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTAt                             <  1:702685/139‑1 (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  <  2:738730/149‑1 (MQ=255)
                                                                                                                                                    tATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTAc  >  1:178879/1‑149 (MQ=255)
                                                                                                                                                    tATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTAc  >  3:164631/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACTTGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTAC  >  NC_000913/2089048‑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: