Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,006 T→A 18.0% intergenic (+166/+50) 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,0060TA18.0% 471.5 / 87.7 188intergenic (+166/+50)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (69/85);  new base A (17/17);  total (86/102)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 202 total aligned reads displayed.

GATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTG  >  NC_000913/720909‑721105
                                                                                                 |                                                                                                   
  tCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTACGATAAACAGCCGTCTtcatta                                                                                                 <  2:3700192/100‑1 (MQ=255)
  tCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttatta                                                                                                 >  2:1450573/1‑100 (MQ=255)
   cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                                                                                <  1:2324404/100‑1 (MQ=255)
       cTTAATTTATTTAATGATGAGCTTTTTACTCAGTACTATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatta                                                                                              >  1:1586773/1‑98 (MQ=255)
          aaTTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                                                                                         >  2:33208/1‑100 (MQ=255)
           atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                        >  1:105318/1‑100 (MQ=255)
           atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                        >  1:431378/1‑100 (MQ=255)
           atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                        >  1:1264349/1‑100 (MQ=255)
            tttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                       <  1:1425489/100‑1 (MQ=255)
            attatttaGTGATGAGCTTTTTACTCAGTAATCTAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                       <  2:2906619/99‑1 (MQ=255)
             ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTtta                                                                                      <  2:3793325/100‑1 (MQ=255)
             ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttn                                                                                      <  1:3179203/100‑2 (MQ=38)
               atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                                                                                    >  1:139479/1‑100 (MQ=255)
               atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                                                                                    <  2:2813866/100‑1 (MQ=255)
               atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                                                                                    >  2:53982/1‑100 (MQ=255)
               atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttttta                                                                                    <  1:1434342/100‑4 (MQ=38)
                tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                   <  2:237276/100‑1 (MQ=255)
                tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                   <  1:2583076/100‑1 (MQ=255)
                tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttat                                                                                   <  1:942061/100‑5 (MQ=37)
                 ttaATGATGAGCTTTTTAATCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTATTTAAGATAATAAACCGTATTTATTATTACTGCTttatt                                                                                    >  2:1657609/1‑96 (MQ=255)
                  taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                 <  1:1829335/100‑7 (MQ=37)
                  taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                 <  1:3562770/100‑7 (MQ=37)
                   aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGACTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                                <  2:2652187/100‑8 (MQ=25)
                   aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                                <  1:1920684/100‑8 (MQ=25)
                    atgatgATTTTTTTACTAAGTAATATAAAATATTGAATTGTTATTTTTTTATATAAATAAAAAAAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                   <  2:1590967/96‑1 (MQ=255)
                    atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                               <  2:2959088/100‑1 (MQ=255)
                    atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                               <  1:803844/100‑1 (MQ=255)
                    atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                               <  2:3095431/100‑1 (MQ=255)
                       atgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa                                                                            <  2:2533597/100‑1 (MQ=255)
                        tgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag                                                                           <  1:2118425/100‑1 (MQ=255)
                        tgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag                                                                           >  1:760766/1‑100 (MQ=255)
                         gaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                          >  1:525207/1‑100 (MQ=255)
                         gaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                          >  2:366907/1‑100 (MQ=255)
                         gaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                          >  1:2456564/1‑100 (MQ=255)
                          aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc                                                                         >  2:101541/1‑100 (MQ=255)
                          aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc                                                                         <  1:3673347/100‑1 (MQ=255)
                          aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                          <  2:2860103/99‑1 (MQ=255)
                           tcTTTTTACTCAGTAATATAAAATATTGAACTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                        >  2:580152/2‑100 (MQ=255)
                             tttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcagg                                                                      >  1:374854/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)
                                   cTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAt                                                                >  1:2097050/1‑100 (MQ=255)
                                    tCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAtt                                                               <  2:1784387/100‑1 (MQ=255)
                                          atatAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAtt                                                         >  1:1493322/1‑100 (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:1358217/100‑1 (MQ=255)
                                               aaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAatat                                                    >  1:1954847/1‑100 (MQ=255)
                                                 atatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTc                                                  <  1:2782246/100‑1 (MQ=255)
                                                  tatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCa                                                 <  2:606808/100‑1 (MQ=255)
                                                     tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaaa                                               <  1:2555378/99‑1 (MQ=255)
                                                     tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                              <  2:700763/100‑1 (MQ=255)
                                                     tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                              <  2:216942/100‑1 (MQ=255)
                                                      gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                             >  2:85457/1‑100 (MQ=255)
                                                      gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                             <  2:1363929/100‑1 (MQ=255)
                                                        aTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCAt                                           <  2:3431202/100‑1 (MQ=255)
                                                             tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                      <  2:413344/100‑1 (MQ=255)
                                                             tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                      <  1:320246/100‑1 (MQ=255)
                                                             tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                      >  1:2282011/1‑100 (MQ=255)
                                                               tttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGa                                    >  2:3594197/1‑100 (MQ=255)
                                                                ttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAt                                   >  2:2779996/1‑100 (MQ=255)
                                                                ttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAt                                   <  2:945083/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)
                                                                    gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                               <  1:2827102/100‑1 (MQ=255)
                                                                     tgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGc                                >  2:1874606/1‑98 (MQ=255)
                                                                     tgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCa                              >  2:2460876/1‑100 (MQ=255)
                                                                      gtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCaa                             >  1:3169489/1‑100 (MQ=255)
                                                                            tttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAcctg                       <  1:2675553/100‑3 (MQ=255)
                                                                             ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                      <  2:301589/100‑1 (MQ=255)
                                                                             ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                      >  2:3413853/1‑100 (MQ=255)
                                                                                ttaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                   >  1:2949306/7‑100 (MQ=37)
                                                                                ttaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                   >  1:3668675/7‑100 (MQ=37)
                                                                                aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                   <  2:1920930/100‑1 (MQ=255)
                                                                                aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                   >  2:1739328/1‑100 (MQ=255)
                                                                                aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                   <  2:3024636/100‑1 (MQ=255)
                                                                                 taattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGATGTATTAAAATTAATATTCAAAACATAAACCGATAGCCAATAATGGttata                  >  2:2704326/6‑97 (MQ=21)
                                                                                 gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                  <  1:3107508/100‑1 (MQ=255)
                                                                                  aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                 <  2:112488/100‑1 (MQ=255)
                                                                                  aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                 <  2:1247968/100‑1 (MQ=255)
                                                                                   attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                >  1:3414804/4‑100 (MQ=37)
                                                                                   attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                >  1:2366180/4‑100 (MQ=37)
                                                                                      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)
                                                                                         gCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaa          >  1:585209/1‑100 (MQ=255)
                                                                                          ccGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaaa         >  1:1355385/1‑100 (MQ=255)
                                                                                            gTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATg       >  1:1537867/1‑100 (MQ=255)
                                                                                            gTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATg       <  1:2757584/100‑1 (MQ=255)
                                                                                            gTATTAATTATTACGGCTTTAATCAATAAAAGGCAGGCGGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAAt        <  2:1127862/99‑1 (MQ=255)
                                                                                              atttattATTACGGCTTTAATTAATAAAAGGCAGTCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCTAATACCGGTTTCAGTAATGAAATGGc     >  2:516018/1‑100 (MQ=255)
                                                                                               tttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCg    >  1:1145353/1‑100 (MQ=255)
                                                                                                ttattattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGt   <  2:2835799/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
GATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTG  >  NC_000913/720909‑721105

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: