Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,009 T→A 15.9% intergenic (+169/+47) 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,0090TA15.9% 367.7 / 54.1 138intergenic (+169/+47)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (44/72);  new base A (10/12);  total (54/84)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 151 total aligned reads displayed.

ATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGC  >  NC_000913/720910‑721108
                                                                                                   |                                                                                                   
aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                                     <  2:3310055/100‑1 (MQ=255)
aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                                     <  2:1602881/100‑1 (MQ=255)
  cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                                                                                   <  2:150535/100‑1 (MQ=255)
  cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                                                                                   >  1:3688690/1‑100 (MQ=255)
  cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                                                                                   <  2:2801879/100‑1 (MQ=255)
      cTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                                                                                               <  2:1151129/100‑1 (MQ=255)
       ttAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACgg                                                                                              <  2:1123376/100‑1 (MQ=255)
        tAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                             <  2:1736581/100‑1 (MQ=255)
          atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                           <  2:2718972/100‑1 (MQ=255)
           tttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttt                                                                                          <  1:1050899/100‑1 (MQ=38)
            ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                         >  1:2635960/1‑99 (MQ=38)
             tatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaa                                                                                        <  2:1679223/100‑1 (MQ=255)
             tatttaATGATGAGCTTTTTACTCAATATTATAAAATATTGAATTGTTATTTTTGTGTGTTGTTAAGATAAAAAGCCGCAATAATAAATACGGCttttta                                                                                       <  1:3573250/100‑4 (MQ=25)
              atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTAGTTTAAGATATAAAGCCGTATTTATTATTATGGCTtta                                                                                         >  1:2980570/1‑98 (MQ=255)
               tttaATGATGATCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttat                                                                                      <  2:2187086/100‑5 (MQ=37)
               tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                      <  2:3321899/100‑1 (MQ=255)
               tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttat                                                                                      <  1:1756991/100‑5 (MQ=38)
                 taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                    <  2:2771065/100‑7 (MQ=37)
                  gatcttgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                                   <  2:1297867/95‑1 (MQ=255)
                   atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                  <  2:3197637/100‑1 (MQ=255)
                   atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                  >  1:3697655/1‑100 (MQ=255)
                   atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                  <  1:1728791/100‑1 (MQ=255)
                   atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                  <  2:2546674/100‑1 (MQ=255)
                   atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                  <  2:3737110/100‑1 (MQ=255)
                    tgatgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTTTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaataa                                                                                 <  1:1196757/100‑1 (MQ=255)
                      atgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa                                                                               >  1:2581732/1‑100 (MQ=255)
                      atgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa                                                                               <  2:1500898/100‑1 (MQ=255)
                          gCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                           <  2:2460327/100‑1 (MQ=255)
                           cTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag                                                                          >  2:290545/1‑100 (MQ=255)
                            tttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTCTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                           >  1:2112695/1‑98 (MQ=255)
                             ttttACTAAGTAATATAAAATATTAAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc                                                                        <  1:2255389/100‑1 (MQ=255)
                               ttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTg                                                                      <  1:3152110/100‑1 (MQ=255)
                                  cTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAt                                                                   <  1:738800/100‑1 (MQ=255)
                                  cTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAt                                                                   >  1:2232432/1‑100 (MQ=255)
                                   tCAGTAATATAAAATATTGAATTTTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAtt                                                                  <  1:3239101/100‑1 (MQ=255)
                                   tCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAtt                                                                  >  2:1245996/1‑100 (MQ=255)
                                    cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                                 >  2:1309675/1‑100 (MQ=255)
                                    cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                                 >  2:2446920/1‑100 (MQ=255)
                                       tAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaaa                                                              <  1:979386/100‑1 (MQ=255)
                                       tAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaaa                                                              <  2:1080694/100‑1 (MQ=255)
                                       tAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAAATAAAAAGCCGTATTTCTTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaaa                                                              <  1:3788131/100‑1 (MQ=255)
                                         atatAAAATATTGAATTGTTATTTTTTTGTTTTGTTTAAAATAAAAAGCTGTATTTATTATTACGCCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                             <  1:1563410/99‑1 (MQ=255)
                                         atatAAAATATTGAATTGTTATTTTTGTGTGTTGCTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAATGAAATAATGTATTaaaact                                                            >  1:3171001/1‑98 (MQ=255)
                                          gccaacatatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAACAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTa                                                           <  1:2447824/94‑1 (MQ=255)
                                           cttccgatctGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTaa                                                          <  2:2054091/91‑1 (MQ=255)
                                            taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                         <  2:2626091/100‑1 (MQ=255)
                                            taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                         <  1:2942348/100‑1 (MQ=255)
                                             aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                        >  1:3246181/1‑100 (MQ=255)
                                             aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                        >  2:2054493/1‑100 (MQ=255)
                                              aaaTATTAAATTTTTATTTTTGTGTTGTGTTTAAAATAAAAAGCCGTATTTAATATAACGGATTTAATTAATAAAAGGCAGTCTGTATTAAAATTAatat                                                       <  1:1713829/100‑1 (MQ=255)
                                               aaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATAtt                                                      >  1:1070043/1‑100 (MQ=255)
                                                 tatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCa                                                    <  1:2100905/100‑1 (MQ=255)
                                                    tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                 <  2:2835220/100‑1 (MQ=255)
                                                    tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                 <  2:2062985/100‑1 (MQ=255)
                                                     gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTaaagatcggaag                                                >  1:2054081/1‑90 (MQ=255)
                                                      aatttctATTTTTGTTTATTGTTAAAGATAAAAAGCCGAATTTACTATTACGACTTTAATTAATAAAAGGCAGTCTGTATTAAAATTAATATTCAAAg                                                 <  1:2939660/92‑1 (MQ=255)
                                                      aaTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                 <  2:3429021/98‑1 (MQ=255)
                                                      aaTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                                <  1:2794293/99‑1 (MQ=255)
                                                      aaTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCa                                               <  2:765839/100‑1 (MQ=255)
                                                          gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                           >  1:899938/1‑100 (MQ=255)
                                                              ttttttTGTGTTGTTTAAAATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                         <  1:2366110/98‑1 (MQ=255)
                                                              tttttGAGTGTTGGTTAAGATAAAAAGGCGTATTTATTATTACGGCTTGAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAACGa                                       <  1:1850343/100‑1 (MQ=255)
                                                               ttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAt                                      <  2:1238362/100‑1 (MQ=255)
                                                               ttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAt                                      <  1:2416018/100‑1 (MQ=255)
                                                                tttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAtn                                     <  2:1898860/100‑2 (MQ=255)
                                                                   gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                                  <  1:1309667/100‑1 (MQ=255)
                                                                   gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                                  <  2:3814963/100‑1 (MQ=255)
                                                                   gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                                  <  2:398689/100‑1 (MQ=255)
                                                                   gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                                  <  1:85930/100‑1 (MQ=255)
                                                                        ttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAc                             >  1:1547291/1‑100 (MQ=255)
                                                                        ttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAc                             <  2:2221638/100‑1 (MQ=255)
                                                                          gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                           >  1:3778573/1‑100 (MQ=255)
                                                                           tttattaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                          >  2:703858/11‑100 (MQ=25)
                                                                            ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                         <  2:2237989/100‑1 (MQ=255)
                                                                             tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                        >  1:3298580/1‑100 (MQ=255)
                                                                              aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                       >  2:512413/1‑100 (MQ=255)
                                                                              aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                       >  2:2754069/1‑100 (MQ=255)
                                                                               aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                      >  2:1599119/1‑100 (MQ=255)
                                                                                gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                     >  2:2303013/1‑100 (MQ=255)
                                                                                gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                     <  1:1245989/100‑1 (MQ=255)
                                                                                gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                     >  2:1371693/1‑100 (MQ=255)
                                                                                 aattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                     >  1:1472886/5‑99 (MQ=37)
                                                                                 aTAAAAAGCCGTATTTATTATTCCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTTAAAGCATAAACCGATAGCCAATACCGGTTTCAg                    <  2:921042/100‑1 (MQ=255)
                                                                                  tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                   <  2:1426109/100‑1 (MQ=255)
                                                                                  tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                   >  2:3385213/1‑100 (MQ=255)
                                                                                  attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                   >  2:3101077/4‑100 (MQ=37)
                                                                                  attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                   >  1:7885/4‑100 (MQ=37)
                                                                                     aaaGCCGTATTTATTATTACGGGTTTAACTACTAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt                <  1:2568015/100‑1 (MQ=255)
                                                                                     aaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt                <  2:197224/100‑1 (MQ=255)
                                                                                     aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt                >  2:2317528/1‑100 (MQ=37)
                                                                                      aaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg               <  2:3733476/100‑1 (MQ=255)
                                                                                      aaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGACTGTATAAAAATTAATATTCAAAACATAAAACCATAACCAATAACGGTTTCAGCAATg               >  1:1315174/1‑100 (MQ=255)
                                                                                         ccGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaaa            >  1:1384691/1‑100 (MQ=255)
                                                                                              tttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTAAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCg       >  1:437336/1‑100 (MQ=255)
                                                                                               ttattattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATGCCGGTTTCAGTAACGAAATAGc        >  1:3114054/1‑98 (MQ=255)
                                                                                                tattattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTg     <  2:529024/100‑1 (MQ=255)
                                                                                                tattattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTg     >  2:1011806/1‑100 (MQ=255)
                                                                                                 attattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgg    <  1:3486777/100‑1 (MQ=255)
                                                                                                   tattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgggc  >  1:3010686/1‑100 (MQ=255)
                                                                                                   tattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgggc  >  1:454490/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGC  >  NC_000913/720910‑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: