Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,004 T→A 16.1% 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,0040TA16.1% 468.4 / 73.6 180intergenic (+164/+52)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (56/95);  new base A (16/13);  total (72/108)
Fisher's exact test for biased strand distribution p-value = 9.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01

Only 100 of 191 total aligned reads displayed.

AAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCG  >  NC_000913/720905‑721103
                                                                                                   |                                                                                                   
    gATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTTTGTGTTGTTTAAGATAAAAAGCCGtatttat                                                                                                 >  2:3111228/1‑100 (MQ=255)
     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)
        aaTATTAATTTATTTAATGATGAGCTTTTTACTAAGTAACATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatt                                                                                             <  1:189970/100‑1 (MQ=255)
           cTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGTTAAAAAGCCGTATTTATTATTACg                                                                                          <  2:108440/100‑1 (MQ=255)
             tAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                        <  2:917596/100‑1 (MQ=255)
              aaTTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                                                                                       <  2:1051176/100‑1 (MQ=255)
               atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGCATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                      >  1:2981048/1‑100 (MQ=255)
                tttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                     <  2:1205571/100‑1 (MQ=255)
                 ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                    <  1:2446325/100‑2 (MQ=38)
                    tttaATGATGAGCTTTTTACTCAGTAATATAAAATATAGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttat                                                                                 <  2:3342550/100‑5 (MQ=37)
                     ttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatta                                                                                <  1:1601029/100‑1 (MQ=255)
                     ttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatta                                                                                <  1:2143353/100‑1 (MQ=255)
                     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)
                       aaTGATGAGCTTTTTATTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                              <  2:2475879/100‑1 (MQ=255)
                       aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                              <  1:3597887/100‑8 (MQ=37)
                       aaTGATGAGCTTTTTACGCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                              <  1:1433842/100‑1 (MQ=255)
                        atgatgAGTTTTTTACTCAATAATATAAAATATTGAATTTTTATTTTTTTGTCTTGTTTAATATAAAAAGCCTTATTTATTATTACGGCTTTAATtaata                                                                             <  2:1489872/100‑1 (MQ=255)
                        atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                             >  2:1597461/1‑100 (MQ=255)
                        atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATGTTTGTGGGTTGTTTAAGATAAAAAGCCGTATTTATTATTACCGCTTTAATtaata                                                                             >  2:1004054/1‑100 (MQ=255)
                          gatgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaa                                                                           >  1:328875/1‑100 (MQ=255)
                            tgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag                                                                         >  1:1399882/1‑100 (MQ=255)
                             gaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                        >  2:2640022/1‑100 (MQ=255)
                              aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc                                                                       <  1:2024450/100‑1 (MQ=255)
                                cTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag                                                                     <  2:3528466/100‑1 (MQ=255)
                                  ttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc                                                                   <  1:606352/100‑1 (MQ=255)
                                    ttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTg                                                                 >  1:1767703/1‑100 (MQ=255)
                                     tACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt                                                                <  2:11318/100‑1 (MQ=255)
                                      aCTCAGTAATATAAAATATTGAATTGTTATTTGTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa                                                               >  1:4054328/1‑100 (MQ=255)
                                         cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                            >  1:2933730/1‑100 (MQ=255)
                                         cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                            >  2:2561740/1‑100 (MQ=255)
                                         cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                            >  1:577910/1‑100 (MQ=255)
                                           gTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaa                                                           >  2:974802/1‑99 (MQ=255)
                                              atatAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAtt                                                       <  1:1834030/100‑1 (MQ=255)
                                               tataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTa                                                      <  1:2191466/100‑1 (MQ=255)
                                                ataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTTAAATTa                                                      >  2:777105/1‑99 (MQ=255)
                                                  aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                   <  2:1806412/100‑1 (MQ=255)
                                                     atatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTc                                                >  1:1063308/1‑100 (MQ=255)
                                                         tGAATTGTTATTTTTGTTTGTTGTTTAAGATAAATAGCCGTTTTTATTATTACGGCTTTAATTAATAAAAGGCAGTCTGTATTAAAATTAATATTCAAAg                                            <  2:221875/100‑1 (MQ=255)
                                                         tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                            <  1:2252533/100‑1 (MQ=255)
                                                         tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                            <  1:1169228/100‑1 (MQ=255)
                                                         tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                            <  1:2051533/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)
                                                             ttGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATa                                        <  1:97850/100‑1 (MQ=255)
                                                             ttGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATa                                        <  1:2096559/100‑1 (MQ=255)
                                                                ttATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                     >  1:799881/1‑100 (MQ=255)
                                                                 tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                    <  2:1181294/100‑1 (MQ=255)
                                                                 tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                    >  2:4107172/1‑100 (MQ=255)
                                                                  aTTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCg                                   >  1:3870495/1‑100 (MQ=255)
                                                                    ttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAt                                 <  2:4159719/100‑1 (MQ=255)
                                                                     tttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAt                                 >  1:433070/1‑99 (MQ=255)
                                                                     tttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                <  1:3766515/100‑1 (MQ=255)
                                                                     tttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                <  2:726559/100‑1 (MQ=255)
                                                                     tttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                <  2:3752914/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)
                                                                          gtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCaa                           <  1:2303004/100‑1 (MQ=255)
                                                                          gtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCaa                           >  2:2867365/1‑100 (MQ=255)
                                                                          gtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCa                            <  2:470797/99‑1 (MQ=255)
                                                                           tgttgtTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAAt                          >  2:2434085/1‑100 (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:3320257/100‑1 (MQ=255)
                                                                                 ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                    <  1:3106583/100‑1 (MQ=255)
                                                                                 ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                    <  2:1747269/100‑1 (MQ=255)
                                                                                 ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAAGTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                    >  2:2227976/1‑100 (MQ=255)
                                                                                 ttAAAATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                    <  2:3755729/100‑1 (MQ=255)
                                                                                  tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                   <  1:760694/100‑1 (MQ=255)
                                                                                  tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                   <  1:3406692/100‑1 (MQ=255)
                                                                                   aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                  <  1:1182479/100‑1 (MQ=255)
                                                                                    aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                 >  1:909157/1‑100 (MQ=255)
                                                                                    aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                 <  2:1335288/100‑1 (MQ=255)
                                                                                     gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                <  1:2561747/100‑1 (MQ=255)
                                                                                      aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTTAATGCATAAACCGATAGCCAATACCGGTTTCAg               <  2:2762395/100‑1 (MQ=255)
                                                                                      aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg               <  2:3639807/100‑1 (MQ=255)
                                                                                      aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg               >  2:3120107/1‑100 (MQ=255)
                                                                                      aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg               <  1:3616429/100‑1 (MQ=255)
                                                                                       tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt              >  2:1960380/1‑100 (MQ=255)
                                                                                       attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt              >  1:15302/4‑100 (MQ=37)
                                                                                        ttaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTa             >  2:451443/3‑100 (MQ=37)
                                                                                          aaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt           <  2:1833511/100‑1 (MQ=255)
                                                                                          aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt           >  1:880081/1‑100 (MQ=37)
                                                                                          aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt           >  2:1307748/1‑100 (MQ=37)
                                                                                           aaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCCTAAACCGATAGCCAATACCGGTTTCAGTAATg          <  1:3976352/100‑1 (MQ=255)
                                                                                             gCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaa        <  1:1741926/100‑1 (MQ=255)
                                                                                             gCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaa        >  1:4078738/1‑100 (MQ=255)
                                                                                                gTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATg     >  1:1674430/1‑100 (MQ=255)
                                                                                                 tatttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATgg    <  2:308742/100‑1 (MQ=255)
                                                                                                   tttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCg  >  1:35234/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCG  >  NC_000913/720905‑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: