Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,011 T→A 16.8% 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,0110TA16.8% 491.1 / 85.8 192intergenic (+171/+45)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (72/87);  new base A (16/16);  total (89/103)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

Only 100 of 201 total aligned reads displayed.

CAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGC  >  NC_000913/720912‑721108
                                                                                                   |                                                                                                 
cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                                                                                   <  1:2324404/100‑1 (MQ=255)
   tCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAc                                                                                                <  1:2549277/100‑1 (MQ=255)
       aaTTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                                                                                            >  2:33208/1‑100 (MQ=255)
        atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                           >  1:1264349/1‑100 (MQ=255)
        atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                           >  1:105318/1‑100 (MQ=255)
         tttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                          <  1:1425489/100‑1 (MQ=255)
         attatttaGTGATGAGCTTTTTACTCAGTAATCTAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                          <  2:2906619/99‑1 (MQ=255)
          ttattttAGTATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTGTGGTGTGTTGTTTCAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                          <  1:3187917/99‑1 (MQ=255)
          ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTtta                                                                                         <  1:280027/100‑1 (MQ=255)
          ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTtta                                                                                         <  2:3793325/100‑1 (MQ=255)
            atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                                                                                       >  2:53982/1‑100 (MQ=255)
            atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                                                                                       >  1:139479/1‑100 (MQ=255)
            atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                                                                                       <  2:2813866/100‑1 (MQ=255)
            atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttttta                                                                                       <  1:1434342/100‑4 (MQ=38)
             tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                      <  1:1866943/100‑1 (MQ=255)
             tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                      <  2:237276/100‑1 (MQ=255)
             tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                      <  1:3464436/100‑1 (MQ=255)
             tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                      <  1:2573171/100‑1 (MQ=255)
             tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttat                                                                                      <  1:942061/100‑5 (MQ=37)
              ttaATGATGAGCTTTTTAATCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTATTTAAGATAATAAACCGTATTTATTATTACTGCTttatt                                                                                       >  2:1657609/1‑96 (MQ=255)
               taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                    <  1:1536596/100‑7 (MQ=37)
                aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGACTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                                   <  2:2652187/100‑8 (MQ=25)
                 atgatgATTTTTTTACTAAGTAATATAAAATATTGAATTGTTATTTTTTTATATAAATAAAAAAAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                      <  2:1590967/96‑1 (MQ=255)
                 atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                  <  1:803844/100‑1 (MQ=255)
                 atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                  <  2:3095431/100‑1 (MQ=255)
                 atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                  <  2:2959088/100‑1 (MQ=255)
                    atgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa                                                                               <  2:2533597/100‑1 (MQ=255)
                     tgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag                                                                              >  1:760766/1‑100 (MQ=255)
                      gaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                             >  1:2832869/1‑100 (MQ=255)
                      gaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                             >  2:366907/1‑100 (MQ=255)
                       aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc                                                                            >  2:101541/1‑100 (MQ=255)
                       aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                             <  2:2860103/99‑1 (MQ=255)
                        tcTTTTTACTCAGTAATATAAAATATTGAACTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                           >  2:580152/2‑100 (MQ=255)
                         cTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag                                                                          >  1:377080/1‑100 (MQ=255)
                           ttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc                                                                        <  2:2252609/100‑1 (MQ=255)
                               aCTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa                                                                    <  2:377080/100‑1 (MQ=255)
                                cTCAGTAATATAAAATATTGAATTGTTATTTTTTTTTTTTTTTTAAGATAAAAAGCCGTATTTATTATTACGTCTTTAATTAATAAAAGGCAGGCTGTAt                                                                   >  2:3904168/1‑100 (MQ=255)
                                 tgaGTAATATAAAATATTGAATTGGTATTTTTGTGTGTTGATCAAGATAAAAAGCCGTATTAATTATTACGGTTTTAATTAATAAAAGGCAGGCAGTa                                                                    >  1:3662521/3‑98 (MQ=255)
                                 tCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAtt                                                                  <  2:1784387/100‑1 (MQ=255)
                                        tataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTa                                                           <  2:1713937/100‑1 (MQ=255)
                                         ataAAATATTGAATTATTATTATTGTGTGATGTTTAAGATACAAAGCCGTATTTATTATTACTGTTTTATTTAATAAAAGGAATGCTGTATAAAAAtt                                                            >  2:387289/1‑98 (MQ=255)
                                          taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                         <  2:3232527/100‑1 (MQ=255)
                                           aaaaTATTGAATTGTTATTTTTTTTTGTTTTTTAAGATAAAAAGCCGTATTTATTATTACGTCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                        >  2:2609057/1‑100 (MQ=255)
                                           aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                        <  2:1428931/100‑1 (MQ=255)
                                           aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                        <  1:3676696/100‑1 (MQ=255)
                                           aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                        <  1:1358217/100‑1 (MQ=255)
                                               tatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCa                                                    <  2:606808/100‑1 (MQ=255)
                                                  tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                 <  2:700763/100‑1 (MQ=255)
                                                  tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                 <  2:216942/100‑1 (MQ=255)
                                                   gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                                >  2:85457/1‑100 (MQ=255)
                                                     aTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCAt                                              <  2:3431202/100‑1 (MQ=255)
                                                        gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                           <  1:1925731/100‑1 (MQ=255)
                                                          tATTTTTGTGTGTTGTTTAAGATAAAAAGTCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                         <  1:1971363/100‑1 (MQ=255)
                                                          tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                         >  1:2282011/1‑100 (MQ=255)
                                                          tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                         <  2:413344/100‑1 (MQ=255)
                                                            tttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGa                                       >  2:3594197/1‑100 (MQ=255)
                                                             ttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAt                                      <  2:945083/100‑1 (MQ=255)
                                                             ttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAt                                      >  2:2779996/1‑100 (MQ=255)
                                                             ttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAt                                      <  1:3040752/100‑1 (MQ=255)
                                                              tttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                     <  1:968911/100‑1 (MQ=255)
                                                              tttGTGTGTTCTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCACGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                     <  2:1307009/100‑1 (MQ=255)
                                                               ttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAg                                    <  2:467601/100‑1 (MQ=255)
                                                                  tgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGc                                   >  2:1874606/1‑98 (MQ=255)
                                                                  tgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCa                                 >  2:2460876/1‑100 (MQ=255)
                                                                   gtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGCTATCCaa                                >  1:3169735/1‑100 (MQ=255)
                                                                   gtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCaa                                >  1:436217/1‑100 (MQ=255)
                                                                    tgttgtTTAAGATAAAAAGCCGTATTTATTAGTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAAt                               >  1:2657916/1‑100 (MQ=255)
                                                                     gttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATa                              >  1:547399/1‑100 (MQ=255)
                                                                       tgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAcc                            >  1:2695180/1‑100 (MQ=255)
                                                                         tttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGCCTGTATTAAAATTAATATTCAAAGCATAAAGCGATAGCCAATACCgg                          >  2:1397240/1‑100 (MQ=255)
                                                                          ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAGTTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                         >  1:2777916/1‑100 (MQ=255)
                                                                          ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                         <  1:2117212/100‑1 (MQ=255)
                                                                          ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                         <  2:301589/100‑1 (MQ=255)
                                                                          ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                         >  2:3413853/1‑100 (MQ=255)
                                                                           tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                        <  1:1772599/100‑1 (MQ=255)
                                                                             ttaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                      >  1:3672521/7‑100 (MQ=37)
                                                                             aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                      <  2:1920930/100‑1 (MQ=255)
                                                                             aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                      <  2:3024636/100‑1 (MQ=255)
                                                                             aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                      >  2:1739328/1‑100 (MQ=255)
                                                                              taattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGATGTATTAAAATTAATATTCAAAACATAAACCGATAGCCAATAATGGttata                     >  2:2704326/6‑97 (MQ=21)
                                                                               aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                    <  2:1247968/100‑1 (MQ=255)
                                                                               aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                    <  2:112488/100‑1 (MQ=255)
                                                                               aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                    <  1:2437462/100‑1 (MQ=255)
                                                                                attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                   >  1:2366180/4‑100 (MQ=37)
                                                                                 aaaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTa                  >  1:3298064/1‑100 (MQ=255)
                                                                                   aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                 >  2:1135821/1‑99 (MQ=37)
                                                                                   aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                 >  2:3763690/1‑99 (MQ=37)
                                                                                   aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                 >  2:3202605/1‑99 (MQ=37)
                                                                                    aaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg               <  2:2421715/100‑1 (MQ=255)
                                                                                    aaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                 >  2:317807/1‑98 (MQ=37)
                                                                                       ccGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaaa            >  1:1355385/1‑100 (MQ=255)
                                                                                        cGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAAt           <  1:3147234/100‑1 (MQ=255)
                                                                                         gTATTAATTATTACGGCTTTAATCAATAAAAGGCAGGCGGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAAt           <  2:1127862/99‑1 (MQ=255)
                                                                                           atttattATTACGGCTTTAATTAATAAAAGGCAGTCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCTAATACCGGTTTCAGTAATGAAATGGc        >  2:516018/1‑100 (MQ=255)
                                                                                           atttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGc        >  1:1564469/1‑100 (MQ=255)
                                                                                            tttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCg       >  1:1145353/1‑100 (MQ=255)
                                                                                             ttattattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGt      <  2:2835799/100‑1 (MQ=255)
                                                                                               attattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgg    >  1:601644/1‑100 (MQ=255)
                                                                                                ttattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTggg   <  2:3722737/100‑1 (MQ=255)
                                                                                                 tattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgggc  >  1:2129639/1‑100 (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: