Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,812,802 G→A A63T (GCC→ACC)  emrB → multidrug efflux pump membrane subunit EmrB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,812,8020GA97.9% 173.0 / ‑4.0 48A63T (GCC→ACC) emrBmultidrug efflux pump membrane subunit EmrB
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (33/14);  total (34/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

GATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATGCCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCCGCGTGATTCAGGGGATTGTCGCCGGGCCGTTGATCCCGCTTTCGCAAAGTCTATTGCTGAATAACTACCCGCCAGCCAAACGC  >  NC_000913/2812566‑2813023
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                             
gATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATAcc                                                                                                                                                                                                                             >  1:574945/1‑239 (MQ=255)
gATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATAcc                                                                                                                                                                                                                             >  1:296367/1‑239 (MQ=255)
gATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCa                                                                                                                                                                                                                            >  2:425284/1‑240 (MQ=255)
gATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCa                                                                                                                                                                                                                            >  2:136286/1‑240 (MQ=255)
          atCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATcc                                                                                                                                                                                                                   <  1:462364/239‑1 (MQ=255)
                   gCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCtt                                                                                                                                                                                                              >  2:721861/1‑235 (MQ=255)
                             gCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCtggctggc                                                                                                                                                                                                >  2:277446/1‑239 (MQ=255)
                                       agGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCg                                                                                                                                                                                     >  1:728777/1‑240 (MQ=255)
                                            cGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGaa                                                                                                                                                                                >  1:535478/1‑240 (MQ=255)
                                            cGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGa                                                                                                                                                                                 >  2:13746/1‑239 (MQ=255)
                                                 gATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGa                                                                                                                                                                            <  2:391339/239‑1 (MQ=255)
                                                         aGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTc                                                                                                                                                                    <  1:553007/239‑1 (MQ=255)
                                                           cAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCt                                                                                                                                                                  <  1:482613/239‑1 (MQ=255)
                                                             aaaaCCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCttt                                                                                                                                                                >  2:187484/1‑239 (MQ=255)
                                                             aaaaCCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCttt                                                                                                                                                                >  1:635757/1‑239 (MQ=255)
                                                                   gCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTccac                                                                                                                                                         >  1:636661/1‑240 (MQ=255)
                                                                                      gTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGcg                                                                                                                                      >  1:236814/1‑240 (MQ=255)
                                                                                         aTTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGcgtcg                                                                                                                                   <  1:13746/240‑1 (MQ=255)
                                                                                           tATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCg                                                                                                                                                    >  1:527728/1‑221 (MQ=255)
                                                                                             tGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTg                                                                                                                                              <  2:636651/225‑1 (MQ=255)
                                                                                                             aCTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCcag                                                                                                               <  1:277449/240‑1 (MQ=255)
                                                                                                                ggCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCcagcag                                                                                                            >  1:716937/1‑240 (MQ=255)
                                                                                                                 gCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCagcagc                                                                                                           <  1:362961/240‑1 (MQ=255)
                                                                                                                     cATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGa                                                                                                       <  1:289231/240‑1 (MQ=255)
                                                                                                                           tGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTg                                                                                                        >  2:265345/1‑233 (MQ=255)
                                                                                                                               ggTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGAt                                                                                             <  1:429824/240‑1 (MQ=255)
                                                                                                                                  gCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGtct                                                                                                                   >  1:162994/1‑215 (MQ=255)
                                                                                                                                         tCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCCCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTccacca                                                                                                                                                       <  2:490553/172‑1 (MQ=255)
                                                                                                                                           caccaTTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCcgcg                                                                                  >  1:74670/1‑239 (MQ=255)
                                                                                                                                                ttGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCACTCAGCCAGGGAAAGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCg                                                                                                                                                                                                          >  2:545118/1‑114 (MQ=255)
                                                                                                                                                  gcgaaCGTGGAGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGAAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATGCCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTACTTTGGTCCAGCATc                                                                                                                                                     >  2:312056/4‑165 (MQ=255)
                                                                                                                                                        gTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCCGCGt                                                                                 >  2:369294/1‑227 (MQ=255)
                                                                                                                                                           gCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCCGCGTGATTCAGGGGATTGt                                                                  >  1:95767/1‑239 (MQ=255)
                                                                                                                                                           gCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCCGCGTGATTCAGGGGATTGt                                                                  >  1:417210/1‑239 (MQ=255)
                                                                                                                                                           gCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCCGAGTGATTCAGGGGAt                                                                     >  2:235619/1‑236 (MQ=255)
                                                                                                                                                                         cGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCCGCGTGATTCAggg                                                                        >  1:372261/1‑219 (MQ=255)
                                                                                                                                                                                     gggCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCCGCGTGATTCAGGGGATTGTCGCCGGGCCGTTGATCCCGCTTTCGCa                                       <  1:540864/240‑1 (MQ=255)
                                                                                                                                                                                           aTCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCAt                                                                                                                                                      >  2:598899/1‑123 (MQ=255)
                                                                                                                                                                                             cGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCCGCGTGATTCAGGGGATTGTCGCCGGGCCGTTGATCCCGCTTTCGCAAAGTCTAt                               >  1:368441/1‑240 (MQ=255)
                                                                                                                                                                                                    ccAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGAtct                                                                                           <  1:721874/173‑1 (MQ=255)
                                                                                                                                                                                                     cAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCCGCGTGATTCAGGGGATTGTCGCCGGGCCGTTGATCCCGCTTTCGCaa                                      >  1:650181/1‑225 (MQ=255)
                                                                                                                                                                                                        gggaCGTGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCCGCGTGATTCAGGGGATTGTCGCCGGGCCGTTGATCCCGCTTTCGCAAAGTCTATTGCTGAATaa                     >  2:643445/4‑239 (MQ=255)
                                                                                                                                                                                                              tGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCCGCGTGATTCAGGGGATTGTCGCCGGGCCGTTGATCCCGCTTTCGCAAAGTCTATTGCTGAATAACTAccc               >  1:22661/1‑239 (MQ=255)
                                                                                                                                                                                                              tGGGTAATCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCCGCGTGATTCAGGGGATTGTCGCCGGGCCGTTGATCCCGCTTTCGCAAAGTCTATTGCTGAATAACTAccc               >  1:445138/1‑239 (MQ=255)
                                                                                                                                                                                                                    aTCACTTCTTTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCCGCGTGATTCAGGGGATTGTCGCCGGGCCGTTGATCCCGCTTTCGCAAAGTCTATTGCTGAATAACTACCCg              >  2:403278/1‑234 (MQ=255)
                                                                                                                                                                                                                         ttcttTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCCGCGTGATTCAGGGGATTGTCGCCGGGCCGTTGATCCCGCTTTCGCAAAGTCTAt                               <  1:100602/212‑1 (MQ=255)
                                                                                                                                                                                                                         ttcttTCGGGGTGGCGAATACCATCCCGATCCCGCTTACCGGCTGGCTGGCCAAACGCGTCg                                                                                                                                                                                     >  2:134455/1‑62 (MQ=255)
                                                                                                                                                                                                                           cttTCGGGGTGGCGAATACCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCCGCGTGATTCAGGGGATTGTCGCCGGGCCGTTGATCCCGCTTTCGCAAAGTCTATTGCTGAATAACTACCCGCCAGCCAAACGc  >  1:186339/1‑239 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                             
GATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATGCCATCTCGATCCCGCTTACCGGCTGGCTGGCAAAGCGCGTCGGGGAAGTGAAACTGTTCCTTTGGTCCACCATCGCCTTTGCTATTGCGTCGTGGGCGTGTGGTGTCTCCAGCAGCCTGAATATGCTGATCTTCTTCCGCGTGATTCAGGGGATTGTCGCCGGGCCGTTGATCCCGCTTTCGCAAAGTCTATTGCTGAATAACTACCCGCCAGCCAAACGC  >  NC_000913/2812566‑2813023

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