Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,006 T→A 18.1% intergenic (+166/+50) 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,0060TA18.1% 367.3 / 65.0 150intergenic (+166/+50)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/11);  total (72/78)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 157 total aligned reads displayed.

TGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGC  >  NC_000913/720908‑721102
                                                                                                  |                                                                                                
tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattta                                                                                                 <  1:3046891/100‑1 (MQ=255)
 gATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttan                                                                                                <  1:3559948/100‑2 (MQ=255)
  aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                               <  1:2323442/100‑1 (MQ=255)
      aTCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatta                                                                                           <  2:1114260/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)
           aaTTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTAAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCt                                                                                      <  2:2051546/100‑1 (MQ=255)
             tttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                    >  1:1768455/1‑100 (MQ=255)
             tttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttt                                                                                    <  1:3159952/100‑1 (MQ=38)
              ttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                   <  2:685491/100‑2 (MQ=38)
                 tttaATGATTAGTTTTTTACTTAGTAATTTAAAATATAGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                   <  2:3268449/97‑2 (MQ=25)
                 tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                <  1:951674/100‑1 (MQ=255)
                 tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttan                                                                                <  2:1559050/100‑5 (MQ=37)
                  ttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatta                                                                               <  2:3393693/100‑1 (MQ=255)
                   taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                                                                              <  2:3681405/100‑1 (MQ=255)
                   taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                              <  1:2045797/100‑7 (MQ=25)
                    aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCTTATTTATTATTACGGCTTtaattaat                                                                             <  1:2513465/100‑1 (MQ=255)
                    aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                             <  1:3155074/100‑1 (MQ=255)
                    aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                             <  2:3321099/100‑8 (MQ=25)
                    aaTGATGAGCTGCTTACTTAGTACTATACAATATTGAATTGTTATTTTTGTGTGTTTTTTAAGATAAAAAGCCGTATTTATTTTTAAGGCTtta                                                                                   >  2:2353628/1‑94 (MQ=255)
                     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)
                       gatgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatcttaaa                                                                          <  2:1751511/100‑11 (MQ=21)
                         tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCAGTATTTATTATTACGGCTTTAATTAATAAAag                                                                        >  2:56313/1‑100 (MQ=255)
                          gaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                       >  1:2474767/1‑100 (MQ=255)
                           aGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTTTTTTTTGTTTAATATAAAAAGCCGTATTTATTATTAATGCTTTAATTAATAAAag                                                                        >  2:2146071/1‑98 (MQ=255)
                           aGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                       <  2:2952240/99‑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)
                                  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)
                                     ttaGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAt                                                             <  2:3965297/97‑1 (MQ=255)
                                      taGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGTCTGTATTa                                                           >  2:1728212/2‑100 (MQ=255)
                                      taGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                           >  1:576273/2‑100 (MQ=255)
                                        gTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaa                                                         <  2:2081992/100‑1 (MQ=255)
                                        gTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaa                                                         <  1:568499/100‑1 (MQ=255)
                                          aaTATAAAATTTTGAATTGTTATTTTTGTGTGTTGTTTAAGTTTATAGGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                       <  1:1309864/100‑1 (MQ=255)
                                          aaTATAAAATATTGAATTGTTATTTTTGTGTGTTCTTTAAGATAAAAAGCCGTATTTATAATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                       <  2:928879/100‑1 (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)
                                                aaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAatat                                                 >  1:725721/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)
                                                       gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                          >  2:1411619/1‑100 (MQ=255)
                                                       gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                          <  1:1457840/100‑1 (MQ=255)
                                                       gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                          <  1:1741658/100‑1 (MQ=255)
                                                       gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCATATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAACATTAATATTCAAAGc                                          >  2:3562876/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)
                                                            gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                     >  2:177697/1‑100 (MQ=255)
                                                             ttATTTTTTTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                    <  2:1983621/100‑1 (MQ=255)
                                                              tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                   <  2:1661200/100‑1 (MQ=255)
                                                              tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGAAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                    >  2:463741/1‑99 (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:542197/1‑100 (MQ=255)
                                                                            gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                     <  2:3099271/100‑1 (MQ=255)
                                                                            gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                     >  1:1362585/1‑100 (MQ=255)
                                                                             tttattaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                    >  2:1249096/11‑100 (MQ=25)
                                                                             tttAAGATAAAAAGCCGTATTTATTATTACGACTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgn                    <  2:305770/100‑2 (MQ=255)
                                                                              ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                   >  2:2341995/1‑100 (MQ=255)
                                                                               tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                  >  1:1096559/1‑100 (MQ=255)
                                                                               tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATAGTCAAAGCATAAACCGATCGCCAATACCGGtt                  <  1:2782372/100‑1 (MQ=255)
                                                                                attaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                 >  2:2627164/8‑100 (MQ=25)
                                                                                 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             >  1:3485727/1‑100 (MQ=255)
                                                                                    tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt             <  2:2610278/100‑1 (MQ=255)
                                                                                    tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt             >  1:3015112/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)
                                                                                      taaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa           >  2:1285925/2‑100 (MQ=37)
                                                                                      taaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa           >  2:2277106/2‑100 (MQ=37)
                                                                                      aaaaGCCGTATTTATTATTACGTCTTTAATTAACAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATCGCCAATACCGATTTCAGt             >  2:2108212/1‑98 (MQ=38)
                                                                                       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)
                                                                                             gTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATg    >  1:3884857/1‑100 (MQ=255)
                                                                                              tatttattATTACGGCTTTTATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATgg   >  2:2122228/1‑100 (MQ=255)
                                                                                                  |                                                                                                
TGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGC  >  NC_000913/720908‑721102

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