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,6640GA96.3% 162.3 / ‑5.8 54G620R (GGA→AGA) yhcDputative fimbrial usher protein YhcD
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (23/29);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGGGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCC  >  NC_000913/3364423‑3364898
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                          
gTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGAcc                                                                                                                                                                                                                                    <  2:340282/250‑1 (MQ=255)
       cGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCAcc                                                                                                                                                                                                                            <  1:344074/251‑1 (MQ=255)
           aTGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGGGAGTGACCTTTTCACCGTgg                                                                                                                                                                                                                        >  2:240793/1‑251 (MQ=255)
             gACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGGTGTCTCAGATGAGATGTTCTACGATGTGAAGGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACTCCCAAGTCAGCCGTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTCCAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTgg                                                                                                                                                                                                                        <  1:282956/249‑1 (MQ=255)
                 aTCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCATGCCGCAGGAGCAAGCGGCAGCGATAGCCGTCCTTACAACGGTACGATGTCAGGTGGGATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTAt                                                                                                                                                                                                                   >  2:175879/1‑250 (MQ=255)
                               tATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACCTGATCAGAGAGTGACCTTTTCAcccgatcggacgagcgtcgtgtcggg                                                                                                                                                                                                    >  2:128823/1‑227 (MQ=255)
                               tATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCAcc                                                                                                                                                                                                                            <  1:128823/227‑1 (MQ=255)
                               tATGGTACTGGAGTCATGGGAGTCGTCTCAGATGAGCTGTCCTACTATGCGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTCTCAGTTCCAATCTCCATTACACCCAAGCCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACCATGTCAGGTGGGATTGCCGGACCTGATCAGAGAGTGACCTTTTCACCGTGGTCTATCAATgcaactctcgc                                                                                                                                                                                                    >  2:180775/1‑241 (MQ=255)
                                   gTACGGGTGTCATGGGAGTGTGCTCAGATGAGATGTCATACTATTTGAATGTTGAACGAAATCACGACGAACGGGAAACGAGCTTGAACGGCAGTATCAGTCCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGATAGTGACCGGGTCACCGTGGACTATCAATGACACTTTCTCCATc                                                                                                                                                                                                <  2:452436/251‑1 (MQ=255)
                                              aTGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTAGTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGa                                                                                                                                                                                     <  1:80405/251‑1 (MQ=255)
                                                 ggAGTCGTCTCAGATGAGATGACATACTAAGTGAATGCTGAACGAGATCACGCCGCACGTGAACCGAGCTTGAACGGCAGAATCAGTTCCAATCTCCATTACACCCAAGTCAACCTTGCCGCAGGAGCAAGCGGCAGTGATATCCGTACTTACAACGGGACGCTGTCATGAGGTATTGCCGAACATGGTCAAAGAGTGCCCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCatccatgtata                                                                                                                                                                                   >  2:400625/1‑240 (MQ=255)
                                                       gtcTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAatat                                                                                                                                                                            >  1:162071/1‑251 (MQ=255)
                                                        gtcccAGAGGAGATGTCCTACTATTTGAATGCTGAACGAGAACACGACGAACGTGAAACGAGCTTGAACGGCAGAATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTTGTATTGCCGTACATGATCAGAGAGTTACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATTGATAACAatat                                                                                                                                                                            <  2:327291/247‑1 (MQ=255)
                                                             agatgagatgTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAg                                                                                                                                                                       <  1:31083/251‑1 (MQ=255)
                                                                   gatgTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCATTGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGTTGTc                                                                                                                                                                  >  2:258117/1‑250 (MQ=255)
                                                                      gTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGa                                                                                                                                                               >  2:91713/1‑250 (MQ=255)
                                                                       tCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAAt                                                                                                                                                             >  1:116848/1‑251 (MQ=255)
                                                                           actaTGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTa                                                                                                                                                           <  1:286229/249‑1 (MQ=255)
                                                                           actaTGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTAc                                                                                                                                                          <  1:224139/250‑1 (MQ=255)
                                                                            ctaTGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCCGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACAt                                                                                                                                                        >  1:385270/1‑251 (MQ=255)
                                                                                 tgAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGCCTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGACTTACAtctc                                                                                                                                                     >  2:466869/1‑249 (MQ=255)
                                                                                    aTGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTcaggc                                                                                                                                                 >  1:300373/1‑250 (MQ=255)
                                                                                            cGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGGCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATAGTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGt                                                                                                                                        <  2:380459/251‑1 (MQ=255)
                                                                                            cGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGGATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGt                                                                                                                                        <  2:383079/251‑1 (MQ=255)
                                                                                                 tCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGATCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGTTGTCAGGTGGTATTGCCGTACATGATCAGAGCGAGACCTTTTCACCGTGGACTATCATTGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAATATgg                                                                                                                                    >  2:500849/1‑250 (MQ=255)
                                                                                                   acgacgAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGGTGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATCTTGCAGTTGTAAAAATTAAATCTCAGGCAGGCCCCGTATggca                                                                                                                                  >  2:287558/1‑248 (MQ=255)
                                                                                                     gacgaACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGa                                                                                                                               <  1:139140/251‑1 (MQ=255)
                                                                                                              gAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCggg                                                                                                                       <  1:98050/250‑1 (MQ=255)
                                                                                                               aaaCGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCAACGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGCTTTCCggggg                                                                                                                     >  2:431795/1‑250 (MQ=255)
                                                                                                                        tGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGt                                                                                                             <  1:499615/250‑1 (MQ=255)
                                                                                                                          aaCGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCAt                                                                                                          <  2:116848/251‑1 (MQ=255)
                                                                                                                          aaCGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCAt                                                                                                          <  1:395047/251‑1 (MQ=255)
                                                                                                                           aCGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTAAGCCTTGCCGCAGCAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGAGGTATTGCCGTACATGATCAGAGAGAGACCTTTTCACCGAGGACTATCAATGACACTTTCGCCATCGCAAAAATGGCTAACAATATTGCAGGCGTCAGAATTACATCTCAGGCAGGCCAAGTACGGACAGATTTTCGGGGTGATGACGTTAtt                                                                                                         >  1:63074/1‑251 (MQ=255)
                                                                                                                             ggCAGTATCAGTTCCAATTTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTcc                                                                                                       <  2:207915/251‑1 (MQ=255)
                                                                                                                                    tCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGAAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATc                                                                                                >  2:460316/1‑251 (MQ=255)
                                                                                                                                            aaTCTCCATCACACCCAAGCCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACATACAACGGTACTATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTGTCACCGTGGACTATCAAGGACACGTTCGCCATCGCAAAAATGGATAACAATATGGCAGGTGTCAGAAGTACATCTCAGGCAGGCCCAGGATGGACAGATTTTCTGGGTAATGCCGTCATTCCATCAATCCAGCCGTg                                                                                        <  2:176757/251‑1 (MQ=255)
                                                                                                                                                 ccATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTgg                                                                                       <  1:442731/247‑1 (MQ=255)
                                                                                                                                                                cAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTgag                                                                     <  1:211045/250‑1 (MQ=255)
                                                                                                                                                                  gCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATc                                                                  <  1:287558/251‑1 (MQ=255)
                                                                                                                                                                                 cAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTg                                                   <  1:129416/251‑1 (MQ=255)
                                                                                                                                                                                        cAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCaaaaa                                            <  1:213314/251‑1 (MQ=255)
                                                                                                                                                                                          gTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCCGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAAt                                          >  1:480318/1‑251 (MQ=255)
                                                                                                                                                                                          gTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCaaaaaa                                           >  2:107996/1‑250 (MQ=255)
                                                                                                                                                                                             aTAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCACTCAAGCCGTGGCGAACATCAGGGGTTGAGATCGATACCGCAAGCTTGCCAAAAAATgcc                                       >  2:202362/1‑249 (MQ=255)
                                                                                                                                                                                                gCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGGGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGTCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGat                                    <  1:240793/251‑1 (MQ=255)
                                                                                                                                                                                                         aCAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAAc                           <  1:258117/251‑1 (MQ=255)
                                                                                                                                                                                                          cAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGTGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACg                          >  1:10442/1‑251 (MQ=255)
                                                                                                                                                                                                              ggTACGATGTCAGGTGGGATTGCCGTACATGATCAGAGAGTGTCCTTTTCACCGTGGACTATCAATGTCACTTTCGCCACCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGcac                      <  2:499630/251‑1 (MQ=255)
                                                                                                                                                                                                                     tGTCAGGCGGGATTGCCGTACATGATCAGAGAGCGACCTTTTCACCGTGGACTATCACTGCCACTTTGGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGACTTACATCTCAGGCAGGCCCAGTATGGACAGATTTACGGGGTAATGCCGTCTTTCCATCACTCCAGACGTGGCGAACATCAAGAGTTGAGATCGATACCGCCAACTTGCCAAAAAATGACGGTATCGGTAACGGCACaacattg               >  2:112313/1‑246 (MQ=255)
                                                                                                                                                                                                                        cAGGTGGTATTGCCGTACATGATCTGAGAGTAACCCTTTTACCGAGGACTATCAATGACCCTTCCGCCATCGCAAAAATGGATAAAAATATTGCAGGTGGCAGAATTACATCGCAGGCAGGCCCAGTATGGACACATTTTCGGGGTAAGTCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCAGTCACGGCACAAAAATGAt             <  2:22498/250‑1 (MQ=255)
                                                                                                                                                                                                                           gtggtATTGCCGTACATGATCAGAGAGTGACCTTTTTACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCaaa         <  1:51494/251‑1 (MQ=255)
                                                                                                                                                                                                                              gtATTGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTGCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACaa      <  1:460316/251‑1 (MQ=255)
                                                                                                                                                                                                                                  tGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGCTATCGGTAACGGCACAACAATGATCAAACAAGGcc  >  1:259464/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tGCCGTACATGATCAGAGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCAAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGcc  >  1:27476/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                          
GTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACATGATCAGGGAGTGACCTTTTCACCGTGGACTATCAATGACACTTTCGCCATCGCAAAAATGGATAACAATATTGCAGGTGTCAGAATTACATCTCAGGCAGGCCCAGTATGGACAGATTTTCGGGGTAATGCCGTCATTCCATCAATCCAGCCGTGGCGAACATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCC  >  NC_000913/3364423‑3364898

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: