Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,009 T→A 17.0% 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,0090TA17.0% 483.5 / 83.5 188intergenic (+169/+47)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (70/86);  new base A (16/16);  total (86/102)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 204 total aligned reads displayed.

TCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGC  >  NC_000913/720911‑721108
                                                                                                  |                                                                                                   
tCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttatta                                                                                                    >  2:1450573/1‑100 (MQ=255)
 cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                                                                                   <  1:2324404/100‑1 (MQ=255)
    tCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAc                                                                                                <  1:2549277/100‑1 (MQ=255)
    tCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAc                                                                                                <  1:951865/100‑1 (MQ=255)
       tAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGATTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                             <  2:1027134/100‑1 (MQ=255)
         atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                           >  1:105318/1‑100 (MQ=255)
         atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                           >  1:1264349/1‑100 (MQ=255)
         atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                           >  1:431378/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)
           ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTtta                                                                                         <  1:840164/100‑1 (MQ=255)
           ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttn                                                                                         <  1:3179203/100‑2 (MQ=38)
             atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                                                                                       >  1:991376/1‑100 (MQ=255)
             atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                                                                                       >  2:53982/1‑100 (MQ=255)
             atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                                                                                       >  2:1125817/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)
             atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttttta                                                                                       <  1:605804/100‑4 (MQ=38)
              tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                      <  1:1866943/100‑1 (MQ=255)
              tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                      <  1:2573171/100‑1 (MQ=255)
              tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                      <  2:237276/100‑1 (MQ=255)
              tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttat                                                                                      <  2:3369058/100‑5 (MQ=37)
               ttaATGATGAGCTTTTTAATCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTATTTAAGATAATAAACCGTATTTATTATTACTGCTttatt                                                                                       >  2:1657609/1‑96 (MQ=255)
                taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                    <  1:3562770/100‑7 (MQ=37)
                taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                    <  1:1536596/100‑7 (MQ=37)
                 aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGACTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                                   <  2:2652187/100‑8 (MQ=25)
                 aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                                   <  2:3077040/100‑8 (MQ=25)
                  atgatgATTTTTTTACTAAGTAATATAAAATATTGAATTGTTATTTTTTTATATAAATAAAAAAAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                      <  2:1590967/96‑1 (MQ=255)
                  atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                  <  2:2959088/100‑1 (MQ=255)
                     atgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa                                                                               <  2:2533597/100‑1 (MQ=255)
                     atgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa                                                                               <  2:1299688/100‑1 (MQ=255)
                      tgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag                                                                              <  2:1129765/100‑1 (MQ=255)
                       gaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                             >  1:525207/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)
                                 cTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGCTGTTTAAGCTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAACGCAGGCTCTAt                                                                   >  2:3600170/1‑100 (MQ=255)
                                  tCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAtt                                                                  <  2:1784387/100‑1 (MQ=255)
                                   cAGTAATATAAAATATTGAATTGTTATTTTTGTTAGTTGTTAACGATAAAAAGGCGTATTTATTGTTACTGCTTTAATTAATAAAACGAAGGCTGTACTa                                                                 >  2:3420973/1‑100 (MQ=255)
                                         tataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTa                                                           <  2:1713937/100‑1 (MQ=255)
                                          ataAAATATTGAATTATTATTATTGTGTGATGTTTAAGATACAAAGCCGTATTTATTATTACTGTTTTATTTAATAAAAGGAATGCTGTATAAAAAtt                                                            >  2:387289/1‑98 (MQ=255)
                                           taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                         >  2:3676668/1‑100 (MQ=255)
                                            aaaaTATTGAATTGTTATTTTTTTTTGTTTTTTAAGATAAAAAGCCGTATTTATTATTACGTCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                        >  2:2609057/1‑100 (MQ=255)
                                            aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                        <  1:1358217/100‑1 (MQ=255)
                                            aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                        <  2:1428931/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:1363929/100‑1 (MQ=255)
                                                    gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                                >  2:85457/1‑100 (MQ=255)
                                                        tGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaa                                            >  2:34671/1‑100 (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)
                                                           tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                         <  1:320246/100‑1 (MQ=255)
                                                              ttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAt                                      >  2:2779996/1‑100 (MQ=255)
                                                              ttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAt                                      <  2:945083/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)
                                                                     tgttgtTTAAGATAAAAAGCCGTATTTATTAGTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAAt                               >  1:2657916/1‑100 (MQ=255)
                                                                        tgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAcc                            >  1:2695180/1‑100 (MQ=255)
                                                                          tttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                          <  2:1211948/100‑1 (MQ=255)
                                                                           ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAGTTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                         >  1:2777916/1‑100 (MQ=255)
                                                                           ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                         <  2:301589/100‑1 (MQ=255)
                                                                           ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                         <  1:2117212/100‑1 (MQ=255)
                                                                           ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                         >  2:3195656/1‑100 (MQ=255)
                                                                            tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                        <  1:1772599/100‑1 (MQ=255)
                                                                              ttaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                      >  1:3668675/7‑100 (MQ=37)
                                                                              ttaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                      >  1:2949306/7‑100 (MQ=37)
                                                                              aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                      >  2:1739328/1‑100 (MQ=255)
                                                                              aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                      <  2:1920930/100‑1 (MQ=255)
                                                                               taattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGATGTATTAAAATTAATATTCAAAACATAAACCGATAGCCAATAATGGttata                     >  2:2704326/6‑97 (MQ=21)
                                                                               gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                     <  1:3107508/100‑1 (MQ=255)
                                                                                aattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                    >  2:3241176/5‑100 (MQ=37)
                                                                                aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                    <  1:2437462/100‑1 (MQ=255)
                                                                                 tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTCAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                   >  2:3217670/1‑100 (MQ=255)
                                                                                 attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                   >  1:3414804/4‑100 (MQ=37)
                                                                                 attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                   >  1:2366180/4‑100 (MQ=37)
                                                                                    aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                 >  2:3763690/1‑99 (MQ=37)
                                                                                     aaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg               <  2:2421715/100‑1 (MQ=255)
                                                                                       gCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaa             >  1:585209/1‑100 (MQ=255)
                                                                                        ccGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaaa            >  1:1355385/1‑100 (MQ=255)
                                                                                            atttattATTACGGCTTTAATTAATAAAAGGCAGTCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCTAATACCGGTTTCAGTAATGAAATGGc        >  2:516018/1‑100 (MQ=255)
                                                                                            atttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGc        >  1:1564469/1‑100 (MQ=255)
                                                                                             tttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACTGATAGCCAATACCGGTTTCAGTAATGAAATGGCg       <  1:794511/100‑1 (MQ=255)
                                                                                             tttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCg       >  1:1145353/1‑100 (MQ=255)
                                                                                              ttattattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGt      <  2:2835799/100‑1 (MQ=255)
                                                                                               tattattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTg     >  2:3097586/1‑100 (MQ=255)
                                                                                                 ttattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTggg   <  2:3722737/100‑1 (MQ=255)
                                                                                                  tattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgggc  >  1:2129639/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGC  >  NC_000913/720911‑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: