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.

CAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCG  >  NC_000913/720912‑721109
                                                                                                   |                                                                                                   
cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAAATGTTATTTTTGTGTATTGTTTAAGATAAAAAGCCGTATttattat                                                                                                     >  1:3435794/1‑100 (MQ=255)
  aTCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatta                                                                                                  >  2:2217342/1‑100 (MQ=255)
    cTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                                                                                                <  1:2019150/100‑1 (MQ=255)
      tAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                              >  1:270237/1‑100 (MQ=255)
       aaTTTATTTAATGATGAGCTTTTTACTCTGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCt                                                                                             <  2:623761/100‑1 (MQ=255)
       aaTTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                                                                                             >  1:332004/1‑100 (MQ=255)
        atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                            >  1:1747400/1‑100 (MQ=255)
         tttatttaATGATGAGCTTTTTACTCATTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                           <  1:465410/100‑1 (MQ=255)
         tttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAATAAGCCGTAATAATAAATACGGCttt                                                                                           <  2:243666/100‑1 (MQ=255)
         tttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                           >  2:1462908/1‑100 (MQ=255)
         tttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttt                                                                                           <  1:3757743/100‑1 (MQ=38)
           tatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttttt                                                                                         <  1:3405804/100‑3 (MQ=38)
            atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttttta                                                                                        <  2:2902671/100‑4 (MQ=38)
             tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                       <  2:552828/100‑1 (MQ=255)
             tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttat                                                                                       <  1:192279/100‑5 (MQ=37)
               taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                                                                                     <  2:1594272/100‑1 (MQ=255)
               taATGATAAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                     <  2:336991/100‑7 (MQ=25)
                aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                                    <  2:122267/100‑1 (MQ=255)
                aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                                    >  2:3040989/1‑100 (MQ=255)
                aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                                    <  2:2200920/100‑8 (MQ=37)
                 atgatgAGTTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGTCGTAATAATAAATACGGCtttttatctta                                                                                   <  2:3479168/100‑9 (MQ=25)
                  tgatgaGTTTTTTATTTAGTAATATTAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCTTAATATTATATACTGCtttttatctt                                                                                    <  2:1340121/98‑8 (MQ=255)
                    atgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa                                                                                <  1:1375903/100‑1 (MQ=255)
                    atgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGATTTAATtaata                                                                                   <  2:2367926/97‑1 (MQ=255)
                      gaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa                                                                                >  2:1999944/1‑98 (MQ=255)
                       aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTAATTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc                                                                             >  2:2083435/1‑100 (MQ=255)
                        gCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                            >  2:144659/1‑100 (MQ=255)
                        gCTTTGTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                            >  1:2307998/1‑100 (MQ=255)
                          tttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcagg                                                                          <  2:1892280/100‑1 (MQ=255)
                           ttttACTCAGTAATATAAAATATTGAATTTTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc                                                                         <  2:2573265/100‑1 (MQ=255)
                           ttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc                                                                         <  1:1057161/100‑1 (MQ=255)
                           ttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc                                                                         >  2:2141197/1‑100 (MQ=255)
                              tACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt                                                                      >  2:1505429/1‑100 (MQ=255)
                               aCTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa                                                                     >  1:1147343/1‑100 (MQ=255)
                               aCTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa                                                                     >  1:2946096/1‑100 (MQ=255)
                                  cAGTACTATAAAATATTGAATTGTTATTATTGTGTGTAATAGAAGATAAAAAGGCGTATTTATTATTACGGCTTTAATTAATAAAAGGAAGGCTGTAtt                                                                   >  2:1038463/1‑99 (MQ=255)
                                  cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                                  >  2:3383649/1‑100 (MQ=255)
                                  cAGTAATATAAAATATTGAATTGTTAATTTTGTGTGTTGTTTAAGATAAAAAACCGTATTTATTATTATGCCTTTAATTAATACTAGACAGGCGGTAt                                                                    >  2:2800357/1‑98 (MQ=255)
                                   aGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaa                                                                 <  2:3236244/100‑1 (MQ=255)
                                     tAATATAAAATATTGAATTGTTATTTTTGTGTGGTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                                  <  1:3656931/97‑1 (MQ=255)
                                        tatataATATTGAATTGTTATTTTTGTGTGTTGTTTATGTTAAAAAGTCGTATTTATTATTACGTTTTTAATTAATAAATGGTAGGTTGTATTAAAAtt                                                             <  2:1782084/99‑1 (MQ=255)
                                         ataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTaa                                                           <  1:773042/100‑1 (MQ=255)
                                         ataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTaa                                                           >  2:839158/1‑100 (MQ=255)
                                          taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                          >  2:314013/1‑100 (MQ=255)
                                          taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                          >  1:2235002/1‑100 (MQ=255)
                                          taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                          <  2:2608063/100‑1 (MQ=255)
                                                atTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaa                                                    >  2:1855888/1‑100 (MQ=255)
                                                atTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaa                                                    >  1:3834428/1‑100 (MQ=255)
                                                 ttGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaaa                                                   <  1:35307/100‑1 (MQ=255)
                                                 ttGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaaa                                                   <  2:2610352/100‑1 (MQ=255)
                                                  ttaaTTGTTATTTTTGTGTGTTGTTTATTTTAAAAATCCTTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                  <  1:1095738/98‑1 (MQ=255)
                                                  tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                  <  1:2499689/100‑1 (MQ=255)
                                                   gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACAGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                  <  2:1194103/99‑1 (MQ=255)
                                                   gAATTGTTATTTTTGTGTGTTATTTAAAATAAAAAGACGTATTTATTAATACGGCTTTAATTAATAAAAGGCATGTTGTATTAAAATTAATATTCAAAgt                                                 <  2:2293505/100‑2 (MQ=255)
                                                       tGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATAAAAATTAATATTCAAAGCATaa                                             >  1:1574023/1‑100 (MQ=255)
                                                        gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                            >  1:818872/1‑100 (MQ=255)
                                                         ttATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                           >  1:3739827/1‑100 (MQ=255)
                                                            tttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGa                                        >  2:3552389/1‑100 (MQ=255)
                                                             ttttttttGTTTTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAATCATAAACCGa                                        >  2:263304/1‑99 (MQ=255)
                                                               ttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAg                                     <  2:1308639/100‑1 (MQ=255)
                                                               ttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAg                                     <  1:2475247/100‑1 (MQ=255)
                                                               ttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAg                                     <  2:3452854/100‑1 (MQ=255)
                                                                 gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTACAATTAATATTCAAAGCATAAACCGATaac                                    >  2:3072462/1‑97 (MQ=255)
                                                                 gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcn                                   <  1:35990/100‑2 (MQ=255)
                                                                 gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                                   <  2:1366723/100‑1 (MQ=255)
                                                                 gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                                   <  1:3656468/100‑1 (MQ=255)
                                                                    tgttgtTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGGTattggct                                >  1:321519/1‑94 (MQ=255)
                                                                     gttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATa                               >  2:3159189/1‑100 (MQ=255)
                                                                        gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                            >  2:2088340/1‑100 (MQ=255)
                                                                        gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATCAATATTCAAAGCATAAACCGATAGCCAATACCg                            >  1:1527865/1‑100 (MQ=255)
                                                                          ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                           >  2:981837/1‑99 (MQ=255)
                                                                          ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                          <  1:3073769/100‑1 (MQ=255)
                                                                           tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                         <  2:2005329/100‑1 (MQ=255)
                                                                           tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                         <  2:1573160/100‑1 (MQ=255)
                                                                              gATAAAAAGCGGTATGTATTATTACGGCTTTAAGTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCAGAAACCGATAGCCAATAGCGGtt                         <  2:2457737/97‑1 (MQ=255)
                                                                               aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGTTGTATTAAAATTAATATTCAAAGCAGAAACCGATAGCCAATACCGGTTTCAg                     <  1:108634/100‑1 (MQ=255)
                                                                                tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                     >  2:2542637/1‑99 (MQ=255)
                                                                                tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                    >  2:2116734/1‑100 (MQ=255)
                                                                                attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                    >  2:405226/4‑100 (MQ=37)
                                                                                  aaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                  >  1:3465132/1‑100 (MQ=255)
                                                                                  aaaaGCCGTAATAATAAATAGGGCTTTAATTAATAAAAGGCAGGCTGCATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAACGGTTTAAGTaa                  >  2:2395994/1‑100 (MQ=25)
                                                                                   aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGTCATGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt                 <  1:1058337/100‑1 (MQ=37)
                                                                                   aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt                 >  1:1340922/1‑100 (MQ=37)
                                                                                   aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt                 >  2:1058326/1‑100 (MQ=37)
                                                                                   aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACAGATAGCCAATACCGGTTTCAGTAAt                 >  1:3385369/1‑100 (MQ=37)
                                                                                    aaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg                <  1:514426/100‑1 (MQ=255)
                                                                                    aaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg                <  1:2621422/100‑1 (MQ=255)
                                                                                     ccgatcttttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGa               <  1:382928/93‑1 (MQ=255)
                                                                                        cGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAAt            >  1:2792736/1‑100 (MQ=255)
                                                                                         gTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTTAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATg           >  2:3808546/1‑100 (MQ=255)
                                                                                          tatttattATTACGGCTTTAATTAATAAAAGGCATGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATgg          <  1:3218632/100‑1 (MQ=255)
                                                                                           atttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATACAACGATAGACACTACCGGTTTCAGTAATTaaagagc         >  2:3304362/1‑96 (MQ=38)
                                                                                           atttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCCATAGCAAATACCGGTTTCAGTAATGAAATg           >  1:691703/1‑98 (MQ=255)
                                                                                            tttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaagatcgg        >  2:382922/1‑94 (MQ=255)
                                                                                              tattattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTg      >  1:1100901/1‑100 (MQ=255)
                                                                                               attattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgg     <  1:796959/100‑1 (MQ=255)
                                                                                               attattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgg     >  2:1499740/1‑100 (MQ=255)
                                                                                               attattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgg     <  2:2189784/100‑1 (MQ=255)
                                                                                                 tattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTggg    <  2:223431/99‑1 (MQ=255)
                                                                                                  attaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCTATAGCCAATACCGGTATCAGTAATAAAATGGCGTgggcg  >  1:775733/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCG  >  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: