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,6640GA97.3% 223.2 / NA 73G620R (GGA→AGA) yhcDputative fimbrial usher protein YhcD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (33/38);  minor base C (2/0);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

ATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGGGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGT  >  NC_000913/3364425‑3364908
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
ataCTCGCCATGACGATCATAAAACCCACTATGGTACTGGTTTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCGTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCttt                                                                                                                                                                                                                                           <  1:373460/251‑1 (MQ=255)
  aCTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAAAGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTc                                                                                                                                                                                                                                         >  2:535627/1‑251 (MQ=255)
  aCTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGGCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTc                                                                                                                                                                                                                                         <  2:63865/251‑1 (MQ=255)
     cGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGTTGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCa                                                                                                                                                                                                                                        <  1:191707/249‑1 (MQ=255)
        cATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGCAACGAGCTTGAACGGCAGCATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTCCAACGGCACGAGGTCAGGTGGTATTGCCGTACATGATCAGAGAGAGACCTTTTCACCGTg                                                                                                                                                                                                                                   >  2:20431/1‑251 (MQ=255)
               aTCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGCCTATc                                                                                                                                                                                                                            <  2:318873/251‑1 (MQ=255)
                  aTAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAAt                                                                                                                                                                                                                         >  2:474504/1‑251 (MQ=255)
                   aaaaaCCCACTATGGTAGTGGGGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAAGGCCGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACTCCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATg                                                                                                                                                                                                                        <  1:514488/250‑1 (MQ=255)
                     aaaCCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCATTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGTACTATCAATGac                                                                                                                                                                                                                      <  1:121023/251‑1 (MQ=255)
                      aaCCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGACTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGGGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGac                                                                                                                                                                                                                      >  2:465369/1‑250 (MQ=255)
                        cccACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGACACGAGCTTGCACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACt                                                                                                                                                                                                                   >  2:310452/1‑251 (MQ=255)
                                   aCTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCg                                                                                                                                                                                                         <  2:310911/250‑1 (MQ=255)
                                   aCTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGc                                                                                                                                                                                                        >  2:470463/1‑251 (MQ=255)
                                           cATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGCACGATATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTCCACCCAAGTCAGCCTTGCCGCAGGCGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCCAAAATgg                                                                                                                                                                                                >  1:304157/1‑251 (MQ=255)
                                                gAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAACATGGATAAc                                                                                                                                                                                           >  2:359151/1‑251 (MQ=255)
                                                    cgtcTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAata                                                                                                                                                                                       <  2:340734/251‑1 (MQ=255)
                                                        tcAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGc                                                                                                                                                                                   >  1:352989/1‑251 (MQ=255)
                                                         cAGATGAGATGTACTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCa                                                                                                                                                                                  <  1:135629/251‑1 (MQ=255)
                                                         cAGATGAGAGGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACGTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGATAGTGACCTTTTCACCGTGGACTATCAATGGCACTTTCGCCAGCGCAAAAATGGATAACAATATTGc                                                                                                                                                                                   <  2:381404/250‑1 (MQ=255)
                                                            atgagatgTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGgt                                                                                                                                                                               <  1:159065/251‑1 (MQ=255)
                                                                gatgTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAAGGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCa                                                                                                                                                                           <  2:560039/251‑1 (MQ=255)
                                                                 atgTCCTACTATGCGAATGATGCACGAGATCACGACGAACGTGAAAAGAACTTGAACGGCAGTATCAGTTCCACTCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCCAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATAGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAACTGGATAACAATAATGCAGGTGTccg                                                                                                                                                                          >  2:578081/1‑249 (MQ=255)
                                                                   ggccTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGa                                                                                                                                                                         <  2:169889/248‑1 (MQ=255)
                                                                      ctactaTGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTg                                                                                                                                                                                                                                   <  2:439240/189‑1 (MQ=255)
                                                                      ctactaTGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTg                                                                                                                                                                                                                                   >  1:439240/1‑189 (MQ=255)
                                                                        actaTGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTa                                                                                                                                                                     <  1:51470/249‑1 (MQ=255)
                                                                        actaTGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTAc                                                                                                                                                                    >  2:157663/1‑250 (MQ=255)
                                                                           aTGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTACGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACAt                                                                                                                                                                  <  2:2243/249‑1 (MQ=255)
                                                                             gtgAATGCTGAACGAGATCAAGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAAGCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTGTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTca                                                                                                                                                              <  2:320010/251‑1 (MQ=255)
                                                                              tgAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACAtctc                                                                                                                                                               >  1:162943/1‑249 (MQ=255)
                                                                                 aTGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGAAGGTGTCAGAATTACATCTcaggc                                                                                                                                                           >  1:482209/1‑250 (MQ=255)
                                                                                  ggCTGAACGATATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTGCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTTCGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCGTTTCACCGTGGACTATCAATGTCACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGGGTCAGAATTACAGCTcaggcag                                                                                                                                                         <  2:304157/250‑1 (MQ=255)
                                                                                                 cgacgaACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTAGCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGTCCAGTATGGACAg                                                                                                                                          <  1:141257/251‑1 (MQ=255)
                                                                                                    cgaACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGCTCAGAAAGTGACCTTTTCACCGTGGACTATCAATGCCACTTTCGCCATCGCAAAACTGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGCCAGAtt                                                                                                                                       >  1:434674/1‑251 (MQ=255)
                                                                                                     gaACGCGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCCGGTGGTATTGCCGTACATGATCAGCGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCCGGCCCAGTATGGACAGAttt                                                                                                                                      >  2:557911/1‑251 (MQ=255)
                                                                                                       aCGTGAAACGAGCTTGAACGGCAGTATCAGGTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTc                                                                                                                                    <  2:46456/251‑1 (MQ=255)
                                                                                                         gTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCgg                                                                                                                                  <  2:30678/251‑1 (MQ=255)
                                                                                                             aaCGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTa                                                                                                                              <  1:470463/251‑1 (MQ=255)
                                                                                                               cGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAACCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAAt                                                                                                                            >  2:26469/1‑251 (MQ=255)
                                                                                                                gAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATg                                                                                                                           >  2:141826/1‑251 (MQ=255)
                                                                                                                  gCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGcc                                                                                                                         <  2:488596/251‑1 (MQ=255)
                                                                                                                           gCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTcc                                                                                                                 >  1:494251/1‑250 (MQ=255)
                                                                                                                             aGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACTATGTCAGGTGGTATTGCCGTACATGATCAGAGATTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATc                                                                                                              <  1:94154/251‑1 (MQ=255)
                                                                                                                                aTCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAAt                                                                                                           >  1:6709/1‑251 (MQ=255)
                                                                                                                                  cAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACGCTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATcc                                                                                                         <  1:315229/251‑1 (MQ=255)
                                                                                                                                   aGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCa                                                                                                        <  2:162943/251‑1 (MQ=255)
                                                                                                                                    gTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGCCACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGACCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATcc                                                                                                         >  2:61786/1‑249 (MQ=255)
                                                                                                                                        cAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGt                                                                                                   <  1:439757/251‑1 (MQ=255)
                                                                                                                                                    cacCCAAGCCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTAAAACGGTACGATGTCCGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGACATCGCAAAACTGGATAACAATATTGCAGGTGTCAGAATTACATCTCCGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATGAATCCAGCCGTGGCGAACAtacg                                                                                       >  2:56558/1‑248 (MQ=255)
                                                                                                                                                      cccAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACTGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAgg                                                                                      >  1:321308/1‑250 (MQ=255)
                                                                                                                                                        cAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCCTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGt                                                                                   >  1:506864/1‑251 (MQ=255)
                                                                                                                                                         aaGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGtt                                                                                  >  1:108335/1‑251 (MQ=255)
                                                                                                                                                             cAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTgag                                                                               <  1:349020/250‑1 (MQ=255)
                                                                                                                                                              aGCCTTGCCGCAGGAGCACGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTg                                                                                 >  2:282455/1‑247 (MQ=255)
                                                                                                                                                                ccTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCg                                                                           <  2:235095/251‑1 (MQ=255)
                                                                                                                                                                    gCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCATGAGTTGAGATCGATa                                                                        <  1:157663/250‑1 (MQ=255)
                                                                                                                                                                        cAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGcc                                                                   >  2:112127/1‑251 (MQ=255)
                                                                                                                                                                         aGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATAcc                                                                      >  2:474827/1‑247 (MQ=255)
                                                                                                                                                                         aGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATAcc                                                                      <  1:474827/247‑1 (MQ=255)
                                                                                                                                                                           gAGCAGGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGATTATCAAGGACACTTTCGCCATCGCAAAAATGGATACCAATATTGCAGGTGTCAGAATTACATCCCAGGCAGGCCCAGTATGGACAGATTTGCGGGGTAATGCCGTCATTCCATCCATCCAGCCGTGGCGAACATCAGGAGTTGAGACCGATTCCGCCAGc                                                                <  2:529239/251‑1 (MQ=255)
                                                                                                                                                                                         gATAGCCGTACTTACAACGGTACGATGTCAGGTGGTGTTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGt                                                  <  2:47972/251‑1 (MQ=255)
                                                                                                                                                                                               cGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGatat                                            <  1:604379/251‑1 (MQ=255)
                                                                                                                                                                                                 tACTTACAACGGTACGGTGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCg                                          >  2:332836/1‑251 (MQ=255)
                                                                                                                                                                                                  aCTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCgg                                         >  1:360159/1‑251 (MQ=255)
                                                                                                                                                                                                        aaCGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTTAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACgg                                   <  1:535627/251‑1 (MQ=255)
                                                                                                                                                                                                        aaCGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACgg                                   >  1:127139/1‑251 (MQ=255)
                                                                                                                                                                                                              aCGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGGGACCCTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCTTCAGTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACaaa                             <  2:29719/251‑1 (MQ=255)
                                                                                                                                                                                                                 aTGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAAt                          <  1:68243/251‑1 (MQ=255)
                                                                                                                                                                                                                    tCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGAt                       >  1:148354/1‑251 (MQ=255)
                                                                                                                                                                                                                      aGGTGGTATTTCCGTACATGATCAGAGAGTGTCCTTTTCACCGTGGACTATCCATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACAGCAGGAGTTGAGATCGATCCCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCCCAAAAATGATCa                     <  2:223773/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ccGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCgt          >  1:108889/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aTGATCAGCGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCTGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGc     >  1:472238/1‑250 (MQ=255)
                                                                                                                                                                                                                                          aTCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAAAGGCACAAAAATGATCAAACAAGGCCGTGGGGCAga  >  1:26159/1‑249 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
ATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGGGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGT  >  NC_000913/3364425‑3364908

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