Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,364,664 G→A G620R (GGA→AGA)  yhcD → putative fimbrial usher protein YhcD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,364,6640GA98.7% 233.1 / ‑4.6 78G620R (GGA→AGA) yhcDputative fimbrial usher protein YhcD
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (40/37);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

ACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGGGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGT  >  NC_000913/3364418‑3364903
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
ccacggTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCGCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGt                                                                                                                                                                                                                                             <  1:470329/246‑1 (MQ=255)
  acTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCagagagagc                                                                                                                                                                                                                                           >  2:492072/1‑248 (MQ=255)
            cGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGACTGCTGAACGAGATCACGACGACCGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCTTACATGATCAGAGAGTGACCTTTTCAcc                                                                                                                                                                                                                                 >  1:242449/1‑251 (MQ=255)
                    cGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTa                                                                                                                                                                                                                         <  2:659027/251‑1 (MQ=255)
                         aTAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAAt                                                                                                                                                                                                                    >  1:363081/1‑251 (MQ=255)
                         aTAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAAt                                                                                                                                                                                                                    <  1:635179/251‑1 (MQ=255)
                         aTAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAAt                                                                                                                                                                                                                    <  1:650212/251‑1 (MQ=255)
                            aaaCCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGACGTCCTACTATGTGAATGCTGAAAGAGATCACGACGAACGTGAAACGAGCTTGAACGGCGGTATAAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGAACTTCCAACGGTCCGATGTCAGGTGGTATTGCAGTACATGATCGGAGAGTGACCTGTTCACCGTGGACTATCAATgcc                                                                                                                                                                                                                 >  2:390613/1‑249 (MQ=255)
                                    tATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTACTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAACCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGc                                                                                                                                                                                                         >  1:399846/1‑251 (MQ=255)
                                     agggTACTGGTGTCATGGGAGTCGTCTCAGATGAGAAGTCCTACTATGTGAATGCTGAACGAGATCGCGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGGTCCAATCTCCATTACAACCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGcc                                                                                                                                                                                                        <  2:204364/249‑1 (MQ=255)
                                     aTGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGcc                                                                                                                                                                                                        <  1:524573/251‑1 (MQ=255)
                                       ggaactggcggcaggggAGACGTCTCAGATGAGATGTCCTTCTAGGTGAATGCTGAACTAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGGTTCAGGTCCAATCTCCATTACACCCAATTCAGCCTTGCCTAAGGAGCAAGAGGCAGTGATGGCCGTACTTGCATAGTTACGATGTCAGGTGGTATTGTAGTACATGATCAGAGAGTTACCTTGTCACCGTGGACTATCTATGACACTTTCGCCa                                                                                                                                                                                                       <  2:92752/236‑1 (MQ=255)
                                       ggTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCa                                                                                                                                                                                                       <  1:21035/250‑1 (MQ=255)
                                       ggTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGGGGACTATCAATGACACTTTCGCCa                                                                                                                                                                                                       >  2:21035/1‑250 (MQ=255)
                                        gTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATc                                                                                                                                                                                                     <  1:145255/251‑1 (MQ=255)
                                               tgtCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGCACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCTGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCaaaa                                                                                                                                                                                               >  1:214333/1‑250 (MQ=255)
                                                 tCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGATGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTCCGATGTCAGGTGGAATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCCATGACACTTACGCAATCGCAAAAATg                                                                                                                                                                                            >  2:580876/1‑251 (MQ=255)
                                                                cAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGCTAACAATATTGc                                                                                                                                                                              >  1:134754/1‑250 (MQ=255)
                                                                   atgagatgTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGgt                                                                                                                                                                          <  1:240547/251‑1 (MQ=255)
                                                                          gtcttactaTGTGAATGCTGAACGAGATAACGACGAACGCTAAAAGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCTGGAGCAAGCGGCAGTGATAGCCGTACTTATAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCTCCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATACCAATAGTGCAGGTGTCAGa                                                                                                                                                                    <  1:521501/246‑1 (MQ=255)
                                                                            ccTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAAt                                                                                                                                                                  <  1:512291/250‑1 (MQ=255)
                                                                               actaTGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACa                                                                                                                                                              <  1:492072/251‑1 (MQ=255)
                                                                                ctaTGTGAATGCTGAACGAGATCACGACGAACGTGAAAAGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGCGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGAGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTAcct                                                                                                                                                             >  2:539149/1‑249 (MQ=255)
                                                                                  aTGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACAtc                                                                                                                                                            <  1:465581/250‑1 (MQ=255)
                                                                                          gCTGACCGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTcaggca                                                                                                                                                     >  2:109712/1‑249 (MQ=255)
                                                                                           cTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGGATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCcggcaggc                                                                                                                                                  >  2:537790/1‑251 (MQ=255)
                                                                                             gAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGccc                                                                                                                                                >  1:560521/1‑251 (MQ=255)
                                                                                             gAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGcc                                                                                                                                                 <  1:544506/250‑1 (MQ=255)
                                                                                                  agaTCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTAt                                                                                                                                           >  1:394052/1‑251 (MQ=255)
                                                                                                  agaTCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGGTCCAAGCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTAt                                                                                                                                           <  2:134754/251‑1 (MQ=255)
                                                                                                        cgacgaACGTGAAACGAGCTTGAACGGCAGTATCAGTTACAATCTCCATTACACCAAAATCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGTTGTCAGGTGGGATTGCCGTACATGATCCGAGAGTGACCTTTTCACCGTGGACTATAAATGACACTTTCGCCATCGCAAAAATGGATAACAATCTTGCAGGTGTCAGAATTACATCTCAGGCAGGccacgaatggaca                                                                                                                                      >  2:420368/1‑239 (MQ=255)
                                                                                                           cgaACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCCGGTGGTATTGCCGTACATGATCAGAGAGAGACCTTTTCACCGTGGACTATCAATGACCCTTTCGCCATc                                                                                                                                                                                                     >  2:260798/1‑184 (MQ=255)
                                                                                                           cgaACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATc                                                                                                                                                                                                     <  1:260798/184‑1 (MQ=255)
                                                                                                             aaCGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGAtttt                                                                                                                                >  1:243548/1‑251 (MQ=255)
                                                                                                                     aCGAGCTTGAACGGCAGTATCCGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCCGTGATAGCCGTACTTACAACGGTACGCTGTCAGGTGGTATTGCCGTACATGATCCGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCCAAAATGGAAAACAATATTTAAGGAGTCAGACTTACATCTCCGGCCAGCCCAGTATGGACCGATTTTCGGGGTac                                                                                                                        >  2:222810/1‑250 (MQ=255)
                                                                                                                                   cAGTATCAGTGCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTATCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCAt                                                                                                          <  1:669657/251‑1 (MQ=255)
                                                                                                                                     gTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCa                                                                                                        <  1:539149/251‑1 (MQ=255)
                                                                                                                                     gTATCAGTTCCAATCTCCATTACACCCAAGTCAGACTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATc                                                                                                         >  1:170365/1‑250 (MQ=255)
                                                                                                                                       aTCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGa                                                                                                                                                                                                                            <  2:453090/133‑1 (MQ=255)
                                                                                                                                       aTCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGa                                                                                                                                                                                                                            >  1:453090/1‑133 (MQ=255)
                                                                                                                                        tCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATc                                                                                                     >  1:245742/1‑251 (MQ=255)
                                                                                                                                                          acacCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGGACATc                                                                                    >  2:182519/1‑250 (MQ=255)
                                                                                                                                                          acacCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCa                                                                                   <  1:209578/251‑1 (MQ=255)
                                                                                                                                                              ccAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGAACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAgg                                                                                 >  1:340384/1‑249 (MQ=255)
                                                                                                                                                                  gTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCTGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTga                                                                           <  1:537239/251‑1 (MQ=255)
                                                                                                                                                                  gTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTga                                                                           <  1:214379/251‑1 (MQ=255)
                                                                                                                                                                       ccTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAACTGGATAACCATATTGCAGGTGTCAGAATTCCATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGCACATCAGGAGTTGAGATcc                                                                      >  1:584668/1‑250 (MQ=255)
                                                                                                                                                                           gCCGCAGGAGCAAGCGGCAGTGATAGCCGGACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCCTTTCCCCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATAc                                                                  <  1:373453/251‑1 (MQ=255)
                                                                                                                                                                            ccGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATAcc                                                                 >  1:457780/1‑251 (MQ=255)
                                                                                                                                                                              gCAGGAGCAAGCGGCAGTGATAGCCGTACTTACCACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACAtc                                                                                                                                                            >  1:540088/1‑158 (MQ=255)
                                                                                                                                                                              gCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACAtc                                                                                                                                                            >  1:49110/1‑158 (MQ=255)
                                                                                                                                                                              gCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGAGAACAATATTGCAGGTGTCAGAATTACAtc                                                                                                                                                            <  2:540088/158‑1 (MQ=255)
                                                                                                                                                                              gCAGGAGCAAGCGGCAGGGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGGCAGAATTACAtc                                                                                                                                                            <  2:49110/158‑1 (MQ=255)
                                                                                                                                                                                 ggAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGcc                                                              >  1:334119/1‑249 (MQ=255)
                                                                                                                                                                                      aaGCGGCAGTGATAGCCGTACTTAAAACGAAACGATGTCACGTGGTAATGCCGTACATGATAAGAGAGTGACCTTTTCACCGTGGACTATCAAAGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCACCGTATGGACAGATTTTCGGGGAAATGCCGTCATTCCATCAATCCAGCAGTGGCGAACATCAGGAGTTGAGATCGATACCGACAgatagc                                                       >  2:481338/1‑246 (MQ=255)
                                                                                                                                                                                       aGCGGCAGTGATAGCCGCAATTACAAAGGTACGATGCCAGGGGGCATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCCAAAATGGATAACAATATTGCAGGTGTCAGTATTACATCTCCGGCAAGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCCGGGGTTGAGGTCGATTCCGCCCGCtatcc                                                      >  2:31652/1‑247 (MQ=255)
                                                                                                                                                                                        gCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCATAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCa                                                     <  2:243548/251‑1 (MQ=255)
                                                                                                                                                                                             aGTGTTAGCAGTACTTACAACGGTTCGATGTCAGGTGGTATGGCCGTACATGATCAGAGAGTGGCCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCCGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGATTTGAGATCGATACCGCCAGCTTGCCaaaaaa                                                <  2:373028/251‑1 (MQ=255)
                                                                                                                                                                                                gATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATg                                              <  1:187435/250‑1 (MQ=255)
                                                                                                                                                                                                       gTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGGGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAACTTGCCAAAAAATGTCGATATCg                                     <  1:180308/251‑1 (MQ=255)
                                                                                                                                                                                                       gTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCACTCCAGCCGTGGCGAACATCAGGAGTTGAGATCGTTACCGCCAGCTTGCCAAAAAATGTCGTTATc                                      >  1:97790/1‑251 (MQ=255)
                                                                                                                                                                                                             acacaggtaccattttagttggtATTGCCGTTCATGTTCAAAGAGTGCACTTTTCACCGTTGACTATCAATGACACTGTCGCAATCGCAAAAATGGATAACAAGATTGCAGGTGTCAGAATTCCATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCTATACCGCCAGCTTGCTAAAAAATGTCGATATCGGTaaa                                <  2:368134/232‑2 (MQ=255)
                                                                                                                                                                                                                  ggTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTa                                  >  2:562782/1‑244 (MQ=255)
                                                                                                                                                                                                                  ggTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTa                                  <  1:562782/244‑1 (MQ=255)
                                                                                                                                                                                                                   gTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTTTCAATGACACTGTCGCCATCGCAAAAATGGATAACAATATGGCAGGTGTCAGAATTTCATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGcaca                          <  2:334119/251‑1 (MQ=255)
                                                                                                                                                                                                                    tACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGcaca                          >  1:82459/1‑250 (MQ=255)
                                                                                                                                                                                                                       gATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTaa                                                                                                                        >  2:323961/1‑153 (MQ=255)
                                                                                                                                                                                                                       gATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTaa                                                                                                                        <  1:323961/153‑1 (MQ=255)
                                                                                                                                                                                                                       gATGTCAGGTGGTATTGACGTACATGATCAGAGAGTGACCTTTTCACCGTGGCCTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCCGGTGTCAGAATTACATCTCAGACAGGCCCAGTATGGACAGCTTTTAGGGGTAATGCCGTCATTCCATCATTCCAGCCGTGGCGAACATCAGGAGTTGAGATAGAGCCCGACCGCTTGCAAAAAAATGTCGATATCGTTAACGGCACaaaa                       >  2:676166/1‑250 (MQ=255)
                                                                                                                                                                                                                          gTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCACAACTGGATAACAATATTGCAGGTGTCAGAATTACATCTCAAGCAGGCCCAGTATGGACAGATTTTAGGGGAAATGCCGTCATTCCATCAGTCCAGCCGTGGCGTACATCAGGAGATGAGATCGATACCGCCAACTTGCCCACAATTGACGATGTCGGTAAAGGCACAAAAATg                    >  2:510823/1‑250 (MQ=255)
                                                                                                                                                                                                                          gTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATg                    >  2:62766/1‑250 (MQ=255)
                                                                                                                                                                                                                          gTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGa                   <  1:109811/251‑1 (MQ=255)
                                                                                                                                                                                                                             aggaggtACTGCCGTACATGATCAGAGCGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGCACAGATTTTCGGGGCAATGCCGTCATTCCATCAATCCTGCCGTGGCGAACATCAGGCGTTGAGATCGATACCGCAAGCTTGCCCAAAAGTGACGTTATCGGTAACGGCAAAAAAATGATc                 >  2:143298/5‑250 (MQ=255)
                                                                                                                                                                                                                              ggtggtATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCa                >  1:162174/1‑250 (MQ=255)
                                                                                                                                                                                                                                 ggtATGGCCGTACATGATCAGAGAGTGACCTTTTCACCGGGGACTATTAATGACCCTTTAGCAAGCGCAAAAATGGATACCAATATTGCAGGTGTTAGAATTACATCTCAGGCAGGCCAAGTATGGACAGATTTTCGGGGTAATGCGGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAAGGTCGATATCGGTAACGGCACAAAAATGATCAAACa            <  2:560521/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aTTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAgg         <  1:345768/251‑1 (MQ=255)
                                                                                                                                                                                                                                           tACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGAt                                                                    >  1:499378/1‑185 (MQ=255)
                                                                                                                                                                                                                                           tACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCgtggg  >  2:218153/1‑250 (MQ=255)
                                                                                                                                                                                                                                           tACATGATCAGAGAGTGACCTTTGCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGAt                                                                    <  2:499378/185‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
ACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGGGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGT  >  NC_000913/3364418‑3364903

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