Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,011 T→A 15.7% 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,0110TA15.7% 362.2 / 53.5 139intergenic (+171/+45)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (45/72);  new base A (10/12);  total (55/84)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 148 total aligned reads displayed.

CAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGC  >  NC_000913/720912‑721108
                                                                                                   |                                                                                                 
cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                                                                                   >  1:3688690/1‑100 (MQ=255)
cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                                                                                   <  2:150535/100‑1 (MQ=255)
cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                                                                                   <  2:2801879/100‑1 (MQ=255)
   tCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAc                                                                                                <  2:2056249/100‑1 (MQ=255)
    cTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                                                                                               <  2:1151129/100‑1 (MQ=255)
    cTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                                                                                               <  1:3815168/100‑1 (MQ=255)
        atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                           <  2:2718972/100‑1 (MQ=255)
        atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                           <  1:175456/100‑1 (MQ=255)
         tttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttt                                                                                          <  1:1050899/100‑1 (MQ=38)
          ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                         <  2:1972432/100‑2 (MQ=38)
           tatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaa                                                                                        <  2:1679223/100‑1 (MQ=255)
            atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTAGTTTAAGATATAAAGCCGTATTTATTATTATGGCTtta                                                                                         >  1:2980570/1‑98 (MQ=255)
             tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttat                                                                                      <  2:540267/100‑5 (MQ=37)
              ttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatta                                                                                     <  2:1185407/100‑1 (MQ=255)
               taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                                                                                    >  1:3155778/1‑100 (MQ=255)
               taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                    <  2:2771065/100‑7 (MQ=37)
               taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                    <  2:3590021/100‑7 (MQ=37)
                attgatgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                                   <  2:90108/98‑1 (MQ=255)
                aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                                   <  2:3342864/100‑1 (MQ=255)
                aaTGAGGAGCTTTTTACTCAGTAATATAAAAAATGGAATTGTTATTTTTGTGTGTGGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                                   <  1:3245187/100‑8 (MQ=21)
                 atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                  <  2:3197637/100‑1 (MQ=255)
                 atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                  >  1:2546662/1‑100 (MQ=255)
                 atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                  <  1:1728791/100‑1 (MQ=255)
                 atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                  <  2:2546674/100‑1 (MQ=255)
                    atgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa                                                                               >  1:2581732/1‑100 (MQ=255)
                     tgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTCTTATTACGGCTTTAATTAATaaaa                                                                               >  1:3754029/1‑99 (MQ=255)
                     tgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaatagatc                                                                              >  1:1297859/1‑96 (MQ=255)
                        gCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                           <  2:2460327/100‑1 (MQ=255)
                         cTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag                                                                          >  2:290545/1‑100 (MQ=255)
                                cTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa                                                                    <  2:1119096/99‑1 (MQ=255)
                                cTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAt                                                                   >  1:2232432/1‑100 (MQ=255)
                                cTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAt                                                                   <  1:738800/100‑1 (MQ=255)
                                 tCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAtt                                                                  >  2:1245996/1‑100 (MQ=255)
                                  cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                                 >  2:2446576/1‑100 (MQ=255)
                                  cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                                 >  2:2446920/1‑100 (MQ=255)
                                  cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                                 >  2:1309675/1‑100 (MQ=255)
                                     tAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaaa                                                              <  2:1080694/100‑1 (MQ=255)
                                     tAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaaa                                                              <  1:979386/100‑1 (MQ=255)
                                      aaTATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                                 <  1:152874/96‑1 (MQ=255)
                                       atatAAAATATTGAATTGTTATTTTTTTGTTTTGTTTAAAATAAAAAGCTGTATTTATTATTACGCCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                             <  1:1563410/99‑1 (MQ=255)
                                       atatAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAtt                                                            <  2:1026794/100‑1 (MQ=255)
                                       atatAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAtt                                                            >  2:1551012/1‑100 (MQ=255)
                                       atatAAAATATTGAATTGTTATTTTTGTGTGTTGCTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAATGAAATAATGTATTaaaact                                                            >  1:3171001/1‑98 (MQ=255)
                                        gccaacatatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAACAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTa                                                           <  1:2447824/94‑1 (MQ=255)
                                         cttccgatctGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTaa                                                          <  2:2054091/91‑1 (MQ=255)
                                         ataAAATATTGACTTCTTATTTTTGTGTGTTGTTTAACATACAAAGCCGTATTTATTATTACGGTTTTAATTAATAAAAGTGAGGCTGTATTAAAATTaa                                                          >  1:1318935/1‑100 (MQ=255)
                                          taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                         <  2:1669943/100‑1 (MQ=255)
                                          taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                         <  2:2626091/100‑1 (MQ=255)
                                           aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGTCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                        <  1:296554/100‑1 (MQ=255)
                                           aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                        >  2:1339915/1‑100 (MQ=255)
                                           aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                        >  1:831636/1‑100 (MQ=255)
                                           aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                        >  1:3246181/1‑100 (MQ=255)
                                             aaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATAtt                                                      >  1:1070043/1‑100 (MQ=255)
                                               tatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCa                                                    <  1:2100905/100‑1 (MQ=255)
                                                 ttGAATTGTTATTTTTTTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTATGTTTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaaa                                                  >  1:199747/1‑100 (MQ=255)
                                                  tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                 <  2:2835220/100‑1 (MQ=255)
                                                  tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                 <  2:2062985/100‑1 (MQ=255)
                                                    aatttctATTTTTGTTTATTGTTAAAGATAAAAAGCCGAATTTACTATTACGACTTTAATTAATAAAAGGCAGTCTGTATTAAAATTAATATTCAAAg                                                 <  1:2939660/92‑1 (MQ=255)
                                                    aaTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                                <  1:2794293/99‑1 (MQ=255)
                                                      ttGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCAt                                              >  1:921035/1‑99 (MQ=255)
                                                           aTTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCg                                        <  2:436590/100‑1 (MQ=255)
                                                            ttttttTGTGTTGTTTAAAATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                         <  1:2366110/98‑1 (MQ=255)
                                                            tttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGa                                       <  2:3804109/100‑1 (MQ=255)
                                                            tttttGAGTGTTGGTTAAGATAAAAAGGCGTATTTATTATTACGGCTTGAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAACGa                                       <  1:1850343/100‑1 (MQ=255)
                                                              tttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAtn                                     <  2:1898860/100‑2 (MQ=255)
                                                                 gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAg                                    <  1:2397361/98‑1 (MQ=255)
                                                                 gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                                  <  2:3814963/100‑1 (MQ=255)
                                                                 gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                                  <  1:85930/100‑1 (MQ=255)
                                                                   gtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCaa                                >  1:503401/1‑100 (MQ=255)
                                                                    tgttgtTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAAt                               <  1:210895/100‑1 (MQ=255)
                                                                      ttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAc                             <  2:2221638/100‑1 (MQ=255)
                                                                        gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                           >  1:3778573/1‑100 (MQ=255)
                                                                         tttattaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                          >  2:703858/11‑100 (MQ=25)
                                                                         tttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                          >  1:2237979/1‑100 (MQ=255)
                                                                          ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                         >  2:221269/1‑100 (MQ=255)
                                                                            aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                       >  1:3015393/1‑100 (MQ=255)
                                                                            aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                       >  2:2754069/1‑100 (MQ=255)
                                                                            aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                       >  2:512413/1‑100 (MQ=255)
                                                                             aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                      <  2:1823829/100‑1 (MQ=255)
                                                                              gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCTAATATCGGTTTCa                     >  1:3471942/1‑100 (MQ=255)
                                                                              gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                     >  2:2303013/1‑100 (MQ=255)
                                                                              gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                     <  1:1245989/100‑1 (MQ=255)
                                                                               aattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                     >  1:1472886/5‑99 (MQ=37)
                                                                               aTAAAAAGCCGTATTTATTATTCCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTTAAAGCATAAACCGATAGCCAATACCGGTTTCAg                    <  2:921042/100‑1 (MQ=255)
                                                                                tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                   >  2:3385213/1‑100 (MQ=255)
                                                                                tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                   <  2:1426109/100‑1 (MQ=255)
                                                                                attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                   >  2:3101077/4‑100 (MQ=37)
                                                                                  taaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                 >  1:165849/2‑100 (MQ=37)
                                                                                  aaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                 <  1:2861600/100‑1 (MQ=255)
                                                                                  aaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                 >  2:1458491/1‑100 (MQ=255)
                                                                                  aaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                 <  2:1289358/100‑1 (MQ=255)
                                                                                   aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGAATAAACCGATAGCCAATACCGGTTTCAGTAAt                >  1:1186892/1‑100 (MQ=37)
                                                                                    aaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg               <  2:3733476/100‑1 (MQ=255)
                                                                                    aaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGACTGTATAAAAATTAATATTCAAAACATAAAACCATAACCAATAACGGTTTCAGCAATg               >  1:1315174/1‑100 (MQ=255)
                                                                                      gCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaa             <  2:2258819/100‑1 (MQ=255)
                                                                                        cGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAAt           <  1:3584641/100‑1 (MQ=255)
                                                                                         gTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATg          <  1:259925/100‑1 (MQ=255)
                                                                                            tttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTAAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCg       >  1:437336/1‑100 (MQ=255)
                                                                                             ttattattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATGCCGGTTTCAGTAACGAAATAGc        >  1:3114054/1‑98 (MQ=255)
                                                                                               attattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgg    <  1:3486777/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
CAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGC  >  NC_000913/720912‑721108

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: