Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,009 T→A 16.3% 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,0090TA16.3% 471.7 / 71.7 177intergenic (+169/+47)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (54/94);  new base A (16/13);  total (70/107)
Fisher's exact test for biased strand distribution p-value = 6.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

ATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGC  >  NC_000913/720910‑721108
                                                                                                   |                                                                                                   
aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                                     <  2:3744298/100‑1 (MQ=255)
  cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTCTTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                                                                                   >  1:1001212/1‑100 (MQ=255)
  cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                                     >  2:1271026/1‑98 (MQ=255)
   aaTCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatt                                                                                                  >  2:3225513/1‑100 (MQ=255)
      cTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGTTAAAAAGCCGTATTTATTATTACg                                                                                               <  2:108440/100‑1 (MQ=255)
        tAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                             <  2:1883737/100‑1 (MQ=255)
        tAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                             <  2:917596/100‑1 (MQ=255)
          atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGCATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                           >  1:2981048/1‑100 (MQ=255)
          atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                           <  2:2454507/100‑1 (MQ=255)
           tttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                          <  2:1205571/100‑1 (MQ=255)
           tttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttt                                                                                          <  1:1290896/100‑1 (MQ=38)
               tttaATGATGAGCTTTTTACTCAGTAATATAAAATATAGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttat                                                                                      <  2:3342550/100‑5 (MQ=37)
                ttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTAGTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatc                                                                                     <  2:3357067/100‑6 (MQ=37)
                 taATGATGAGCTTTTTATTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                    <  2:15337/100‑7 (MQ=25)
                 taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                    <  2:3268696/100‑7 (MQ=37)
                 taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGAATAAAATAAAAAGCTGTATTTATTATTACGGCTttaatt                                                                                      <  2:295079/98‑1 (MQ=255)
                  aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                                   <  1:3597887/100‑8 (MQ=37)
                  aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                    <  1:2244991/99‑7 (MQ=37)
                   atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                                  <  2:1811948/100‑1 (MQ=255)
                   atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATGTTTGTGGGTTGTTTAAGATAAAAAGCCGTATTTATTATTACCGCTTTAATtaata                                                                                  >  2:1004054/1‑100 (MQ=255)
                     gatgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaa                                                                                >  1:328875/1‑100 (MQ=255)
                        gaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                             >  2:2640022/1‑100 (MQ=255)
                         aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc                                                                            <  1:1411929/100‑1 (MQ=255)
                             ttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc                                                                        <  1:237569/100‑1 (MQ=255)
                             ttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc                                                                        <  1:606352/100‑1 (MQ=255)
                                tACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt                                                                     >  2:3521534/1‑100 (MQ=255)
                                tACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt                                                                     <  2:11318/100‑1 (MQ=255)
                                 aCTCAGTAATATAAAATATTGAATTGTTATTTGTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa                                                                    >  1:4054328/1‑100 (MQ=255)
                                   tCAGTAATATAAAATATTGAATTTTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAtt                                                                  <  2:1744573/100‑1 (MQ=255)
                                    cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                                 >  1:2933730/1‑100 (MQ=255)
                                    cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                                 >  1:577910/1‑100 (MQ=255)
                                      gTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaa                                                                >  2:974802/1‑99 (MQ=255)
                                          tataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAttn                                                           <  1:1867405/100‑2 (MQ=255)
                                           ataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTTAAATTa                                                           >  2:777105/1‑99 (MQ=255)
                                             aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                        >  1:1791345/1‑100 (MQ=255)
                                                atatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTc                                                     >  1:1063308/1‑100 (MQ=255)
                                                 tatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCa                                                    <  1:1574386/100‑1 (MQ=255)
                                                    tGAATTGTTATTTTTGTTTGTTGTTTAAGATAAATAGCCGTTTTTATTATTACGGCTTTAATTAATAAAAGGCAGTCTGTATTAAAATTAATATTCAAAg                                                 <  2:221875/100‑1 (MQ=255)
                                                    tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                 <  1:2051533/100‑1 (MQ=255)
                                                    tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                 <  1:1169228/100‑1 (MQ=255)
                                                    tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAAGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaaa                                                  <  1:2742936/99‑1 (MQ=255)
                                                     gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                                <  2:2703289/100‑1 (MQ=255)
                                                      aaTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCa                                               <  2:3175012/100‑1 (MQ=255)
                                                        ttGTTATTTTTGTGTGTTGTTTAATATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATa                                             <  2:572783/100‑1 (MQ=255)
                                                        ttGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCAt                                              <  2:1332561/99‑1 (MQ=255)
                                                        ttGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATa                                             <  1:97850/100‑1 (MQ=255)
                                                        ttGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATa                                             <  1:2096559/100‑1 (MQ=255)
                                                            tATTTTTTTTTTTTTTTTAAGATAAAAAGCCGTATTTATTATTATGTTTTTAATTAATAAAAGGCAGGCTGTTTTTAAATTAATAtttatttttt                                              >  2:1428815/1‑87 (MQ=255)
                                                            tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                         >  2:4107172/1‑100 (MQ=255)
                                                            tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                         <  2:1181294/100‑1 (MQ=255)
                                                             aTTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCg                                        >  1:3870495/1‑100 (MQ=255)
                                                               ttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAt                                      <  2:4159719/100‑1 (MQ=255)
                                                                tttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAt                                      >  1:433070/1‑99 (MQ=255)
                                                                tttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                     <  2:3752914/100‑1 (MQ=255)
                                                                tttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                     <  2:726559/100‑1 (MQ=255)
                                                                tttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                     <  1:3766515/100‑1 (MQ=255)
                                                                 ttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAg                                    <  1:2181777/100‑1 (MQ=255)
                                                                   gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                                  <  2:2200371/100‑1 (MQ=255)
                                                                    tgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCa                                 <  1:2646232/100‑1 (MQ=255)
                                                                    tgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCa                                 >  1:1234035/1‑100 (MQ=255)
                                                                    tgagttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCa                                 >  2:1778088/4‑100 (MQ=255)
                                                                     gtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCaa                                >  2:2867365/1‑100 (MQ=255)
                                                                     gtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCa                                 <  2:470797/99‑1 (MQ=255)
                                                                      tgttgtTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAAt                               >  1:203069/1‑100 (MQ=255)
                                                                      tgttgtTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAAt                               <  2:2295466/100‑1 (MQ=255)
                                                                       gttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATa                              <  1:694976/100‑1 (MQ=255)
                                                                         tgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAcc                            <  1:1099705/100‑1 (MQ=255)
                                                                          gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATAATCCAAGCATAAACCGAGAGCCAATAc                             >  2:3882185/1‑98 (MQ=255)
                                                                           tttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACGGATAGCCAATa                              >  2:3434455/1‑96 (MQ=255)
                                                                           tttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                          >  1:1527822/1‑100 (MQ=255)
                                                                           tttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                          <  1:3320257/100‑1 (MQ=255)
                                                                            ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                         <  1:3106583/100‑1 (MQ=255)
                                                                             tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                        <  1:3406692/100‑1 (MQ=255)
                                                                             tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                        <  1:760694/100‑1 (MQ=255)
                                                                              aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                       <  1:1182479/100‑1 (MQ=255)
                                                                               aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                      <  1:2254589/100‑1 (MQ=255)
                                                                               aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                      >  1:909157/1‑100 (MQ=255)
                                                                                gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                     <  1:2561747/100‑1 (MQ=255)
                                                                                 aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTTAATGCATAAACCGATAGCCAATACCGGTTTCAg                    <  2:2762395/100‑1 (MQ=255)
                                                                                 aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                    >  2:3120107/1‑100 (MQ=255)
                                                                                 aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                    <  1:3616429/100‑1 (MQ=255)
                                                                                 aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                    <  2:3639807/100‑1 (MQ=255)
                                                                                  tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                    <  2:2477405/99‑1 (MQ=255)
                                                                                  attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCCAAGCATAAACCGATAGCCATGACCGGGTTCAGt                   >  1:1989235/4‑100 (MQ=25)
                                                                                   ttaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTa                  >  2:451443/3‑100 (MQ=37)
                                                                                    aaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                 >  1:1732070/1‑100 (MQ=255)
                                                                                     aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt                >  1:880081/1‑100 (MQ=37)
                                                                                      aaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCCTAAACCGATAGCCAATACCGGTTTCAGTAATg               <  1:3976352/100‑1 (MQ=255)
                                                                                      aaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg               <  1:1603834/100‑1 (MQ=255)
                                                                                      aaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg               >  2:1076633/1‑100 (MQ=37)
                                                                                        gCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaa             >  1:4078738/1‑100 (MQ=255)
                                                                                            tatttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATgg         >  1:808295/1‑100 (MQ=255)
                                                                                            tatttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATgg         <  2:308742/100‑1 (MQ=255)
                                                                                              tttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGc        >  2:3857989/1‑99 (MQ=255)
                                                                                              tttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCg       >  1:35234/1‑100 (MQ=255)
                                                                                               ttattattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGt      <  2:1595569/100‑1 (MQ=255)
                                                                                                tattattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTTAGTAATGAAATGTCGTg     >  2:1976509/1‑100 (MQ=255)
                                                                                                  ttattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTggg   <  2:2217778/100‑1 (MQ=255)
                                                                                                   tattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgggc  <  2:3554940/100‑1 (MQ=255)
                                                                                                   tattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgggc  >  1:1747722/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: