Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,009 T→A 18.5% intergenic (+169/+47) 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,0090TA18.5% 371.5 / 68.0 150intergenic (+169/+47)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (56/66);  new base A (16/12);  total (72/78)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 157 total aligned reads displayed.

ATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGC  >  NC_000913/720910‑721108
                                                                                                   |                                                                                                   
aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                                     <  1:2323442/100‑1 (MQ=255)
      cTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                                                                                               <  1:2804534/100‑1 (MQ=255)
      cTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                                                                                               <  2:866270/100‑1 (MQ=255)
        tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                             >  1:2082005/1‑100 (MQ=255)
        tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                             >  1:1322392/1‑100 (MQ=255)
          atttatttAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                           <  1:874371/100‑1 (MQ=255)
           tttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                          >  1:1768455/1‑100 (MQ=255)
            ttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                         <  2:685491/100‑2 (MQ=38)
               tttaATGATTAGTTTTTTACTTAGTAATTTAAAATATAGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                         <  2:3268449/97‑2 (MQ=25)
               tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                      <  1:951674/100‑1 (MQ=255)
               tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                      <  2:2072970/100‑1 (MQ=255)
                ttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatta                                                                                     <  2:3393693/100‑1 (MQ=255)
                 taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                                                                                    <  2:3681405/100‑1 (MQ=255)
                 taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                    <  1:2045797/100‑7 (MQ=25)
                 taATCATGTGCTTTTTACATAGTAATATAAAATATTTATTTGTTGTTTTTGTGTGTTGTGTCTGATACAGAGCCCTACTTTTTATTACGGCTttaattga                                                                                    <  2:3024111/100‑3 (MQ=255)
                  aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCTTATTTATTATTACGGCTTtaattaat                                                                                   <  1:2513465/100‑1 (MQ=255)
                  aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                                   <  1:3155074/100‑1 (MQ=255)
                  aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                                   <  2:3094922/100‑8 (MQ=25)
                   atgattAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTTTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAGTACGGCTTTAATtagta                                                                                  <  2:3910556/100‑1 (MQ=255)
                   atgatgAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                  <  1:3841850/100‑1 (MQ=255)
                   atgatgAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                  >  1:1155716/1‑100 (MQ=255)
                    tgatgaGCTTTTTCCTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaataa                                                                                 >  1:3258969/1‑100 (MQ=255)
                       tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCAGTATTTATTATTACGGCTTTAATTAATAAAag                                                                              >  2:56313/1‑100 (MQ=255)
                        gaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                             >  1:2474767/1‑100 (MQ=255)
                         aGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTTTTTTTTGTTTAATATAAAAAGCCGTATTTATTATTAATGCTTTAATTAATAAAag                                                                              >  2:2146071/1‑98 (MQ=255)
                          gCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                           <  1:1217243/100‑1 (MQ=255)
                           cTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag                                                                          >  1:1553393/1‑100 (MQ=255)
                           cTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag                                                                          <  2:1645189/100‑1 (MQ=255)
                            tttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcagg                                                                         <  2:779298/100‑1 (MQ=255)
                             ttttACTTAGTAATATAAAATATTGCATCGCTATTTTTGTGTTCTTTTAAAGATAAAAAGCCGTCTTTATTATTACTGCTTTAATTAATAAAaggcagtc                                                                        >  2:2573159/1‑98 (MQ=255)
                             ttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc                                                                        <  1:3302847/100‑1 (MQ=255)
                                tACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt                                                                     <  1:1031481/100‑1 (MQ=255)
                                 aCTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa                                                                    >  2:3994749/1‑100 (MQ=255)
                                   ttaGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAt                                                                   <  2:3965297/97‑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)
                                        aaTATAAAATTTTGAATTGTTATTTTTGTGTGTTGTTTAAGTTTATAGGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                             <  1:1309864/100‑1 (MQ=255)
                                        aaTATAAAATATTGAGTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCATCCTGTATTAAAAt                                                             >  2:1576925/1‑100 (MQ=255)
                                        aaTATAAAATATTGAATTGTTATTTTTGTGTGTTCTTTAAGATAAAAAGCCGTATTTATAATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                             <  2:928879/100‑1 (MQ=255)
                                          tataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAATATTa                                                           <  2:1385387/100‑1 (MQ=255)
                                            taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTATATTAAAATTAat                                                         >  2:1504539/1‑100 (MQ=255)
                                             aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                        <  1:2744427/100‑1 (MQ=255)
                                             aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAACATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                        >  2:3480445/1‑100 (MQ=255)
                                              aaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAatat                                                       <  1:183699/100‑1 (MQ=255)
                                                  atTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaa                                                   <  1:2730547/100‑1 (MQ=255)
                                                   ttGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaa                                                   >  1:1377283/1‑99 (MQ=255)
                                                    tGAATTGTTATTTTTGTGTGTTGTTTAAAATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                 <  2:2667603/100‑1 (MQ=255)
                                                     gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                                <  1:1741658/100‑1 (MQ=255)
                                                     gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                                <  1:1457840/100‑1 (MQ=255)
                                                     gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAACGCAGGCTGTATTAAAATTAATATTCAAAGc                                                >  2:1369333/1‑100 (MQ=255)
                                                      aaTTTTTATTTTTGTGTGTTGTTTAAGATATAAAGCTGTATTTATTATTAGGGCTTTAATTAATAAAGGGCAGGCGGTTTTAATATTAATATTCAAAg                                                 <  2:848374/98‑1 (MQ=255)
                                                      aaTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCa                                               >  1:2975562/1‑100 (MQ=255)
                                                      aaTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCa                                               >  1:2445185/1‑100 (MQ=255)
                                                       aTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCAt                                              >  2:2046988/1‑100 (MQ=255)
                                                        ttGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATa                                             >  1:570909/1‑100 (MQ=255)
                                                          gttgttTTTGTGTGTTGTGTAAAATAAAAAGCCGTATTTATTATTACGTCTTTAATTAAAAAAAGGCAGGCTGTATTAAAATTAATATTCaaa                                                  <  2:3699688/93‑1 (MQ=255)
                                                          gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                           >  2:177697/1‑100 (MQ=255)
                                                          gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                           >  2:1990484/1‑100 (MQ=255)
                                                           ttATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                          <  2:2430702/100‑1 (MQ=255)
                                                           ttATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                          <  2:3379354/100‑1 (MQ=255)
                                                            tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                         >  1:437589/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)
                                                                         tgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAc                             >  2:1189246/1‑99 (MQ=255)
                                                                          gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                           >  1:542197/1‑100 (MQ=255)
                                                                          gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                           >  1:1362585/1‑100 (MQ=255)
                                                                          gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                           <  2:3099271/100‑1 (MQ=255)
                                                                           tttattaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                          >  1:3583058/11‑100 (MQ=25)
                                                                           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                        >  2:1763361/1‑100 (MQ=255)
                                                                             tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                        >  1:1096559/1‑100 (MQ=255)
                                                                             tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATAGTCAAAGCATAAACCGATCGCCAATACCGGtt                        <  1:2782372/100‑1 (MQ=255)
                                                                               aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                      <  2:2694598/100‑1 (MQ=255)
                                                                               aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                      <  2:1102960/100‑1 (MQ=255)
                                                                               aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                      <  1:2343850/100‑1 (MQ=255)
                                                                                gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                      <  2:618896/99‑1 (MQ=255)
                                                                                 aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                    <  2:4124461/100‑1 (MQ=255)
                                                                                  tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                   <  2:2610278/100‑1 (MQ=255)
                                                                                  tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                   >  1:3015112/1‑100 (MQ=255)
                                                                                  attaaaGCCGTAATAATAAATACGGCTTTAATTAATACAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                   >  2:2216471/4‑100 (MQ=37)
                                                                                  attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                   >  1:2999198/4‑100 (MQ=37)
                                                                                   aaaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTa                  >  2:229280/1‑100 (MQ=255)
                                                                                    taaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                 >  2:2277106/2‑100 (MQ=37)
                                                                                     aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt                >  1:2206954/1‑100 (MQ=37)
                                                                                         ccGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaaa            >  2:3897541/1‑100 (MQ=255)
                                                                                          cGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaaa            <  2:4086079/99‑1 (MQ=255)
                                                                                           gTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCATTACTGaaacc          <  1:3604395/100‑3 (MQ=255)
                                                                                            tatttattATTACGGCTTTTATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATgg         >  2:2122228/1‑100 (MQ=255)
                                                                                                  ttattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATAGCGt      >  2:3520834/1‑97 (MQ=255)
                                                                                                   tattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgggc  >  2:3810730/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGC  >  NC_000913/720910‑721108

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: