Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,011 T→A 17.8% intergenic (+171/+45) ybfK → / ← kdpE uncharacterized protein/response regulator in two‑component regulatory system with KdpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913721,0110TA17.8% 364.2 / 68.5 152intergenic (+171/+45)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (58/66);  new base A (16/11);  total (74/78)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

Only 100 of 157 total aligned reads displayed.

ATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGG  >  NC_000913/720914‑721110
                                                                                                 |                                                                                                   
aTCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatta                                                                                                   <  2:1114260/100‑1 (MQ=255)
  cTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                                                                                                 <  1:2804534/100‑1 (MQ=255)
  cTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                                                                                                 <  2:866270/100‑1 (MQ=255)
    tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                               >  1:1322392/1‑100 (MQ=255)
    tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                               >  1:2082005/1‑100 (MQ=255)
       tttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                            >  1:1768455/1‑100 (MQ=255)
       tttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttt                                                                                            <  1:3159952/100‑1 (MQ=38)
        ttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                           <  2:685491/100‑2 (MQ=38)
           tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                        <  1:951674/100‑1 (MQ=255)
           tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                        <  2:2072970/100‑1 (MQ=255)
            ttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatta                                                                                       <  2:3393693/100‑1 (MQ=255)
             taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                                                                                      >  2:2606144/1‑100 (MQ=255)
             taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                      <  1:2045797/100‑7 (MQ=25)
              aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                                     <  2:3321099/100‑8 (MQ=25)
              aaTGATGAGCTGCTTACTTAGTACTATACAATATTGAATTGTTATTTTTGTGTGTTTTTTAAGATAAAAAGCCGTATTTATTTTTAAGGCTtta                                                                                           >  2:2353628/1‑94 (MQ=255)
               atgattAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTTTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAGTACGGCTTTAATtagta                                                                                    <  2:3910556/100‑1 (MQ=255)
               atgatgAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                    >  1:1155716/1‑100 (MQ=255)
                   tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCAGTATTTATTATTACGGCTTTAATTAATAAAag                                                                                >  2:56313/1‑100 (MQ=255)
                    gaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                               >  1:2474767/1‑100 (MQ=255)
                    gaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                               >  1:2505501/1‑100 (MQ=255)
                     aGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                               <  2:2952240/99‑1 (MQ=255)
                      gCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                             <  1:3636910/100‑1 (MQ=255)
                      gCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                             <  1:1217243/100‑1 (MQ=255)
                       cTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag                                                                            <  2:1645189/100‑1 (MQ=255)
                       cTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag                                                                            >  1:1553393/1‑100 (MQ=255)
                        tttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcagg                                                                           <  2:779298/100‑1 (MQ=255)
                         ttttACTTAGTAATATAAAATATTGCATCGCTATTTTTGTGTTCTTTTAAAGATAAAAAGCCGTCTTTATTATTACTGCTTTAATTAATAAAaggcagtc                                                                          >  2:2573159/1‑98 (MQ=255)
                         ttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc                                                                          <  1:3302847/100‑1 (MQ=255)
                            tACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTGTTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt                                                                       <  2:4380/100‑1 (MQ=255)
                            tACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt                                                                       <  1:1031481/100‑1 (MQ=255)
                             aCTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa                                                                      >  2:3994749/1‑100 (MQ=255)
                              cttaGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTATGGCTTTAATTAATAAAAGGAAGGCTGTAt                                                                     >  2:3739287/4‑100 (MQ=255)
                               ttaGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAtt                                                                    <  2:2291399/98‑1 (MQ=255)
                                taGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGTCTGTATTa                                                                   >  2:1728212/2‑100 (MQ=255)
                                taGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                                   >  1:576273/2‑100 (MQ=255)
                                taGTAATATAAAATATTGAAGTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAGTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAtt                                                                    <  2:1664735/98‑1 (MQ=255)
                                  gTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaa                                                                 <  2:2081992/100‑1 (MQ=255)
                                  gTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaa                                                                 <  1:568499/100‑1 (MQ=255)
                                    aaTATAAAATTTTGAATTGTTATTTTTGTGTGTTGTTTAAGTTTATAGGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                               <  1:1309864/100‑1 (MQ=255)
                                    aaTATAAAATATTGAGTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCATCCTGTATTAAAAt                                                               >  2:1576925/1‑100 (MQ=255)
                                    aaTATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                               <  1:2419440/100‑1 (MQ=255)
                                    aaTATAAAATATTGAATTGTTATTTTTGTGTGTTCTTTAAGATAAAAAGCCGTATTTATAATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                               <  2:928879/100‑1 (MQ=255)
                                     atatAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAtt                                                              <  1:2917685/100‑1 (MQ=255)
                                     atatAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAtt                                                              >  1:301744/1‑100 (MQ=255)
                                      tataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAATATTa                                                             <  2:1385387/100‑1 (MQ=255)
                                       ataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTaa                                                            <  2:3094434/100‑1 (MQ=255)
                                        taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTATATTAAAATTAat                                                           >  2:1504539/1‑100 (MQ=255)
                                         aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                          <  1:2744427/100‑1 (MQ=255)
                                         aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                          <  1:980437/100‑1 (MQ=255)
                                          aaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAatat                                                         >  1:725721/1‑100 (MQ=255)
                                          aaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAatat                                                         <  1:183699/100‑1 (MQ=255)
                                               ttGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaa                                                     >  1:1377283/1‑99 (MQ=255)
                                                 gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                                  <  1:1457840/100‑1 (MQ=255)
                                                 gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                                  <  1:1741658/100‑1 (MQ=255)
                                                 gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAACGCAGGCTGTATTAAAATTAATATTCAAAGc                                                  >  2:1369333/1‑100 (MQ=255)
                                                 gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCATATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAACATTAATATTCAAAGc                                                  >  2:3562876/1‑100 (MQ=255)
                                                  aaTTTTTATTTTTGTGTGTTGTTTAAGATATAAAGCTGTATTTATTATTAGGGCTTTAATTAATAAAGGGCAGGCGGTTTTAATATTAATATTCAAAg                                                   <  2:848374/98‑1 (MQ=255)
                                                   aTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCAt                                                >  2:2046988/1‑100 (MQ=255)
                                                      gttgttTTTGTGTGTTGTGTAAAATAAAAAGCCGTATTTATTATTACGTCTTTAATTAAAAAAAGGCAGGCTGTATTAAAATTAATATTCaaa                                                    <  2:3699688/93‑1 (MQ=255)
                                                      gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                             >  2:1990484/1‑100 (MQ=255)
                                                      gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                             >  2:177697/1‑100 (MQ=255)
                                                        tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGAAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                            >  2:463741/1‑99 (MQ=255)
                                                             ttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGATTTAATTAATAAAAGGAAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAg                                      <  2:1439336/100‑1 (MQ=255)
                                                               gtgtgtTGTTTAAGATAAAAAGTCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                                    <  2:1883244/100‑1 (MQ=255)
                                                                 gtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCaa                                  >  1:1929098/1‑100 (MQ=255)
                                                                  tgttgtTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAAt                                 <  2:1553381/100‑1 (MQ=255)
                                                                   gttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATa                                >  1:1536184/1‑100 (MQ=255)
                                                                      gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                             >  1:1362585/1‑100 (MQ=255)
                                                                      gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                             >  1:542197/1‑100 (MQ=255)
                                                                      gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                             <  2:3099271/100‑1 (MQ=255)
                                                                       tttattaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                            >  2:1249096/11‑100 (MQ=25)
                                                                       tttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                            <  1:3158601/100‑1 (MQ=255)
                                                                       tttAAGATAAAAAGCCGTATTTATTATTACGACTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgn                            <  2:305770/100‑2 (MQ=255)
                                                                        ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                           >  2:2341995/1‑100 (MQ=255)
                                                                        ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                           >  2:2116149/1‑100 (MQ=255)
                                                                         tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                          >  1:1096559/1‑100 (MQ=255)
                                                                         tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                          >  2:1763361/1‑100 (MQ=255)
                                                                          attaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                         >  2:2627164/8‑100 (MQ=25)
                                                                           aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                        <  2:2694598/100‑1 (MQ=255)
                                                                           aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                        <  1:2343850/100‑1 (MQ=255)
                                                                           aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAACGGTTTc                        >  1:2755882/1‑100 (MQ=255)
                                                                            gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                        <  2:618896/99‑1 (MQ=255)
                                                                             aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                      <  2:4124461/100‑1 (MQ=255)
                                                                              tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                     >  1:3485727/1‑100 (MQ=255)
                                                                              attaaaGCCGTAATAATAAATACGGCTTTACTTAATAAAAGGCAGGCTGTAATAAAATTAATATTCAAAGCATAAACCGATAGCCAAAACCTGTTTCact                     >  2:3418106/4‑98 (MQ=25)
                                                                              attaaaGCCGTAATAATAAATACGGCTTTAATTAATACAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                     >  2:2213128/4‑100 (MQ=37)
                                                                              attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                     >  1:2999198/4‑100 (MQ=37)
                                                                               ttaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTa                    >  2:385352/3‑100 (MQ=37)
                                                                               aaaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAACGGTTTCAGTa                    >  2:318112/1‑100 (MQ=255)
                                                                                taaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                   >  2:2277106/2‑100 (MQ=37)
                                                                                taaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                   >  2:1285925/2‑100 (MQ=37)
                                                                                taaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                   >  1:2593702/2‑100 (MQ=37)
                                                                                taaaGCCGTAATAATAAATAAGGCTTTAATTAATAAAAGGCAGTCTGCAATAAAATTAAGATTCAAAGCATACACAGATAGACAATACCGGtttacggaa                   >  2:3942283/2‑94 (MQ=255)
                                                                                 aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt                  >  1:2206954/1‑100 (MQ=37)
                                                                                       gTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCATTACTGaaacc            <  1:3604395/100‑3 (MQ=255)
                                                                                       gTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATg            <  2:4039762/100‑1 (MQ=255)
                                                                                       gTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATg            >  1:3884857/1‑100 (MQ=255)
                                                                                        tatttattATTACGGCTTTTATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATgg           >  2:2122228/1‑100 (MQ=255)
                                                                                         atttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGc          >  1:3002130/1‑100 (MQ=255)
                                                                                              ttattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATAGCGt        >  2:3520834/1‑97 (MQ=255)
                                                                                                 |                                                                                                   
ATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGG  >  NC_000913/720914‑721110

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