Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,155,197 G→A Q394Q (CAG→CAA mdtA → multidrug efflux pump membrane fusion protein MdtA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,155,1970GA98.5% 235.8 / NA 66Q394Q (CAG→CAAmdtAmultidrug efflux pump membrane fusion protein MdtA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (28/36);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGTCGACACCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGTGGACTATCCGACCATTCAGGTGGTCACGCTCTA  >  NC_000913/2154957‑2155441
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
aGTCGACACCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAATTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCact                                                                                                                                                                                                                                             <  2:117681/250‑1 (MQ=255)
   cGACACCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTc                                                                                                                                                                                                                                         <  1:198355/251‑1 (MQ=255)
    gACACCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTcc                                                                                                                                                                                                                                        >  1:90651/1‑251 (MQ=255)
    gACACCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTcc                                                                                                                                                                                                                                        >  1:504578/1‑251 (MQ=255)
    gACACCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTcc                                                                                                                                                                                                                                        <  1:356601/251‑1 (MQ=255)
    gACACCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGGGGTGGAAGCCCAAAGCGCCACTACTcc                                                                                                                                                                                                                                        <  2:325548/251‑1 (MQ=255)
      cacCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCgg                                                                                                                                                                                                                                      <  2:374830/251‑1 (MQ=255)
             cAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGaa                                                                                                                                                                                                                               <  2:220988/251‑1 (MQ=255)
                   gCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGccac                                                                                                                                                                                                                         <  1:43139/251‑1 (MQ=255)
                        aGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCTGAAGAGAAAGCCACCAGc                                                                                                                                                                                                                     <  2:372624/250‑1 (MQ=255)
                         gTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCcg                                                                                                                                                                                                                   <  1:139695/251‑1 (MQ=255)
                            aTCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGa                                                                                                                                                                                                                <  1:305814/251‑1 (MQ=255)
                                    aGCCGCCCTGCAAATGGGCAATGAAGGCCATGTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGa                                                                                                                                                                                                        <  1:10727/251‑1 (MQ=255)
                                             gCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAgg                                                                                                                                                                                                 <  1:324941/249‑1 (MQ=255)
                                                aaggggCAATGAAGGACAGTTTGTCTGGGGGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCAATCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCAc                                                                                                                                                                                            <  2:586026/248‑1 (MQ=255)
                                                       aaTGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGCCAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGGGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGACAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTg                                                                                                                                                                                     >  2:276932/1‑251 (MQ=255)
                                                       aaTGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTg                                                                                                                                                                                     <  2:216942/251‑1 (MQ=255)
                                                                cATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCACAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGgtg                                                                                                                                                                            >  2:553366/1‑251 (MQ=255)
                                                                    ttGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTa                                                                                                                                                                         <  1:356980/250‑1 (MQ=255)
                                                                         ggggTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCg                                                                                                                                                                   <  2:530736/250‑1 (MQ=255)
                                                                                 gAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGagcagc                                                                                                                                                             <  2:137369/249‑1 (MQ=255)
                                                                                  aaTAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTTGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGcaca                                                                                                                                                          <  1:342650/251‑1 (MQ=255)
                                                                                     aGCGAAAACAAGGTCAGCAAGCATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCGCCTGAGGCAGGTGTTACCCCCGAGCAGCACAggc                                                                                                                                                       <  2:349033/251‑1 (MQ=255)
                                                                                      gCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAggcg                                                                                                                                                      >  2:454479/1‑251 (MQ=255)
                                                                                      gCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAggcg                                                                                                                                                      <  2:571778/251‑1 (MQ=255)
                                                                                       cGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAggcgg                                                                                                                                                     >  1:65009/1‑251 (MQ=255)
                                                                                       cGAAAACAAGGTCAGCAAACATCTGGTGACGCAGGGCATTCAGGACAGTCAGACAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGCCCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAggcgg                                                                                                                                                     >  2:475469/1‑251 (MQ=255)
                                                                                           aaCAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCg                                                                                                                                                 >  1:212569/1‑251 (MQ=255)
                                                                                             cAAGGTCAGCAAACAGCCGGTGACGCCGGGCAGTCAGGACAGTCAGAAAGTGGTGAGCCGTGCAGGTATTCCTGCGGGCTATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGCAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAGTACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTc                                                                                                                                               <  1:16609/251‑1 (MQ=255)
                                                                                                ggTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTcgcg                                                                                                                                            >  2:131091/1‑251 (MQ=255)
                                                                                                    aGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTg                                                                                                                                        <  1:143906/251‑1 (MQ=255)
                                                                                                     gCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGt                                                                                                                                       >  1:64308/1‑251 (MQ=255)
                                                                                                           aTCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTttatt                                                                                                                                  >  1:317412/1‑250 (MQ=255)
                                                                                                                gTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGt                                                                                                                            >  1:375879/1‑251 (MQ=255)
                                                                                                                  gACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTcc                                                                                                                          >  2:580508/1‑251 (MQ=255)
                                                                                                                    cGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCtg                                                                                                                        <  1:458670/251‑1 (MQ=255)
                                                                                                                       cGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTgg                                                                                                                     <  1:562609/251‑1 (MQ=255)
                                                                                                                        gggCATTCAGGACAGTCAGAAAGTGGTAATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATTCGTCCTGTGGc                                                                                                                    <  2:51397/251‑1 (MQ=255)
                                                                                                                                 ggACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACgcc                                                                                                           >  2:583034/1‑250 (MQ=255)
                                                                                                                                              gtggtgATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCg                                                                                              >  1:55730/1‑251 (MQ=255)
                                                                                                                                              gtggtgATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGAAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCg                                                                                              >  1:154193/1‑251 (MQ=255)
                                                                                                                                               tggtgATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCg                                                                                              >  1:489322/1‑250 (MQ=255)
                                                                                                                                                 gtgATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCCCGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATc                                                                                           <  1:309011/251‑1 (MQ=255)
                                                                                                                                                       cGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTc                                                                                     >  1:357488/1‑251 (MQ=255)
                                                                                                                                                       agTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCATGACCGCAGGGGCGAAAGTGGACGTGGTGGAAGCCCAAAGCGCCAATACTCCGGAAGAGAAAGCCACCAGCAGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCACGTCGCGCCTGTTTATTATGCGTCCTGTGGTCACAACGCTGCTGATGGTGGCGATCTTACTc                                                                                     >  1:357502/2‑251 (MQ=255)
                                                                                                                                                        gTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGTAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCg                                                                                    <  2:296365/251‑1 (MQ=255)
                                                                                                                                                        gTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTcc                                                                                    >  2:15039/1‑250 (MQ=255)
                                                                                                                                                          gCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGc                                                                                   <  1:385717/250‑1 (MQ=255)
                                                                                                                                                               tATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGAt                                                                             >  1:395953/1‑251 (MQ=255)
                                                                                                                                                                         ggCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTAt                                                                   <  2:99066/251‑1 (MQ=255)
                                                                                                                                                                               cgggtggtgACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGcc                                                             <  2:16870/248‑1 (MQ=255)
                                                                                                                                                                                  gtggtgACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCATGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGCGGCGATCTTACTCGCCGGGCTTATCGGTTATCGCGCCCTg                                                          >  2:174284/1‑251 (MQ=255)
                                                                                                                                                                                     gggACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACACCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACTCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGccc                                                       <  2:478212/249‑1 (MQ=255)
                                                                                                                                                                                                   tGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCGGCCGCGAATACGCGAACAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACCGGCGGCCCGTCGCGCCTGCTTATTATGCGTCCTGTGGCCACCACGCTGCTGTTGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGcc                                         >  2:289430/1‑251 (MQ=255)
                                                                                                                                                                                                       cGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGa                                      >  1:426449/1‑250 (MQ=255)
                                                                                                                                                                                                         ccTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGt                                   <  1:580686/251‑1 (MQ=255)
                                                                                                                                                                                                           tGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGTgg                                 >  1:118123/1‑251 (MQ=255)
                                                                                                                                                                                                               cGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGTGGACt                              >  2:13304/1‑250 (MQ=255)
                                                                                                                                                                                                                gAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAATCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCGGTTTATTATGCGGCCTGTGGCCACCACGCTGCGGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGTGGACTAt                            <  1:42544/251‑1 (MQ=255)
                                                                                                                                                                                                                         aaaGTGGAAGTGGTGGAAGCCCAAAGCGCCACTAAGCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTAGTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGTGGACTATCCGACCAtt                   <  1:18751/251‑1 (MQ=255)
                                                                                                                                                                                                                          aagtggaagtggTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGTGGACTATCCGACCATTc                  >  1:227119/1‑251 (MQ=255)
                                                                                                                                                                                                                           agtggaagtggTGGAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGTGGACTATCCGACCATTc                  >  1:105241/1‑250 (MQ=255)
                                                                                                                                                                                                                                aaggggtggAAGCCCAAAGCGCCACTACTCAGGAAGAGAAAGCCACCAGACGCGAAGACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGACGGAAGTGGACTATCCGACCATTCAggtgg            <  2:464217/247‑1 (MQ=255)
                                                                                                                                                                                                                                       ggAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGTGGACTATCCGACCATTCAGGTGGTCACGct     <  1:395262/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ggAAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGTGGACTATCCGACCATTCAGGTGGTCACGct     >  2:577324/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGGCCTGTTGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTAGCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGGGGACTATCCGACCATTCAGGTGGTCACGCTCTa  <  1:265687/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
AGTCGACACCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGTGGACTATCCGACCATTCAGGTGGTCACGCTCTA  >  NC_000913/2154957‑2155441

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: