Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,004 T→A 15.5% intergenic (+164/+52) ybfK → / ← kdpE uncharacterized protein/response regulator in two‑component regulatory system with KdpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913721,0040TA15.5% 371.8 / 54.2 141intergenic (+164/+52)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (45/74);  new base A (10/12);  total (55/86)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 150 total aligned reads displayed.

ATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCG  >  NC_000913/720907‑721103
                                                                                                 |                                                                                                   
aTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattt                                                                                                   <  2:2255296/100‑1 (MQ=255)
 tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattta                                                                                                  <  2:1252555/100‑1 (MQ=255)
  gATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttat                                                                                                 <  2:2232970/100‑1 (MQ=255)
   aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                                <  2:1602881/100‑1 (MQ=255)
   aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                                >  2:3200469/1‑100 (MQ=255)
     cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                                                                              <  2:150535/100‑1 (MQ=255)
     cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                                                                              >  1:3688690/1‑100 (MQ=255)
         cTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                                                                                          <  1:3815168/100‑1 (MQ=255)
          ttAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACgg                                                                                         <  2:1123376/100‑1 (MQ=255)
           tAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                        <  2:1736581/100‑1 (MQ=255)
             atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtt                                                                                      <  2:3085011/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)
                 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)
                  tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttat                                                                                 <  2:540267/100‑5 (MQ=37)
                   ttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatta                                                                                <  2:1185407/100‑1 (MQ=255)
                    taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                                                                               >  1:3155778/1‑100 (MQ=255)
                    taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                               <  2:3590021/100‑7 (MQ=37)
                    taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                               <  2:2771065/100‑7 (MQ=37)
                    taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                               <  2:2756019/100‑7 (MQ=37)
                    taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                               <  2:2635973/100‑7 (MQ=37)
                     gatcttgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                              <  2:1297867/95‑1 (MQ=255)
                     attgatgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                              <  2:90108/98‑1 (MQ=255)
                     aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                              <  2:3342864/100‑1 (MQ=255)
                     aaTGAGGAGCTTTTTACTCAGTAATATAAAAAATGGAATTGTTATTTTTGTGTGTGGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                              <  1:3245187/100‑8 (MQ=21)
                      atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                             <  1:1728791/100‑1 (MQ=255)
                      atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                             <  2:3737110/100‑1 (MQ=255)
                      atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                             >  1:2546662/1‑100 (MQ=255)
                      atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                             <  2:3197637/100‑1 (MQ=255)
                       tgatgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTTTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaataa                                                                            <  1:1196757/100‑1 (MQ=255)
                         atgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa                                                                          <  2:1500898/100‑1 (MQ=255)
                          tgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTCTTATTACGGCTTTAATTAATaaaa                                                                          >  1:3754029/1‑99 (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)
                                     cTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa                                                               <  2:1119096/99‑1 (MQ=255)
                                     cTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAt                                                              >  1:2232432/1‑100 (MQ=255)
                                     cTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAt                                                              <  1:738800/100‑1 (MQ=255)
                                      tCAGTAATATAAAATATTGAATTTTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAtt                                                             <  1:3239101/100‑1 (MQ=255)
                                      tCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAtt                                                             >  2:1245996/1‑100 (MQ=255)
                                       cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                            >  2:2446576/1‑100 (MQ=255)
                                       cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                            >  2:1309675/1‑100 (MQ=255)
                                          tAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaaa                                                         <  1:979386/100‑1 (MQ=255)
                                            atatAAAATATTGAATTGTTATTTTTTTGTTTTGTTTAAAATAAAAAGCTGTATTTATTATTACGCCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                        <  1:1563410/99‑1 (MQ=255)
                                            atatAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAtt                                                       <  2:1026794/100‑1 (MQ=255)
                                             gccaacatatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAACAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTa                                                      <  1:2447824/94‑1 (MQ=255)
                                              cttccgatctGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTaa                                                     <  2:2054091/91‑1 (MQ=255)
                                               taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                    <  2:2626091/100‑1 (MQ=255)
                                                aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGTCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                   <  1:296554/100‑1 (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: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)
                                                             gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                      >  1:899938/1‑100 (MQ=255)
                                                                aTTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCg                                   <  2:436590/100‑1 (MQ=255)
                                                                 ttttttTGTGTTGTTTAAAATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                    <  1:2366110/98‑1 (MQ=255)
                                                                 tttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGa                                  >  2:2478376/1‑100 (MQ=255)
                                                                 tttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGa                                  <  2:3804109/100‑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                             <  2:398689/100‑1 (MQ=255)
                                                                      gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                             <  1:1309667/100‑1 (MQ=255)
                                                                      gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                             <  1:85930/100‑1 (MQ=255)
                                                                          gttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATa                         >  2:2828960/1‑100 (MQ=255)
                                                                           ttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAc                        <  2:2221638/100‑1 (MQ=255)
                                                                           ttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAc                        >  1:1547291/1‑100 (MQ=255)
                                                                             gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                      >  1:3778573/1‑100 (MQ=255)
                                                                              tttattaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                     >  2:703858/11‑100 (MQ=25)
                                                                                tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                   >  1:3298580/1‑100 (MQ=255)
                                                                                 aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                  >  2:2754069/1‑100 (MQ=255)
                                                                                  aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                 >  2:1599119/1‑100 (MQ=255)
                                                                                   gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCTAATATCGGTTTCa                >  1:3471942/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)
                                                                                     attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt              >  1:7885/4‑100 (MQ=37)
                                                                                       aaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa            <  1:2861600/100‑1 (MQ=255)
                                                                                        aaaGCCGTATTTATTATTACGGGTTTAACTACTAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt           <  1:2568015/100‑1 (MQ=255)
                                                                                        aaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt           <  2:197224/100‑1 (MQ=255)
                                                                                        aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt           >  2:2317528/1‑100 (MQ=37)
                                                                                         aaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGACTGTATAAAAATTAATATTCAAAACATAAAACCATAACCAATAACGGTTTCAGCAATg          >  1:1315174/1‑100 (MQ=255)
                                                                                           gCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaa        <  2:2258819/100‑1 (MQ=255)
                                                                                            ccGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaaa       >  1:1384691/1‑100 (MQ=255)
                                                                                             cGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAAt      <  1:3584641/100‑1 (MQ=255)
                                                                                              gTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATg     <  1:259925/100‑1 (MQ=255)
                                                                                                 tttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTAAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCg  >  1:437336/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
ATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCG  >  NC_000913/720907‑721103

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: