Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,006 T→A 16.4% 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,0060TA16.4% 465.3 / 74.2 177intergenic (+166/+50)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (54/94);  new base A (16/13);  total (70/107)
Fisher's exact test for biased strand distribution p-value = 6.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 193 total aligned reads displayed.

TGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTG  >  NC_000913/720908‑721105
                                                                                                  |                                                                                                   
tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattta                                                                                                    <  2:1506985/100‑1 (MQ=255)
    cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTCTTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                                                                                >  1:1001212/1‑100 (MQ=255)
     aaTCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatt                                                                                               >  2:3225513/1‑100 (MQ=255)
       tCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAc                                                                                             <  2:2919337/100‑1 (MQ=255)
          tAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                          <  2:917596/100‑1 (MQ=255)
            atttctttATTGTTGATCTTTTTAGTGAGTAATACAAAATATTGAATTGTTATTTTTTTGTGTTTTTTAAGATAAAAAGCCGTATTTAGTATGACGGCtt                                                                                        <  2:2200680/100‑1 (MQ=255)
              ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                      <  2:2183808/100‑2 (MQ=38)
              ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                      <  1:2446325/100‑2 (MQ=38)
                atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaa                                                                                     >  2:3581553/1‑99 (MQ=255)
                atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttttta                                                                                    <  2:2711162/100‑4 (MQ=38)
                 tttaATGATGAGCTTTTTACTCAGTAATATAAAATATAGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttat                                                                                   <  2:3342550/100‑5 (MQ=37)
                  ttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatta                                                                                  <  1:1601029/100‑1 (MQ=255)
                  ttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatta                                                                                  <  1:2143353/100‑1 (MQ=255)
                  ttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTAGTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatc                                                                                  <  2:3357067/100‑6 (MQ=37)
                   taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                                                                                 <  1:619804/100‑1 (MQ=255)
                   taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                 <  2:3268696/100‑7 (MQ=37)
                    aaTGATGAGCTTTTTATTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                                <  2:2475879/100‑1 (MQ=255)
                    aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                                <  1:3597887/100‑8 (MQ=37)
                    aaTGATGAGCTTTTTACGCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                                <  1:1433842/100‑1 (MQ=255)
                     atgatgAGTTTTTTACTCAATAATATAAAATATTGAATTTTTATTTTTTTGTCTTGTTTAATATAAAAAGCCTTATTTATTATTACGGCTTTAATtaata                                                                               <  2:1489872/100‑1 (MQ=255)
                     atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                               >  2:1597461/1‑100 (MQ=255)
                     atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATGTTTGTGGGTTGTTTAAGATAAAAAGCCGTATTTATTATTACCGCTTTAATtaata                                                                               >  2:1004054/1‑100 (MQ=255)
                       gatgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaa                                                                             >  1:328875/1‑100 (MQ=255)
                         tgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag                                                                           >  1:1399882/1‑100 (MQ=255)
                          gaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                          >  2:2640022/1‑100 (MQ=255)
                           aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc                                                                         <  1:2024450/100‑1 (MQ=255)
                            gCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                        >  2:1182478/1‑100 (MQ=255)
                             cTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag                                                                       <  2:3528466/100‑1 (MQ=255)
                             cTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag                                                                       >  1:3731820/1‑100 (MQ=255)
                               ttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc                                                                     >  1:4149860/1‑100 (MQ=255)
                                 ttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTg                                                                   <  2:3050654/100‑1 (MQ=255)
                                  tACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt                                                                  <  2:1571448/100‑1 (MQ=255)
                                      cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                              >  2:2561740/1‑100 (MQ=255)
                                        gTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaa                                                             >  2:974802/1‑99 (MQ=255)
                                            tataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAttn                                                        <  1:1867405/100‑2 (MQ=255)
                                            tataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTa                                                        <  1:2191466/100‑1 (MQ=255)
                                             ataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTTAAATTa                                                        >  2:777105/1‑99 (MQ=255)
                                              taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                      >  1:394316/1‑100 (MQ=255)
                                               aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                     >  1:1791345/1‑100 (MQ=255)
                                               aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                     <  2:1806412/100‑1 (MQ=255)
                                                  atatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTc                                                  >  1:1063308/1‑100 (MQ=255)
                                                   tatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCa                                                 <  1:4036278/100‑1 (MQ=255)
                                                      tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                              <  1:1169228/100‑1 (MQ=255)
                                                      tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                              <  1:2252533/100‑1 (MQ=255)
                                                      tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                              <  1:2051533/100‑1 (MQ=255)
                                                        aaTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                             <  2:3751319/99‑1 (MQ=255)
                                                        aaTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCa                                            <  2:3175012/100‑1 (MQ=255)
                                                          ttGTTATTTTTGTGTGTTGTTTAATATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATa                                          <  2:572783/100‑1 (MQ=255)
                                                          ttGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATa                                          <  1:2096559/100‑1 (MQ=255)
                                                          ttGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATa                                          <  1:97850/100‑1 (MQ=255)
                                                              tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                      >  2:4107172/1‑100 (MQ=255)
                                                              tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                       >  2:1142299/1‑99 (MQ=255)
                                                                tttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGa                                    >  1:2741822/1‑100 (MQ=255)
                                                                 ttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAt                                   <  2:4159719/100‑1 (MQ=255)
                                                                  tttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                  <  1:3767696/100‑1 (MQ=255)
                                                                  tttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                  <  2:726559/100‑1 (MQ=255)
                                                                    tgtgtgTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGc                                <  1:3231979/100‑1 (MQ=255)
                                                                     gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                               <  2:123046/100‑1 (MQ=255)
                                                                      tgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCa                              >  1:1234035/1‑100 (MQ=255)
                                                                       gtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCaa                             <  1:2303004/100‑1 (MQ=255)
                                                                       gtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCa                              <  2:470797/99‑1 (MQ=255)
                                                                        tgttgtTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAAt                            >  2:2434085/1‑100 (MQ=255)
                                                                         gttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATa                           <  1:694976/100‑1 (MQ=255)
                                                                           tgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAcc                         <  1:1099705/100‑1 (MQ=255)
                                                                            gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATAATCCAAGCATAAACCGAGAGCCAATAc                          >  2:3882185/1‑98 (MQ=255)
                                                                             tttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACGGATAGCCAATa                           >  2:3434455/1‑96 (MQ=255)
                                                                             tttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                       <  1:3320257/100‑1 (MQ=255)
                                                                              ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                      <  2:1747269/100‑1 (MQ=255)
                                                                              ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAAGTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                      >  2:2227976/1‑100 (MQ=255)
                                                                              ttAAAATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                      <  2:3755729/100‑1 (MQ=255)
                                                                               tattaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                     >  1:2580404/9‑100 (MQ=25)
                                                                               tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                     <  1:4063129/100‑1 (MQ=255)
                                                                               tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                     <  2:2688891/100‑1 (MQ=255)
                                                                               tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                     <  1:3406692/100‑1 (MQ=255)
                                                                               tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                     <  1:760694/100‑1 (MQ=255)
                                                                                attaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                    >  1:2861011/8‑100 (MQ=25)
                                                                                aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                    <  1:1182479/100‑1 (MQ=255)
                                                                                 aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                   <  2:1335288/100‑1 (MQ=255)
                                                                                 aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                   >  1:909157/1‑100 (MQ=255)
                                                                                 aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                   >  2:3650647/1‑100 (MQ=255)
                                                                                   aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                 >  2:3120107/1‑100 (MQ=255)
                                                                                   aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                 <  1:2979591/100‑1 (MQ=255)
                                                                                   aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                 <  1:3085023/100‑1 (MQ=255)
                                                                                    tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                >  2:1960380/1‑100 (MQ=255)
                                                                                    attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCCAAGCATAAACCGATAGCCATGACCGGGTTCAGt                >  1:1989235/4‑100 (MQ=25)
                                                                                    attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                >  1:15302/4‑100 (MQ=37)
                                                                                    attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                >  2:1128429/4‑100 (MQ=37)
                                                                                     ttaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTa               >  2:451443/3‑100 (MQ=37)
                                                                                      taaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATACAAAGAATAAACCGGTAGCCAATACCGGTTTCAGTa               >  2:2822549/2‑99 (MQ=37)
                                                                                      aaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa              >  1:579478/1‑100 (MQ=255)
                                                                                      aaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa              >  1:470398/1‑100 (MQ=255)
                                                                                       aaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt             <  2:1833511/100‑1 (MQ=255)
                                                                                       aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt             >  1:880081/1‑100 (MQ=37)
                                                                                       aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt             >  2:1307748/1‑100 (MQ=37)
                                                                                       aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt             >  2:3100731/1‑100 (MQ=37)
                                                                                        aaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg            >  2:1076633/1‑100 (MQ=37)
                                                                                          gCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaa          <  1:1741926/100‑1 (MQ=255)
                                                                                             gTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATg       >  1:1674430/1‑100 (MQ=255)
                                                                                              tatttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATgg      >  1:808295/1‑100 (MQ=255)
                                                                                                tttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCg    >  1:35234/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTG  >  NC_000913/720908‑721105

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: