Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,004 T→A 15.0% intergenic (+164/+52) 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,0040TA15.0% 530.1 / 73.3 192intergenic (+164/+52)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (60/103);  new base A (19/10);  total (79/113)
Fisher's exact test for biased strand distribution p-value = 6.92e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 205 total aligned reads displayed.

AATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCG  >  NC_000913/720906‑721103
                                                                                                  |                                                                                                   
attGATTAATTTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAATATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatt                                                                                                    >  1:1765470/3‑100 (MQ=255)
 aTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattt                                                                                                   <  1:2410900/100‑1 (MQ=255)
  tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattta                                                                                                  <  2:722867/100‑1 (MQ=255)
  tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattta                                                                                                  <  1:655239/100‑1 (MQ=255)
  tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattta                                                                                                  <  1:88461/100‑1 (MQ=255)
    aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                                >  2:3475751/1‑100 (MQ=255)
     tCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttatta                                                                                               <  1:2733042/100‑1 (MQ=255)
         tCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAc                                                                                           <  1:810327/100‑1 (MQ=255)
          cTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                                                                                          >  2:3395616/1‑100 (MQ=255)
            tAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                        <  1:1850367/100‑1 (MQ=255)
             aaTTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                                                                                       <  1:1758549/100‑1 (MQ=255)
             aaTTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                                                                                       <  2:324204/100‑1 (MQ=255)
              atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                      <  1:427635/100‑1 (MQ=255)
               tttatttaATGATGAGTTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                      <  1:3347476/99‑1 (MQ=255)
                ttatttaATGATGAGTTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                     <  1:2689409/99‑1 (MQ=255)
                ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTtta                                                                                    <  1:2939502/100‑1 (MQ=255)
                ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTtta                                                                                    <  1:2304414/100‑1 (MQ=255)
                ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                    <  2:1726321/100‑2 (MQ=38)
                 tatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaa                                                                                   <  1:2634441/100‑1 (MQ=255)
                  atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                                                                                  >  2:1560693/1‑100 (MQ=255)
                  atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                                                                                  <  1:1531088/100‑1 (MQ=255)
                  atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttttta                                                                                  <  2:2012280/100‑4 (MQ=38)
                   tttaATGATGAGTTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                 <  1:2549216/100‑1 (MQ=255)
                   tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                 <  1:2158638/100‑1 (MQ=255)
                    ttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatta                                                                                <  2:853977/100‑1 (MQ=255)
                     taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                                                                               >  1:1536380/1‑100 (MQ=255)
                     taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                                                                               >  2:1669077/1‑100 (MQ=255)
                     taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                               <  2:3846225/100‑7 (MQ=37)
                     taATGATGAGCTTTTTACTCAGTAAGATAAAATATTGAATTGTTATTTTTGTGGGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                                                                               <  1:1154126/100‑1 (MQ=255)
                      aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTTTTATTTTTGTGTGTTGTTTAAAATAAAAAGCCGTATTTATTATTACGGCTttaattaa                                                                               <  1:1990963/99‑1 (MQ=255)
                      aaTGATAAGCTTTTTACTCAGAAAGATAAAATATTGAATAGTTATTATTGTGTGTTATTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttat                                                                                 <  1:2974784/97‑5 (MQ=21)
                        tgatgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                             <  1:501410/99‑1 (MQ=255)
                           tgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag                                                                         >  1:3795198/1‑100 (MQ=255)
                            gaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                        <  1:491806/100‑1 (MQ=255)
                             aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc                                                                       <  1:1033302/100‑1 (MQ=255)
                             aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc                                                                       <  2:2785441/100‑1 (MQ=255)
                             aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc                                                                       >  1:1021181/1‑100 (MQ=255)
                             aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTCTAATTAATAAAaggc                                                                       >  2:3813923/1‑100 (MQ=255)
                              gCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGAGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                      >  2:362427/1‑100 (MQ=255)
                                tttttACTCAGTAATATAAAATATTGAATTGTTATTTTTTTTTTTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag                                                                     >  1:1166532/1‑99 (MQ=255)
                                   ttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTg                                                                 <  1:1798970/100‑1 (MQ=255)
                                   ttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTg                                                                 <  1:135811/100‑1 (MQ=255)
                                    tACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt                                                                >  2:3326696/1‑100 (MQ=255)
                                     aCTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa                                                               <  2:3427493/100‑1 (MQ=255)
                                      ctaaGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAGAGCGGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTgga                                                               <  1:3567694/96‑3 (MQ=255)
                                           tAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaaa                                                         <  1:2895244/100‑1 (MQ=255)
                                            aaTATAAAATATTGAATTGTTATTTTTTTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                        <  1:1199778/100‑1 (MQ=255)
                                            aaTATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                        <  2:3344263/100‑1 (MQ=255)
                                             atatAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAtt                                                       <  1:3737405/100‑1 (MQ=255)
                                                 aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAACTTAata                                                   >  1:3826201/1‑100 (MQ=255)
                                                 aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                   >  1:1396755/1‑100 (MQ=255)
                                                  aaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTCTTAAAATTAatat                                                  >  2:1286437/1‑100 (MQ=255)
                                                  aaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAatat                                                  <  2:3540501/100‑1 (MQ=255)
                                                    atatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATAtt                                                 >  1:546989/1‑99 (MQ=255)
                                                    atatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTc                                                <  2:2431802/100‑1 (MQ=255)
                                                      atTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaa                                              <  2:1130074/100‑1 (MQ=255)
                                                       ttGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaaa                                             >  2:1752414/1‑100 (MQ=255)
                                                         gAATTGTTATTTTTTTGTTTTTTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTTAAAGc                                           >  1:3303588/1‑100 (MQ=255)
                                                         gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                           >  2:2476369/1‑100 (MQ=255)
                                                         gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                           <  2:46107/100‑1 (MQ=255)
                                                          aaTTGTTATTTTTGTGTGTTGTTTAATATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCa                                          <  2:3739495/100‑1 (MQ=255)
                                                          aaTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCa                                          >  2:3411987/1‑100 (MQ=255)
                                                               ttATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                     >  1:3213147/1‑100 (MQ=255)
                                                               ttATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                     <  1:316234/100‑1 (MQ=255)
                                                                tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                    <  1:1309310/100‑1 (MQ=255)
                                                                tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                    <  1:1053584/100‑1 (MQ=255)
                                                                    tttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCCGGCTGTATTAAAATTAATATTCAACGCATAAACCGATa                                <  1:2389714/100‑1 (MQ=255)
                                                                    tttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                <  2:272607/100‑1 (MQ=255)
                                                                      tgtgtgTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGc                              <  1:1722634/100‑1 (MQ=255)
                                                                      tgtgtgTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGc                              <  1:3783294/100‑1 (MQ=255)
                                                                       gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                             <  2:1881122/100‑1 (MQ=255)
                                                                           gttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAAt                          >  2:2571745/1‑99 (MQ=255)
                                                                              gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                      >  2:2942252/1‑100 (MQ=255)
                                                                              gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                      >  1:2090377/1‑100 (MQ=255)
                                                                                ttattaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                    >  2:2743897/10‑100 (MQ=25)
                                                                                 tattaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                   >  2:1447965/9‑100 (MQ=25)
                                                                                 tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                   >  2:3587638/1‑100 (MQ=255)
                                                                                  attaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                  >  2:2277500/8‑100 (MQ=25)
                                                                                  aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                  >  2:179818/1‑100 (MQ=255)
                                                                                   aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                 <  1:3428671/100‑1 (MQ=255)
                                                                                   aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                 <  1:325117/100‑1 (MQ=255)
                                                                                     aattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATATCCAATACTGGTTTTAg               >  1:2262545/5‑100 (MQ=25)
                                                                                     aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg               <  2:2141032/100‑1 (MQ=255)
                                                                                      tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt              <  2:54928/100‑1 (MQ=255)
                                                                                      tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt              <  1:1507825/100‑1 (MQ=255)
                                                                                      tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt              >  1:595374/1‑100 (MQ=255)
                                                                                      attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt              >  2:1963358/4‑100 (MQ=37)
                                                                                      attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt              >  1:939516/4‑100 (MQ=37)
                                                                                      attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTAAAAAAATAAACGGATAGCCAATACCGCttt                  >  1:1883509/4‑96 (MQ=25)
                                                                                         aaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa            >  1:757266/1‑99 (MQ=255)
                                                                                         aaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt           <  2:390423/100‑1 (MQ=255)
                                                                                         aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt           >  1:3067794/1‑100 (MQ=37)
                                                                                          aaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtttatgctttg          <  2:1597561/100‑9 (MQ=255)
                                                                                          aaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg          <  2:2174956/100‑1 (MQ=255)
                                                                                             ccGTATTTATTATTACGGCTTTAATTAATATAAGGCAGGCTATATTAAAATCAATATTCAAAGCATAACCCGATAGCCAATACCGGTTTCAGTAATTaaa       <  1:3798102/100‑1 (MQ=37)
                                                                                             ccGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaaa       <  2:3049275/100‑1 (MQ=255)
                                                                                              cGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATCCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAAt      >  1:2482679/1‑100 (MQ=255)
                                                                                                tatttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATgg    <  1:1479524/100‑1 (MQ=255)
                                                                                                 atttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGc   >  2:3265333/1‑100 (MQ=255)
                                                                                                  tttattATTACGGCTTTAATTAATAAAAGGCAGGCTTTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCATTACTGaaaccggt  <  1:2139288/100‑6 (MQ=38)
                                                                                                  |                                                                                                   
AATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCG  >  NC_000913/720906‑721103

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