Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 46,448 C→A I214I (ATC→ATA yaaU → putative transporter YaaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091346,4480CA96.2% 151.6 / ‑3.8 53I214I (ATC→ATAyaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (21/29);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGCGACCTGTGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTG  >  NC_000913/46216‑46693
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
cccGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCg                                                                                                                                                                                                                                     <  2:12692/251‑1 (MQ=255)
cccGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGAGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGACATGATGATAAAACTGTTTGGAGAACCg                                                                                                                                                                                                                                     >  2:469152/1‑251 (MQ=255)
         gggCGTTTTCCATCTGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGTTCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGATATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAACATGTTTGGCGAACAGGAGGCtttg                                                                                                                                                                                                                            >  2:542839/1‑250 (MQ=255)
              ttttCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTCGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAACGAGGGCGAAGAGATGATGATAAAACTGTTTGGCGAACGGGTGGCTTTCgttg                                                                                                                                                                                                                        >  2:218096/1‑247 (MQ=255)
                    aTCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGagc                                                                                                                                                                                                                  >  1:357437/1‑250 (MQ=255)
                       aGCTTTATCGCCGCCATGTGGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCACTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGAGTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGcc                                                                                                                                                                                                              <  1:652487/251‑1 (MQ=255)
                                        gtgGTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGtttt                                                                                                                                                                                              <  1:461956/250‑1 (MQ=255)
                                          ggTATGTCGGCGCGACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTcg                                                                                                                                                                                            <  1:504944/250‑1 (MQ=255)
                                                gcGGCGCGGCCTGCGCCGATCTTGTCGGTTACTGGCTTTAGGATTGGGAAGGCGGCTGGCGATGGAGGCTGGGTATCGTGGCGATCCCCTGTTTGTTGATTTGGATTGGTCGATTCGAACTGCCTGAATCGCCCCGCTGGTTATTACGAATAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTGTGGCGATCCTGTGGCTTTTGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTg                                                                                                                                                                                     <  2:408958/250‑1 (MQ=255)
                                                  ggCGCGACCTGCGCCGATCTGGTCGGCTACTGGTTGTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATGCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGTGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGtt                                                                                                                                                                                   <  1:284279/251‑1 (MQ=255)
                                                      cgACCTGCGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTaa                                                                                                                                                                                >  2:602993/1‑250 (MQ=255)
                                                            gcgCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATcgccg                                                                                                                                                                          <  1:469152/248‑1 (MQ=255)
                                                                      gttCGGGTACTGGGTTTATGATGGGGAAGGCGGCTGGCGGTGGATGCTGGGTAGCGCGGCGATCCCCGGTTTGTTGATTTGGATTGGGCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTTGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCtt                                                                                                                                                               <  2:415308/249‑1 (MQ=255)
                                                                              attGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTGGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCt                                                                                                                                                        <  1:50130/248‑1 (MQ=255)
                                                                                     ttATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCtgtttgtt                                                                                                                                                 <  1:627187/250‑1 (MQ=255)
                                                                                          atgTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGTGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTgccgcc                                                                                                                                           <  1:265377/251‑1 (MQ=255)
                                                                                                      gcTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGc                                                                                                                               >  1:372452/1‑251 (MQ=255)
                                                                                                           cgcTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGt                                                                                                                          <  1:602993/251‑1 (MQ=255)
                                                                                                              tGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGAt                                                                                                                       >  1:484239/1‑251 (MQ=255)
                                                                                                                      ggttAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCAATTCGAACTGCCGGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCGTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAAt                                                                                                                 <  1:610323/246‑1 (MQ=255)
                                                                                                                      gggTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTAATCCCAATGt                                                                                                               >  1:218386/1‑251 (MQ=255)
                                                                                                                       ggTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGtt                                                                                                              >  1:644146/1‑251 (MQ=255)
                                                                                                                              cggctATCCCCTGTTTTTTGATTTTGATTGGTCGATTCGAACTGCCTGAACCTCCCCGCTGGTTATTCCGCAAAGGGCGAGTAAAAGAGTGCGAAGATATGATGATAAAACTGTTTGGCGAACCGGTGTCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCAtt                                                                                                       <  2:328457/251‑1 (MQ=255)
                                                                                                                                 cgATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTac                                                                                                    >  1:440414/1‑251 (MQ=255)
                                                                                                                                  gATCCACTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGCGTGCGAAGAGATGATGATCAAACTGTTTGGCGACCCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGGGATCCCAATGTTCTCCATTTacc                                                                                                   >  2:618116/1‑250 (MQ=255)
                                                                                                                                     ccccTGTTTTTGGATTTTGATTGGTCGAGTCGAACTGCATGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACAGTTTGGCGAACCGGGTGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCt                                                                                                <  2:357437/251‑1 (MQ=255)
                                                                                                                                       ccTGTTTGTTGCTTTTGATTGGGCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGACCCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCGGGACCTGCCAGGTGACCCCAATGTTCGCCATTTACACCttt                                                                                              <  1:484827/251‑1 (MQ=255)
                                                                                                                                       ccTGTTTGATGATTTTGATTGGTCGATTAGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTACGAAGAGAAGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGCAGAGCAGCAGCAGCAAACACGTTTTAGAGCTCTGTTTAATCTCCGCCATTTTCCTTTTGTTCGGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCACCATTTACCCCttt                                                                                              >  2:569073/1‑251 (MQ=255)
                                                                                                                                        cTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTg                                                                                             >  1:18799/1‑251 (MQ=255)
                                                                                                                                           tttgttGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCAATTTACACCTTTGGcc                                                                                          >  1:25942/1‑251 (MQ=255)
                                                                                                                                                  attttgattGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAAt                                                                                   <  1:630090/251‑1 (MQ=255)
                                                                                                                                                     ttgattGGTCGATTCGACCTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATTATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAATCCCGTTTTCGCGATCTTTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCg                                                                                 <  2:79230/250‑1 (MQ=255)
                                                                                                                                                       gattGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGtt                                                                               <  1:77043/250‑1 (MQ=255)
                                                                                                                                                          tGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGt                                                                            <  1:542839/250‑1 (MQ=255)
                                                                                                                                                          aggtagttgcaaattGCCTGAATTTCCCCTCTGGTTATTACGCAAAGGGCAAGTAAAATAGTGCGAAGGGATGATGAGAAAACTGTTTGGCGAACCGGTGGCTTCCGATGAAGAGCAGCCGCCGCAAACCCGTTTTCGCGATCTGTTTAAGCGCCGCCATTTTCCTTTTTTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTgg                                                                             <  2:603062/235‑1 (MQ=255)
                                                                                                                                                           ggTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGTTCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGttt                                                                          <  1:610630/251‑1 (MQ=255)
                                                                                                                                                             tcgtttCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATTAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTttgt                                                                        >  1:502694/5‑251 (MQ=255)
                                                                                                                                                               gATTCGAACTGCCGGAATTTCCCCGCGGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTGCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTttgttg                                                                      <  2:372452/251‑1 (MQ=255)
                                                                                                                                                                    gAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGAtt                                                                 >  1:560773/1‑251 (MQ=255)
                                                                                                                                                                     aaCTGCCTGAATCTCCCCGCTGGTTATTACGCGAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTg                                                                >  1:171072/1‑251 (MQ=255)
                                                                                                                                                                             gAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTgg                                                        >  1:90185/1‑251 (MQ=255)
                                                                                                                                                                                        gtgggtattaCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACGGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCgg                                             <  2:198851/246‑1 (MQ=255)
                                                                                                                                                                                        cTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGGTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTTCCAGGTGATCCCAATGTTCGCCATTTACGCCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCgg                                             >  1:545222/1‑251 (MQ=255)
                                                                                                                                                                                          ggTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGcg                                           <  1:27274/251‑2 (MQ=255)
                                                                                                                                                                                              attaCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTTGGTGTTGGCAAAAACGCGGcgctg                                        <  1:243356/250‑5 (MQ=255)
                                                                                                                                                                                                     aaggggCGAGTAAAAGAGGGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGaacgt                                <  2:484239/248‑4 (MQ=255)
                                                                                                                                                                                                     aaaGGGCGAGTAAAAGAGTGCGAAGAGATCATGATAAAACTGTTTGGCGAACAGGAGGCATACGATGAAGAGCAGACGCAGCAAGCCCGTGTTCGGGATCTGTTTAATCGCCGCCCTTCACCTGTTGTTCTGTTTGTAGCCGCAATCTGGCCATGTCAGGTGATCCCATAGTTCGCAATTTACACTTGTGGCCCGCGACTCGTCTGTATGTTGGGCTTGGGTGTTGGGAACAACGCggtgcgggggaacgg                                >  2:23149/1‑238 (MQ=255)
                                                                                                                                                                                                          gctaGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTCCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGa                           <  1:224068/248‑1 (MQ=255)
                                                                                                                                                                                                                         cGAAGAGATGATGATAAAACTGTTTGGCGAACCGTTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCACTGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCACAAACGCGGCGCTGGGGAACGTGATGATTAGCCTGTTCttt             >  1:196033/1‑250 (MQ=255)
                                                                                                                                                                                                                         cGAAAAGATGATGAAAAAACCGTTTTGCGAACCGGTGGCTTTCGATGAAAAGCGGCAGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCCTTTTCCTTTTGTTCTGTTTGTTGACGCCATCTGGACCTGCCATGTGATCCCAATGTTCTCCATTTACACCTTTGGCCCGCAAATCGTTAGTTTGTTGGGATTGAGTGTTGGCAAAAACGCGGCGCTGGGGTACGTGGTGAGTAGCCTGTTGttt             >  2:632653/1‑250 (MQ=255)
                                                                                                                                                                                                                          gAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAATCGCGGCGCTGGGGCACGTGGTGATTAGCCTGTTCTTTAt           >  2:341949/1‑251 (MQ=255)
                                                                                                                                                                                                                                atgatgatAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCAATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGACTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCgg     <  1:4070/251‑1 (MQ=255)
                                                                                                                                                                                                                                   atgatAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTg  >  1:90872/1‑251 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
CCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGCGACCTGTGCCGATCTGGTCGGCTACTGGCTTTATGATGTGGAAGGCGGCTGGCGCTGGATGCTGGGTAGCGCGGCGATCCCCTGTTTGTTGATTTTGATTGGTCGATTCGAACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTG  >  NC_000913/46216‑46693

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