Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,011 T→A 15.0% 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.0% 386.4 / 48.9 145intergenic (+171/+45)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (60/63);  new base A (6/16);  total (66/79)
Fisher's exact test for biased strand distribution p-value = 6.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 155 total aligned reads displayed.

CAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCG  >  NC_000913/720912‑721109
                                                                                                   |                                                                                                   
cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAAATGTTATTTTTGTGTATTGTTTAAGATAAAAAGCCGTATttattat                                                                                                     >  1:3435794/1‑100 (MQ=255)
  aTCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat‑ta                                                                                                  >  2:2217342/1‑100 (MQ=255)
    cTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACg                                                                                                <  1:2019150/100‑1 (MQ=255)
      tAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGc                                                                                              >  1:270237/1‑100 (MQ=255)
       aaTTTATTTAATGATGAGCTTTTTACTCTGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAA‑TACGGCt                                                                                             <  2:623761/100‑1 (MQ=255)
       aaTTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCt                                                                                             >  1:332004/1‑100 (MQ=255)
        atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCtt                                                                                            >  1:1747400/1‑100 (MQ=255)
         tttatttaATGATGAGCTTTTTACTCATTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCttt                                                                                           <  1:465410/100‑1 (MQ=255)
         tttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAATAAGCCGTAATAATAAA‑TACGGCttt                                                                                           <  2:243666/100‑1 (MQ=255)
         tttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCttt                                                                                           >  2:1462908/1‑100 (MQ=255)
         tttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAA‑TACGGCttt                                                                                           <  1:3757743/100‑1 (MQ=38)
           tatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAA‑TACGGCttttt                                                                                         <  1:3405804/100‑3 (MQ=38)
            atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAA‑TACGGCttttta                                                                                        <  2:2902671/100‑4 (MQ=38)
             tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTttaatt                                                                                       <  2:552828/100‑1 (MQ=255)
             tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAA‑TACGGCtttttat                                                                                       <  1:192279/100‑5 (MQ=37)
               taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTttaattaa                                                                                     <  2:1594272/100‑1 (MQ=255)
               taATGATAAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAA‑TACGGCtttttatct                                                                                     <  2:336991/100‑7 (MQ=25)
                aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTtaattaat                                                                                    <  2:122267/100‑1 (MQ=255)
                aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTtaattaat                                                                                    >  2:3040989/1‑100 (MQ=255)
                aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAA‑TACGGCtttttatctt                                                                                    <  2:2200920/100‑8 (MQ=37)
                 atgatgAGTTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGTCGTAATAATAAA‑TACGGCtttttatctta                                                                                   <  2:3479168/100‑9 (MQ=25)
                  tgatgaGTTTTTTATTTAGTAATATTAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCTTA‑ATATTATATACTGCtttttatctt                                                                                    <  2:1340121/98‑8 (MQ=255)
                    atgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATaaaa                                                                                <  1:1375903/100‑1 (MQ=255)
                    atgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGATTTAATtaata                                                                                   <  2:2367926/97‑1 (MQ=255)
                      gaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATaaaa                                                                                >  2:1999944/1‑98 (MQ=255)
                       aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTAATTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAaggc                                                                             >  2:2083435/1‑100 (MQ=255)
                        gCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAaggca                                                                            >  2:144659/1‑100 (MQ=255)
                        gCTTTGTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAaggca                                                                            >  1:2307998/1‑100 (MQ=255)
                          tttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAaggcagg                                                                          <  2:1892280/100‑1 (MQ=255)
                           ttttACTCAGTAATATAAAATATTGAATTTTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAaggcaggc                                                                         <  2:2573265/100‑1 (MQ=255)
                           ttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAaggcaggc                                                                         <  1:1057161/100‑1 (MQ=255)
                           ttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAaggcaggc                                                                         >  2:2141197/1‑100 (MQ=255)
                              tACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGt                                                                      >  2:1505429/1‑100 (MQ=255)
                               aCTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTa                                                                     >  1:1147343/1‑100 (MQ=255)
                               aCTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTa                                                                     >  1:2946096/1‑100 (MQ=255)
                                  cAGTACTATAAAATATTGAATTGTTATTATTGTGTGTAATAGAAGATAAAAAGGCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGAAGGCTGTAtt                                                                   >  2:1038463/1‑99 (MQ=255)
                                  cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                                  >  2:3383649/1‑100 (MQ=255)
                                  cAGTAATATAAAATATTGAATTGTTAATTTTGTGTGTTGTTTAAGATAAAAAACCGTATTTATTAT‑TATGCCTTTAATTAATACTAGACAGGCGGTAt                                                                    >  2:2800357/1‑98 (MQ=255)
                                   aGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTaa                                                                 <  2:3236244/100‑1 (MQ=255)
                                     tAATATAAAATATTGAATTGTTATTTTTGTGTGGTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                                  <  1:3656931/97‑1 (MQ=255)
                                        tatataATATTGAATTGTTATTTTTGTGTGTTGTTTATGTTAAAAAGTCGTATTTATTAT‑TACGTTTTTAATTAATAAATGGTAGGTTGTATTAAAAtt                                                             <  2:1782084/99‑1 (MQ=255)
                                         ataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTaa                                                           <  1:773042/100‑1 (MQ=255)
                                         ataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTaa                                                           >  2:839158/1‑100 (MQ=255)
                                          taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                          >  2:314013/1‑100 (MQ=255)
                                          taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                          >  1:2235002/1‑100 (MQ=255)
                                          taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                          <  2:2608063/100‑1 (MQ=255)
                                                atTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaa                                                    >  2:1855888/1‑100 (MQ=255)
                                                atTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaa                                                    >  1:3834428/1‑100 (MQ=255)
                                                 ttGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaaa                                                   <  1:35307/100‑1 (MQ=255)
                                                 ttGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaaa                                                   <  2:2610352/100‑1 (MQ=255)
                                                  ttaaTTGTTATTTTTGTGTGTTGTTTATTTTAAAAATCCTTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                  <  1:1095738/98‑1 (MQ=255)
                                                  tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                  <  1:2499689/100‑1 (MQ=255)
                                                   gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACAGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                  <  2:1194103/99‑1 (MQ=255)
                                                   gAATTGTTATTTTTGTGTGTTATTTAAAATAAAAAGACGTATTTATTAA‑TACGGCTTTAATTAATAAAAGGCATGTTGTATTAAAATTAATATTCAAAgt                                                 <  2:2293505/100‑2 (MQ=255)
                                                       tGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATAAAAATTAATATTCAAAGCATaa                                             >  1:1574023/1‑100 (MQ=255)
                                                        gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                            >  1:818872/1‑100 (MQ=255)
                                                         ttATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                           >  1:3739827/1‑100 (MQ=255)
                                                            tttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGa                                        >  2:3552389/1‑100 (MQ=255)
                                                             ttttttttGTTTTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAATCATAAACCGa                                        >  2:263304/1‑99 (MQ=255)
                                                               ttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAg                                     <  2:1308639/100‑1 (MQ=255)
                                                               ttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAg                                     <  1:2475247/100‑1 (MQ=255)
                                                               ttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAg                                     <  2:3452854/100‑1 (MQ=255)
                                                                 gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTACAATTAATATTCAAAGCATAAACCGATaac                                    >  2:3072462/1‑97 (MQ=255)
                                                                 gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcn                                   <  1:35990/100‑2 (MQ=255)
                                                                 gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                                   <  2:1366723/100‑1 (MQ=255)
                                                                 gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                                   <  1:3656468/100‑1 (MQ=255)
                                                                    tgttgtTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGGTattggct                                >  1:321519/1‑94 (MQ=255)
                                                                     gttgttTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATa                               >  2:3159189/1‑100 (MQ=255)
                                                                        gttTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                            >  2:2088340/1‑100 (MQ=255)
                                                                        gttTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATCAATATTCAAAGCATAAACCGATAGCCAATACCg                            >  1:1527865/1‑100 (MQ=255)
                                                                          ttAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                           >  2:981837/1‑99 (MQ=255)
                                                                          ttAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                          <  1:3073769/100‑1 (MQ=255)
                                                                           tAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                         <  2:2005329/100‑1 (MQ=255)
                                                                           tAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                         <  2:1573160/100‑1 (MQ=255)
                                                                              gATAAAAAGCGGTATGTATTAT‑TACGGCTTTAAGTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCAGAAACCGATAGCCAATAGCGGtt                         <  2:2457737/97‑1 (MQ=255)
                                                                               aTAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGTTGTATTAAAATTAATATTCAAAGCAGAAACCGATAGCCAATACCGGTTTCAg                     <  1:108634/100‑1 (MQ=255)
                                                                                tAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                     >  2:2542637/1‑99 (MQ=255)
                                                                                tAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                    >  2:2116734/1‑100 (MQ=255)
                                                                                attaaaGCCGTAATAATAAA‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                    >  2:405226/4‑100 (MQ=37)
                                                                                  aaaaGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                  >  1:3465132/1‑100 (MQ=255)
                                                                                  aaaaGCCGTAATAATAAA‑TAGGGCTTTAATTAATAAAAGGCAGGCTGCATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAACGGTTTAAGTaa                  >  2:2395994/1‑100 (MQ=25)
                                                                                   aaaGCCGTAATAATAAA‑TACGGCTTTAATTAATAAAAGTCATGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt                 <  1:1058337/100‑1 (MQ=37)
                                                                                   aaaGCCGTAATAATAAA‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt                 >  1:1340922/1‑100 (MQ=37)
                                                                                   aaaGCCGTAATAATAAA‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt                 >  2:1058326/1‑100 (MQ=37)
                                                                                   aaaGCCGTAATAATAAA‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACAGATAGCCAATACCGGTTTCAGTAAt                 >  1:3385369/1‑100 (MQ=37)
                                                                                    aaGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg                <  1:514426/100‑1 (MQ=255)
                                                                                    aaGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg                <  1:2621422/100‑1 (MQ=255)
                                                                                     ccgatcttttattAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGa               <  1:382928/93‑1 (MQ=255)
                                                                                        cGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAAt            >  1:2792736/1‑100 (MQ=255)
                                                                                         gTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTTAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATg           >  2:3808546/1‑100 (MQ=255)
                                                                                          tatttattAT‑TACGGCTTTAATTAATAAAAGGCATGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATgg          <  1:3218632/100‑1 (MQ=255)
                                                                                           atttattAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATACAACGATAGACACTACCGGTTTCAGTAATTaaagagc         >  2:3304362/1‑96 (MQ=38)
                                                                                           atttattAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCCATAGCAAATACCGGTTTCAGTAATGAAATg           >  1:691703/1‑98 (MQ=255)
                                                                                            tttattAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaagatcgg        >  2:382922/1‑94 (MQ=255)
                                                                                              tattat‑taCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTg      >  1:1100901/1‑100 (MQ=255)
                                                                                               attat‑taCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgg     <  1:796959/100‑1 (MQ=255)
                                                                                               attat‑taCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgg     >  2:1499740/1‑100 (MQ=255)
                                                                                               attat‑taCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgg     <  2:2189784/100‑1 (MQ=255)
                                                                                                 tat‑taCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTggg    <  2:223431/99‑1 (MQ=255)
                                                                                                  at‑taCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCTATAGCCAATACCGGTATCAGTAATAAAATGGCGTgggcg  >  1:775733/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAT‑TACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCG  >  NC_000913/720912‑721109

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