Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 720,946 C→T 100% intergenic (+106/+110) 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_000913720,9460CT98.9% 301.7 / ‑4.4 94intergenic (+106/+110)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (48/45);  total (49/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

AAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTA  >  NC_000913/720849‑721042
                                                                                                 |                                                                                                
aaGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAg                                                                                                <  1:2225517/100‑1 (MQ=255)
 aGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGt                                                                                               >  1:556113/1‑100 (MQ=255)
 aGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGt                                                                                               >  2:1750145/1‑100 (MQ=255)
  ggTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTa                                                                                              >  1:2141095/1‑100 (MQ=255)
   gTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTaa                                                                                             <  2:1257091/100‑1 (MQ=255)
   gTGAAGATATTCTATCCAAGAAATACTGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTaa                                                                                             <  2:1843633/100‑1 (MQ=255)
    tGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAat                                                                                            >  1:73650/1‑100 (MQ=255)
    tGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAat                                                                                            >  1:2019329/1‑100 (MQ=255)
     gAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAata                                                                                           >  2:1043421/1‑100 (MQ=255)
     gAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAAGTTATTTAATGATGAGCTTTTTACTTAGTAata                                                                                           <  1:2426384/100‑1 (MQ=255)
           atTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAat                                                                                     <  1:750143/100‑1 (MQ=255)
             tCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAatat                                                                                   >  2:541832/1‑100 (MQ=255)
              cTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTGATGATGAGCTTTTTACTTAGTAATATAAAATAtt                                                                                  >  1:1245415/1‑100 (MQ=255)
               tATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTg                                                                                 >  2:1334299/1‑100 (MQ=255)
               tATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTg                                                                                 >  1:335599/1‑100 (MQ=255)
                 nccAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGaa                                                                               >  1:847869/2‑100 (MQ=255)
                  ccAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGaa                                                                               <  1:1456810/99‑1 (MQ=255)
                    aagAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTg                                                                            <  2:1872188/98‑1 (MQ=255)
                      gAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTg                                                                            <  2:1482981/98‑1 (MQ=255)
                        aaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAt                                                                        >  2:586503/1‑100 (MQ=255)
                          tACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttt                                                                      >  2:850034/1‑100 (MQ=255)
                            cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttttt                                                                    >  1:2260905/1‑100 (MQ=255)
                             ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   >  2:2319122/1‑100 (MQ=255)
                             ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   >  1:659794/1‑100 (MQ=255)
                             ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   <  1:2187905/100‑1 (MQ=255)
                             ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   <  2:1122014/100‑1 (MQ=255)
                             ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   <  1:1289949/100‑1 (MQ=255)
                               ttGAACTGTGAATGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtg                                                                 <  2:179415/100‑1 (MQ=255)
                               ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtg                                                                 <  2:186877/100‑1 (MQ=255)
                                 gAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtgtg                                                               >  1:2129645/1‑100 (MQ=255)
                                 gAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtgtg                                                               >  2:891496/1‑100 (MQ=255)
                                  aaCTGTGAAGGTTATCAATAAAAAATTATCAATCTAAATTTATTTAATGATGATCTTTTTACTTAGTAATATAAAATATTGAATTGTTAtttt                                                                     <  2:516421/93‑1 (MQ=255)
                                  aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTgtgtgt                                                              >  2:953810/1‑100 (MQ=255)
                                   aCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgtt                                                             >  1:248086/1‑100 (MQ=255)
                                     tgtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgttgt                                                           >  2:1288219/1‑100 (MQ=255)
                                     tgtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgttgt                                                           <  1:2036568/100‑1 (MQ=255)
                                      ctgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTgttgtt                                                          <  1:2071778/99‑1 (MQ=255)
                                        gAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTgttgtta                                                         >  2:2071764/1‑98 (MQ=255)
                                        gAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTa                                                        >  1:505224/1‑100 (MQ=255)
                                          aGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAg                                                      >  2:802469/1‑100 (MQ=255)
                                            gTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATTAAATATTGAATTGTTAGTTTTGTGTGTTGTTTAAGa                                                     <  1:1804097/99‑1 (MQ=255)
                                            gTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGAt                                                    <  2:324999/100‑1 (MQ=255)
                                             tGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATa                                                   >  1:1014530/1‑100 (MQ=255)
                                             tGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATa                                                   >  1:186874/1‑100 (MQ=255)
                                             tGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATa                                                   <  1:1144516/100‑1 (MQ=255)
                                               aaaaaaaaaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATACAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATa                                                   <  2:1363352/92‑1 (MQ=255)
                                               aTCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaa                                                 >  1:1534039/1‑100 (MQ=255)
                                                tCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaa                                                >  1:2531430/1‑100 (MQ=255)
                                                 cAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                               >  1:2341510/1‑100 (MQ=255)
                                                   aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGc                                             <  1:2107307/100‑1 (MQ=255)
                                                     aaaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCg                                           <  1:1144260/100‑1 (MQ=255)
                                                     aaaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCg                                           <  1:400898/100‑1 (MQ=255)
                                                      aaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGt                                          <  1:1150493/100‑1 (MQ=255)
                                                       aaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGta                                         >  2:1313475/1‑100 (MQ=255)
                                                        aaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtat                                        >  1:2482116/1‑100 (MQ=255)
                                                        aaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtat                                        >  1:2482350/1‑100 (MQ=255)
                                                         aaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatt                                       <  1:210975/100‑1 (MQ=255)
                                                           tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattta                                     <  2:2111325/100‑1 (MQ=255)
                                                            gATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTTTTTTTTTTTTAAGATAAAAAGCCGtatttgt                                    >  2:1555654/1‑98 (MQ=255)
                                                            gATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAAATATAAAGCTGtatttat                                    >  2:2482552/1‑100 (MQ=255)
                                                             aTCACTCTTAATTTATTTAATGATGAGCTTTTTACTTATTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTAAAGATAAAAAGTCGtatttatt                                   >  2:725571/1‑100 (MQ=255)
                                                               cAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                 >  1:1527365/1‑100 (MQ=255)
                                                                 aTCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAATATATTGAATTGTTATTTTTTTGTGTTTATTATGATAAAAAGTCGTATttattatt                                >  2:347108/1‑99 (MQ=255)
                                                                 aTCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatta                               <  1:1832904/100‑1 (MQ=255)
                                                                  tCTTAATTTATTTAATGATGAGCTTTTTACTTATTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAc                              <  2:2486479/100‑1 (MQ=255)
                                                                   cTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                             <  2:2260891/100‑1 (MQ=255)
                                                                   cTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGGTTAAGATAAAAAGCCGTATTTAGTATTACg                             >  1:378159/1‑100 (MQ=255)
                                                                    ttAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTTTTTTTTGTTTAAGATAAAAAGCTGTATTTATTATTActt                            >  2:442800/1‑98 (MQ=255)
                                                                      aaTTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                          >  1:369598/1‑100 (MQ=255)
                                                                      aaTTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                          >  1:367681/1‑100 (MQ=255)
                                                                      aaTTTATTTAATGATGAGATTTTTACTTCGTAATATAAAATATTAAATTGTTATTATTGTGTGTTGTTTAATATAAAAAGCCGTATTTATTATTACGGCt                          <  2:2545556/100‑1 (MQ=255)
                                                                       atttatttAATGATGAGCTTTTTACTTAGTAATATAAAATATTGTATTGTTATTTTTGTGTTTTGTTTAAGATAAAAAGCCGTATTTATTATTGCGGc                           <  2:210192/98‑1 (MQ=255)
                                                                       atttatttAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTTTGTATTATTTAAGATAAAAAGCTGTATTTTTTGTTATGGttt                         >  2:2308609/1‑100 (MQ=255)
                                                                       atttatttAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                         <  2:2477691/100‑1 (MQ=255)
                                                                        tttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                        >  1:2052928/1‑100 (MQ=255)
                                                                         ttatttaATGATTAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTtta                       <  2:2162130/100‑1 (MQ=255)
                                                                         ttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAACCGTATTTATTATTACTGCTtta                       >  2:1144196/1‑100 (MQ=255)
                                                                           atttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttttta                     <  1:1673178/100‑4 (MQ=37)
                                                                             ttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatta                   <  1:850041/100‑1 (MQ=255)
                                                                              taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                  <  1:1293389/100‑1 (MQ=255)
                                                                              taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                  >  1:1507718/1‑100 (MQ=255)
                                                                               aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                 >  1:477628/1‑100 (MQ=255)
                                                                                atgatgAGCTTTTTACTTAGTAATATAAAATATTGATTTGTTATTTTTGTGTGATGTTTAAGATAAAAAAACGTATTTATTATTATGGCTttaatta                   >  2:2446885/1‑97 (MQ=255)
                                                                                atgatgAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                <  1:1178817/100‑1 (MQ=255)
                                                                                 tgatgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCAGTATTTATTATTACGGTTTTAATtaataa               >  2:2151592/1‑100 (MQ=255)
                                                                                  gatgaGCTTTTTACTTAGTAATATAAAATCTTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAACGCCGTATTTATTATTACGGCTTTAATTAATaaa              <  2:1537998/100‑1 (MQ=255)
                                                                                  gatgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaa              <  1:1163941/100‑1 (MQ=255)
                                                                                  gatgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaa              <  1:442573/100‑1 (MQ=255)
                                                                                  gatgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatcttaaa              <  2:2141081/100‑11 (MQ=21)
                                                                                    tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag            <  1:1043429/100‑1 (MQ=255)
                                                                                     gaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg           <  2:2052914/100‑1 (MQ=255)
                                                                                      aGCTTTTTACTTAGTAATATAAAATCTTGAATTCTTATTTTTTTGTCTTTTTTAAAATAACAAGCCGTATTTATTATTCCGGCTTTAATTAATAAAag            >  2:2023758/1‑98 (MQ=255)
                                                                                       gCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca         <  2:967485/100‑1 (MQ=255)
                                                                                              aCTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGAAAAAAAAAAAAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa  <  2:1995850/100‑1 (MQ=255)
                                                                                                 |                                                                                                
AAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTA  >  NC_000913/720849‑721042

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: