Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,187,314 C→G 5.1% intergenic (+16/+66) rcnB → / ← yehA periplasmic modulator of Ni and Co efflux/putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,187,3140CG 5.1% 275.9 / 5.8 79intergenic (+16/+66)rcnB/yehAperiplasmic modulator of Ni and Co efflux/putative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base C (36/39);  new base G (2/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-03

ACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGG  >  NC_000913/2187044‑2187571
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                 
aCCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGa                                                                                                                                                                                                                                                               >  4:137794/1‑275 (MQ=255)
               aTTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGa                                                                                                                                                                                                                                                 <  2:468921/274‑1 (MQ=255)
               aTTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGa                                                                                                                                                                                                                                                 <  2:137795/274‑1 (MQ=255)
                                      tttCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGCCGGTGAGATTTTTTATCATCGCGAAAAAAAGCCCCCTCATCATGAGGGGGACATGCAGACCCCTTGTTATTTTTTCTTATTAg                                                                                                                                                                                                                         >  4:366280/1‑275 (MQ=255)
                                      tttCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGCCAGTTGCGTAACCTCCTCGCGCCTCATGCCGGGACGCACTGGAACTATCTGGGGGGCGCGTACGTGTTGATCAGCGACCCCGACGGTAAAAACATTAAAGCCTACGATGGTGAGATTTTTTATCATCTCTAAAAAAAACCCCCCCACCAAGGGGGGGAAATTCAAACTCCTTGGTactcttctcttata                                                                                                                                                                                                                          >  4:120634/1‑261 (MQ=255)
                                        tCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCCTGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGc                                                                                                                                                                                                                        >  2:70860/1‑274 (MQ=255)
                                        tCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCCGACACCTTGTTATTTTTTATTATTAGc                                                                                                                                                                                                                        >  2:7503/1‑274 (MQ=255)
                                          aTTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGa                                                                                                                                                                                                                                                 <  2:120635/247‑1 (MQ=255)
                                          aTTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTtta                                                                                                                                                                                                                                <  2:366281/264‑1 (MQ=255)
                                                                      gggACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTa                                                                                                                                                                                         >  2:2929/1‑275 (MQ=255)
                                                                               tACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACCCCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCAACATGAGGGGGAAATGAAGACCCCTTGTTATTTTTTATTATTAGCCACTTTCTCGGTTTTCTTGTTATTAGTTGTAtgttacttt                                                                                                                                                                                >  4:528680/1‑267 (MQ=255)
                                                                                        tGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACTGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATg                                                                                                                                                                        <  4:332226/274‑1 (MQ=255)
                                                                                         gtaccgtaccGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGGGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCg                                                                                                                                                                      <  4:202776/275‑1 (MQ=255)
                                                                                               taccgATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGcc                                                                                                                                                                >  2:303371/1‑275 (MQ=255)
                                                                                               taccgATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGcc                                                                                                                                                                >  4:362067/1‑275 (MQ=255)
                                                                                                   caTGGGTAAAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGGGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGGGCGTACGTGGTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCa                                                                                                                                                            <  4:61584/274‑1 (MQ=255)
                                                                                                         tACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTgcgc                                                                                                                                                      <  2:528681/275‑1 (MQ=255)
                                                                                                         tACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTgcg                                                                                                                                                       <  2:270553/274‑1 (MQ=255)
                                                                                                                   aGCAGTGGCCTTTTCGTTACCTGCCCGTGCCTGCTGCCCGGACTTACTGGACCTATAGTGGTGTCGCTTACGTTTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTTAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAGATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAAc                                                                                                                                            <  4:208854/275‑1 (MQ=255)
                                                                                                                   aGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTGAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATg                                                                                                                                                                        >  4:256417/1‑247 (MQ=255)
                                                                                                                   aGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAAc                                                                                                                                            >  2:454404/1‑275 (MQ=255)
                                                                                                                       gTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTGAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGt                                                                                                                                                                    <  2:256418/247‑1 (MQ=255)
                                                                                                                       gTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGttat                                                                                                                                                                                                 <  1:58653/218‑1 (MQ=255)
                                                                                                                          gCAGTTTCGTAACTGGCCCGGGCCTGATGCCGGGACGCACGTGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACtttgttt                                                                                                                                     <  4:397527/275‑1 (MQ=255)
                                                                                                                               tGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAgg                                                                                                                                                                                                                                                             >  4:284530/1‑150 (MQ=255)
                                                                                                                                   tAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGa                                                                                                                                                                                                                                                         <  2:284531/150‑1 (MQ=255)
                                                                                                                                               gcCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACt                                                                                                                 <  4:454403/274‑1 (MQ=255)
                                                                                                                                                    aTGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGt                                                                                                           >  2:526854/1‑275 (MQ=255)
                                                                                                                                                         ggtaCGCACTGGACCTATATGGGTGGCGCTTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATc                                                                                                       <  2:342519/271‑1 (MQ=255)
                                                                                                                                                         gggACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCCCGTTGATTAATGCGGTTGCCTCCAGTTCGCCAGATTTAACTTTTTTTGTATCGTAGACGTAGTAACTGGCGTTTATCg                                                                                                      >  4:524823/1‑275 (MQ=255)
                                                                                                                                                                       ctttatGGGTTGGGCGTACGTTGTGATCCGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGtt                                                                                        <  4:526853/272‑1 (MQ=255)
                                                                                                                                                                            tGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGTCTGTTATCGGAATTGAAACCGTTTTATc                                                                                   >  2:25018/1‑275 (MQ=255)
                                                                                                                                                                               gTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTTCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTg                                                                                 >  1:57603/1‑274 (MQ=255)
                                                                                                                                                                                tGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGt                                                                                >  2:438567/1‑274 (MQ=255)
                                                                                                                                                                                  gcgcGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCTTATACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTTTaa                                                                              >  1:118735/1‑274 (MQ=255)
                                                                                                                                                                                  gcgcGGACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTAATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCCCTTGCTCGTCTGGCTTGTTATTATTCGTATTTTACGTGGATTAATGCGTTTg                                                                                                                                                                  >  3:82050/1‑190 (MQ=255)
                                                                                                                                                                                    gcgTACGTGTTGATCAGCGACACCCACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTttattatta                                                                                                                                                                                                                          >  2:264884/1‑132 (MQ=255)
                                                                                                                                                                                        aCGTGTTGATCAGCGACACCAACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCa                                                                                                                                                                                                                      <  4:264883/132‑1 (MQ=255)
                                                                                                                                                                                            gtTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTc                                                                   >  4:233009/1‑275 (MQ=255)
                                                                                                                                                                                              tGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATGATGAggg                                                                                                                                                                                                                                                            >  6:249115/1‑88 (MQ=255)
                                                                                                                                                                                              tGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATGATGAggg                                                                                                                                                                                                                                                            <  5:249115/88‑1 (MQ=255)
                                                                                                                                                                                                          ccGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGAtt                                                                                                                                                                            >  3:115729/1‑156 (MQ=255)
                                                                                                                                                                                                           cGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTAttt                                                                                                                                                                                      >  4:135121/1‑145 (MQ=255)
                                                                                                                                                                                                              cGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCATGTTGATTAATg                                                                                                                                                                        <  1:115729/156‑1 (MQ=255)
                                                                                                                                                                                                               ggTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACg                                                                                                                                                                                  <  2:135122/145‑1 (MQ=255)
                                                                                                                                                                                                                      ttaTTTAAGCCTACGGCGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCTTGGGGGGGAAATGCAGACACCTTGTTATTTTTTATTTTTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTATCTGGCTGTTATCGTAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGtata                                         <  4:25018/273‑1 (MQ=255)
                                                                                                                                                                                                                      tCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCt                                                                                                                                                                                                                >  7:672930/1‑108 (MQ=255)
                                                                                                                                                                                                                      tCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCt                                                                                                                                                                                                                <  8:672930/108‑1 (MQ=255)
                                                                                                                                                                                                                      tCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGtata                                         <  4:303370/275‑1 (MQ=255)
                                                                                                                                                                                                                      tCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGtat                                          <  2:524824/274‑1 (MQ=255)
                                                                                                                                                                                                                        aTTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATc                                       <  2:233010/275‑1 (MQ=255)
                                                                                                                                                                                                                         ttAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCGGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTAttt                                                                                                                                                                                      >  2:97409/1‑131 (MQ=255)
                                                                                                                                                                                                                         ttAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTg                                                                                                                                                         >  4:212408/1‑160 (MQ=255)
                                                                                                                                                                                                                          tAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGa                                                                                                                                                                                                                                                               <  6:1008342/57‑1 (MQ=255)
                                                                                                                                                                                                                          tAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGa                                                                                                                                                                                                                                                               >  5:1008342/1‑57 (MQ=255)
                                                                                                                                                                                                                           aaaGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTg                                                                                                                                                                               >  1:80246/1‑136 (MQ=255)
                                                                                                                                                                                                                             aGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACg                                                                                                                                                                                  <  4:97408/131‑1 (MQ=255)
                                                                                                                                                                                                                             aGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGcc                                                                                                                                                     <  2:212409/160‑1 (MQ=255)
                                                                                                                                                                                                                              gCCTACTACTGTGGGATTTTTTTTCATCGCTCAAAAAATCCCCCTCATCATGAGGGTGAAATGCAGACTCCTTGTNTTTTTTTATTATTGGCCACTTGTTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCttt                                   <  3:118735/273‑1 (MQ=255)
                                                                                                                                                                                                                               ccTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTcttg                                                                                                                                                                                                         >  6:648253/1‑106 (MQ=255)
                                                                                                                                                                                                                               ccTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTcttg                                                                                                                                                                                                         <  5:648253/106‑1 (MQ=255)
                                                                                                                                                                                                                               ccTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTa                                                                                                                                                                           <  3:80246/136‑1 (MQ=255)
                                                                                                                                                                                                                                cTACGACGGGTGGATGTTTTGTCATCGCTAAAAAAAGCCCCCTCCACATGGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTActg                               <  4:438566/275‑1 (MQ=255)
                                                                                                                                                                                                                                 tACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTActgc                              >  2:210292/1‑275 (MQ=255)
                                                                                                                                                                                                                                               tttttATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGAttt                                                                                                                                               >  5:248515/1‑148 (MQ=255)
                                                                                                                                                                                                                                                    atcatcGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAg                                                                                                                                                                                              >  6:362847/1‑96 (MQ=255)
                                                                                                                                                                                                                                                    atcatcGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAg                                                                                                                                                                                              <  5:362847/96‑1 (MQ=255)
                                                                                                                                                                                                                                                       atcGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTAt                                                                                    >  3:53180/1‑199 (MQ=255)
                                                                                                                                                                                                                                                       atcGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTg                                             >  4:176002/1‑238 (MQ=255)
                                                                                                                                                                                                                                                         cGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACtttgttt                                                                                                                                     <  6:248515/148‑1 (MQ=255)
                                                                                                                                                                                                                                                           cTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGt                                                                                <  1:53180/199‑1 (MQ=255)
                                                                                                                                                                                                                                                           cTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGtata                                         <  2:176003/238‑1 (MQ=255)
                                                                                                                                                                                                                                                             aaaaaaaGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTgg  <  1:71239/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTaa                                                                                                                                                                          >  5:654184/1‑101 (MQ=255)
                                                                                                                                                                                                                                                                   aGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTaa                                                                                                                                                                          <  6:654184/101‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cccTCATGATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGt                                                                                                                                                                                             >  6:323151/1‑78 (MQ=255)
                                                                                                                                                                                                                                                                       cccTCATGATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGt                                                                                                                                                                                             <  5:323151/78‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cTCATCATGAGGGGGAAATGCAGACACCTTGTTAttt                                                                                                                                                                                                                                    >  5:702014/1‑37 (MQ=255)
                                                                                                                                                                                                                                                                         cTCATCATGAGGGGGAAATGCAGACACCTTGTTAttt                                                                                                                                                                                                                                    <  6:702014/37‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tcatcaTGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATcc                                      >  2:534919/1‑226 (MQ=255)
                                                                                                                                                                                                                                                                              catGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACgt                                                                                                                        <  5:639898/140‑1 (MQ=255)
                                                                                                                                                                                                                                                                              catGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACgt                                                                                                                        >  6:639898/1‑140 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                 
ACCAATACCTCGGGAATTTCTAAGTATGAGTTAAGTAGTTTCATTGCTGACTTTAAGCATTTCAAACCAGGGGACACCGTACCAGAAATGTACCGTACCGATGAGTACAACATTAAGCAGTGGCAGTTGCGTAACCTGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGG  >  NC_000913/2187044‑2187571

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: