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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,089,8540AT 6.5% 423.2 / 12.9 136intergenic (‑141/‑142)yefM/hisLYefM antitoxin of the YoeB‑YefM toxin‑antitoxin pair and DNA binding transcriptional repressor/his operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (69/58);  new base T (5/4);  total (74/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 136 total aligned reads displayed.

TACGCATGTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACAC  >  NC_000913/2089707‑2090002
                                                                                                                                                   |                                                                                                                                                    
tACGCATGTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGAGTTAATTAACAATTTATTAATTAGTTTGTAGATCAAGCTATTGACAGTGAGACGAAAATAAAGGCTTCGCTATTTATCGTCACATCAGCTAAGAGGACAGt                                                                                                                                                    >  1:903177/1‑150 (MQ=255)
 aCGCATGTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCCGCTAAGAGGACAGTc                                                                                                                                                   >  2:1278416/1‑150 (MQ=255)
    cATGTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCctc                                                                                                                                                >  2:176059/1‑150 (MQ=255)
     aTGTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGCGCGTTAATTAACTATTTATTAATTAGTTTGTAGCTCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGCCAGTCctct                                                                                                                                               >  2:1322616/1‑150 (MQ=255)
     aTGTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCctct                                                                                                                                               <  2:825155/150‑1 (MQ=255)
       gTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTa                                                                                                                                             <  2:1001553/150‑1 (MQ=255)
       gTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTCTCAGTAACAATAAAAAAAATAATTCGCTATTTTTGGTGCACTCAGCTAAGAGGACAGTCCTCTTa                                                                                                                                             >  2:621146/1‑150 (MQ=255)
        tCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAg                                                                                                                                            >  1:1022216/1‑150 (MQ=255)
         cAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGc                                                                                                                                           <  2:251121/150‑1 (MQ=255)
         cAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGc                                                                                                                                           >  2:1444581/1‑150 (MQ=255)
         cAATCTCCTCTTTTGTACAGATCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGc                                                                                                                                           <  1:1066308/150‑1 (MQ=255)
          aaTCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACTGTCCTCTTAGct                                                                                                                                          <  1:993052/150‑2 (MQ=255)
          aaTCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTcc                                                                                                                                                  <  1:1143881/142‑1 (MQ=255)
              ctcatcattGTACAATGAACTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCcct                                                                                                                                      >  2:1188298/8‑150 (MQ=255)
               cctcTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCcctc                                                                                                                                     <  1:1537159/150‑1 (MQ=255)
                   ttttGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGccccc                                                                                                                                       >  1:266728/1‑144 (MQ=255)
                    tttGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCttta                                                                                                                                >  1:770901/1‑149 (MQ=255)
                       gTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTcccc                                                                                                                             >  2:112301/1‑150 (MQ=255)
                       gTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTcccc                                                                                                                             >  2:728599/1‑150 (MQ=255)
                        tACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCg                                                                                                                            >  1:1447697/1‑150 (MQ=255)
                           aGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtc                                                                                                                         <  2:186847/150‑1 (MQ=255)
                           aGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCtc                                                                                                                         >  2:1080110/1‑150 (MQ=255)
                                      cAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAAc                                                                                                              <  2:310454/150‑1 (MQ=255)
                                         tgatgaGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACaaa                                                                                                           >  1:27756/1‑150 (MQ=255)
                                             gaGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCa                                                                                                       >  2:643199/1‑150 (MQ=255)
                                                cGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTg                                                                                                    <  1:1230154/150‑1 (MQ=255)
                                                  ttaattaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGcc                                                                                                  <  1:1045730/150‑1 (MQ=255)
                                                      ttaaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATaa                                                                                              >  1:467553/1‑150 (MQ=255)
                                                       taaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATaaa                                                                                             >  1:206551/1‑150 (MQ=255)
                                                        aaCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCAAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATaaaa                                                                                            >  1:851004/1‑150 (MQ=255)
                                                         aCTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAat                                                                                           <  1:86878/150‑1 (MQ=255)
                                                          cTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAata                                                                                          <  1:1234318/150‑1 (MQ=255)
                                                          cTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAata                                                                                          <  2:1031423/150‑1 (MQ=255)
                                                            atttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAatata                                                                                        <  2:734550/150‑1 (MQ=255)
                                                            atttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAatata                                                                                        >  2:516546/1‑150 (MQ=255)
                                                             tttattAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCAAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAatatat                                                                                       <  2:851004/150‑1 (MQ=255)
                                                               tattaATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATaa                                                                                     >  1:56770/1‑150 (MQ=255)
                                                                attaattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATaaa                                                                                    <  1:1444581/150‑1 (MQ=255)
                                                                 ttaattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATaaaa                                                                                   >  2:752418/1‑150 (MQ=255)
                                                                  taattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATaaaaa                                                                                  <  1:736012/150‑1 (MQ=255)
                                                                   aattaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAAt                                                                                                          >  2:1214831/1‑124 (MQ=255)
                                                                    attaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAg                                                                                <  1:112301/150‑1 (MQ=255)
                                                                     ttaGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGc                                                                               >  1:520610/1‑150 (MQ=255)
                                                                     ttaGTTTGTAGATCAAGGTATTGTCAGTGACACGAAAATCCAGGCTGCGCTATTTTTGGTGCCATCAGCTACGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGc                                                                               <  2:887053/150‑1 (MQ=255)
                                                                        gTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCt                                                                            <  1:994392/150‑1 (MQ=255)
                                                                         tttGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCtt                                                                           >  1:404562/1‑150 (MQ=255)
                                                                         tttGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGACATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTAATAAAACAAATCCATTGCCATAAAATATATTAAAAAGCCCtt                                                                           >  2:702830/1‑150 (MQ=255)
                                                                          ttGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTg                                                                          <  2:1022216/150‑1 (MQ=255)
                                                                           tGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGc                                                                         <  2:404562/150‑1 (MQ=255)
                                                                              aGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCttt                                                                      >  2:1105300/1‑150 (MQ=255)
                                                                                aTCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCt                                                                    >  2:955552/1‑150 (MQ=255)
                                                                                     ggTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGt                                                               <  2:1068450/150‑1 (MQ=255)
                                                                                      gTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTg                                                              <  2:27756/150‑1 (MQ=255)
                                                                                          tGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCtt                                                                       <  2:1112315/137‑1 (MQ=255)
                                                                                          tGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCtt                                                                       >  1:1112315/1‑137 (MQ=255)
                                                                                          tGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAg                                                          >  2:720916/1‑150 (MQ=255)
                                                                                          tGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAg                                                          <  2:39149/150‑1 (MQ=255)
                                                                                           gTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGt                                                         >  2:139514/1‑150 (MQ=255)
                                                                                             cAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCGATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTgg                                                       <  1:1223128/150‑1 (MQ=255)
                                                                                              aGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGt                                                      >  1:1128528/1‑150 (MQ=255)
                                                                                               gTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGtt                                                     >  1:1142415/1‑150 (MQ=255)
                                                                                                tGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCGCTTGCTTTCTAACGTGAAAGTGGgtt                                                    >  2:1076512/1‑150 (MQ=255)
                                                                                                tGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGttt                                                    <  1:955552/150‑1 (MQ=255)
                                                                                                 gagaCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTa                                                   >  1:528032/1‑150 (MQ=255)
                                                                                                   gccGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGCTTTAgg                                                 >  1:280132/3‑150 (MQ=255)
                                                                                                    actaaaaTCCAGTCTTCTCTATTTTTTTTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGt                                                <  1:1272962/147‑1 (MQ=255)
                                                                                                    aCGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGt                                                <  1:720916/150‑1 (MQ=255)
                                                                                                     cGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGtt                                               >  2:618907/1‑150 (MQ=255)
                                                                                                     cGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGtt                                               <  1:643199/150‑1 (MQ=255)
                                                                                                         aaTCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTaaaa                                           <  1:176059/150‑1 (MQ=255)
                                                                                                            ccAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGAc                                        >  2:882333/1‑150 (MQ=255)
                                                                                                            ccAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGAc                                        >  2:888159/1‑150 (MQ=255)
                                                                                                             cAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAATACa                                       <  2:233807/150‑1 (MQ=255)
                                                                                                                gCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCa                                    >  2:191740/1‑150 (MQ=255)
                                                                                                                    cGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTg                                <  2:1612329/150‑1 (MQ=255)
                                                                                                                        aTTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATa                            <  1:313943/150‑1 (MQ=255)
                                                                                                                        aTTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATa                            >  1:1419899/1‑150 (MQ=255)
                                                                                                                        aTTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCACCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATa                            <  1:1297725/150‑1 (MQ=255)
                                                                                                                          ttttGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATaaa                          >  2:1351166/1‑150 (MQ=255)
                                                                                                                           tttGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAAc                         <  2:271936/150‑1 (MQ=255)
                                                                                                                           tttGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAAc                         >  2:664254/1‑150 (MQ=255)
                                                                                                                            ttGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACa                        <  1:882333/150‑1 (MQ=255)
                                                                                                                              ggTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACAtt                      >  1:452157/1‑150 (MQ=255)
                                                                                                                              ggTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACAtt                      >  1:1629292/1‑150 (MQ=255)
                                                                                                                              ggTCCCATCAGCTAAGAGGACACTCCTCTTAGACCCCTCCTTTACCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTATAACGAGAAAGTGGTTTAGGTTAAAAGACATAAGTTGAATAAACAtt                      >  2:1596577/1‑150 (MQ=255)
                                                                                                                               gTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTc                     <  1:1188298/150‑1 (MQ=255)
                                                                                                                                 gCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCt                                                                    <  1:1564068/101‑1 (MQ=255)
                                                                                                                                 gCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCt                                                                    >  2:1564068/1‑101 (MQ=255)
                                                                                                                                 gCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTcac                   <  2:1629292/150‑1 (MQ=255)
                                                                                                                                   cATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACag                 <  2:496564/150‑1 (MQ=255)
                                                                                                                                       ggCTAAGAGGACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTc                     <  1:1572270/141‑1 (MQ=255)
                                                                                                                                       ggCTAAGAGGACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTc                     >  2:1572270/2‑142 (MQ=255)
                                                                                                                                          tAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACttt          >  1:954769/1‑150 (MQ=255)
                                                                                                                                          tAAGAGGACAGTCCTCTTAGCCACCTCCTTTCCCCGATAATTCATTAAAAAAAAACATTGCAATAAAAAATATAAAAAAGAAAATGAATTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACttt          >  2:982652/1‑150 (MQ=255)
                                                                                                                                            agagGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTa        <  2:1447697/150‑1 (MQ=255)
                                                                                                                                              agGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATg      <  2:1546664/150‑1 (MQ=255)
                                                                                                                                               ggACTGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTTACAGAGACTTTTATGa     >  1:757259/1‑150 (MQ=255)
                                                                                                                                               ggACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGa     >  1:488646/1‑150 (MQ=255)
                                                                                                                                               ggACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATACACAGAGACTTTTATGa     >  2:578887/1‑150 (MQ=255)
                                                                                                                                                  cAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGacac  >  2:570650/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TACGCATGTCAATCTCCTCTTTTGTACAGTTCATTGTACAATGATGAGCGTTAATTAACTATTTATTAATTAGTTTGTAGATCAAGGTATTGTCAGTGAGACGAAAATCCAGGCTTCGCTATTTTTGGTGCCATCAGCTAAGAGGACAGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACAC  >  NC_000913/2089707‑2090002

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