Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,252,464 T→C I54T (ATC→ACC)  yhaH → putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,252,4640TC96.4% 182.3 / NA 55I54T (ATC→ACC) yhaHputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (27/26);  minor base A (1/1);  total (28/27)
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.78e-01

GAACGAATATTTACTGCCATTTCATTCACGTCTATTCTTAATTTGCTGCTTTAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTT  >  NC_000913/3252219‑3252708
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
gAACGAATATTTACTGCCATTTCATTCACGTCTATTCTTAATTTGCTGCTTTAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTg                                                                                                                                                                                                                                                   >  1:259480/1‑250 (MQ=255)
      atatTTACTGCCATTTCATTCACGTCTATTCTTAATTTGCTGCTTTAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACa                                                                                                                                                                                                                                             >  2:137764/1‑250 (MQ=255)
        atTTACTTCCATTTCATTCACGTCTATTCTTAATTTGCTGCTTTAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAAt                                                                                                                                                                                                                                           <  1:298166/250‑1 (MQ=255)
        atTTACTGCCATTTCATTCACGTCTATTCTTAATTTGATGCTTTAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGTGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTCCGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACACTGACGACAAt                                                                                                                                                                                                                                           <  1:53583/250‑1 (MQ=255)
                     tcattcaCGTCTATTCTTAATTTGCTGCTTTAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTgg                                                                                                                                                                                                                             >  1:121919/1‑251 (MQ=255)
                     tcattcaCGTCTATTCTTAATTTGCTGCTTTAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTg                                                                                                                                                                                                                              <  2:357319/250‑1 (MQ=255)
                            cGTCTATTCTTAATTTGCTGCTTTAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGttttt                                                                                                                                                                                                                      <  1:156223/251‑1 (MQ=255)
                                tATTCTTAATTTGCTGCTTTAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCg                                                                                                                                                                                                                  >  1:446484/1‑251 (MQ=255)
                                tATTCTTAATTTGCTGCTTTAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCg                                                                                                                                                                                                                  >  1:149081/1‑251 (MQ=255)
                                    cTTAATTTGCTGCTTTAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCCTCTTTACGTTCGCGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCgtggt                                                                                                                                                                                                              >  2:32051/1‑251 (MQ=255)
                                         tttGCTGCTTTAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTAc                                                                                                                                                                                                                                       >  2:49709/1‑221 (MQ=255)
                                         tttGCTGCTTTAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTAc                                                                                                                                                                                                                                       <  1:49709/221‑1 (MQ=255)
                                         tttGCTGCTTTAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCg                                                                                                                                                                                                         >  1:383849/1‑251 (MQ=255)
                                          ttGCTGCTTTAAACACACGTGATTGAGAAAATTTGAGGCAAAACTGGACTGGTATCTGATAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGTTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTTGATAAAATGTTTGGCATGCAACGCGCTGGTGGCGAAGGCCCCCTGGAGACAATCTCCGGTATTCTGTTGTTTTTAccaatgtgggcgg                                                                                                                                                                                                        >  2:71752/1‑240 (MQ=255)
                                                ctTTAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGt                                                                                                                                                                                                  >  2:273057/1‑251 (MQ=255)
                                                  ttAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGAACTATAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGAAATGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGGTAAAATGTTAGGCTGGCAACGCGCTGGGGGCCAAGGCACCCTGACGCCAATCTACGGTATTCTGGTGTTTTTACCGTGGCGGGCGGTTTAGTTc                                                                                                                                                                                                >  2:490942/1‑251 (MQ=255)
                                                   tAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTAAGTTcc                                                                                                                                                                                               >  2:447476/1‑251 (MQ=255)
                                                               aTTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTAAGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCAACCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACg                                                                                                                                                                                   >  1:74707/1‑251 (MQ=255)
                                                                                      tGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGGGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGCCCGCTCGGCGtggtgg                                                                                                                                                            >  2:294064/1‑251 (MQ=255)
                                                                                            ggTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCAc                                                                                                                                                        <  2:149081/249‑1 (MQ=255)
                                                                                            ggTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACt                                                                                                                                                       <  1:107477/250‑1 (MQ=255)
                                                                                             gTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGc                                                                                                                                                     <  2:259480/251‑1 (MQ=255)
                                                                                                 cTGAAAGTACTAAACAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTAtt                                                                                                                                                 >  2:50220/1‑251 (MQ=255)
                                                                                                   gAAAGTACTAAAAAATTATGTCGGTGTCCGGGGTCGTGCGCGGCGCAAGGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCGGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGTGCACTGCTATTc                                                                                                                                                <  2:293921/250‑1 (MQ=255)
                                                                                                    aaaGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGGGGGCACTGCTATTCtt                                                                                                                                              >  1:128497/1‑251 (MQ=255)
                                                                                                                   tATGTCGGCTACCGGCGTCGTGCGCGCCGCAATGCGCACTGGTTGTTTATTCTGGCCACCATTATCTTTACGCTCGTGCTGGGGCTGCTGGATAACCTGTTAGGCTGGCCACGCGCTGGTGGCGAAGGCACCCTGCAGACAATCTACGGTATTCTGGGGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGAACGACAACGACCGCTCGGCGTGGAGGGCACAACTATTCTTAATCCCGTTTACCgg                                                                                                                               >  2:182125/1‑251 (MQ=255)
                                                                                                                         ggTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCACCATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTAAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCggctggct                                                                                                                         >  2:216365/1‑251 (MQ=255)
                                                                                                                              ccGGGGTCGTGCGCGGCGCAAAGAGTACGGGATGTTTATTCTGGGCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCttatta                                                                                                                    <  1:87264/251‑1 (MQ=255)
                                                                                                                               cGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTtattat                                                                                                                   <  1:50220/251‑1 (MQ=255)
                                                                                                                                 gggTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCt                                                                                                                                                                                                                               <  2:401073/141‑1 (MQ=255)
                                                                                                                                 gggTCGCGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCt                                                                                                                                                                                                                               >  1:401073/1‑141 (MQ=255)
                                                                                                                                  gggcGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCAGCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCCGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCGGGCTTATTatcat                                                                                                                <  2:211890/248‑1 (MQ=255)
                                                                                                                                     cGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGt                                                                                                             <  1:17448/251‑1 (MQ=255)
                                                                                                                                      gTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTTTTATCATCgcc                                                                                                            >  2:482167/1‑249 (MQ=255)
                                                                                                                                           cggcgCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTa                                                                                                        >  2:174422/1‑250 (MQ=255)
                                                                                                                                                  aaGAGTACTGGAGGTGTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTAGCGGCTGGCTTATTATCATCGTCTTTAACTGCCAg                                                                                                <  2:125626/251‑1 (MQ=255)
                                                                                                                                                               gTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCa                                                                                    <  1:82208/250‑1 (MQ=255)
                                                                                                                                                                     tCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGaa                                                                             <  2:446484/250‑1 (MQ=255)
                                                                                                                                                                      cTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTACTCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGaaa                                                                            >  1:396969/1‑250 (MQ=255)
                                                                                                                                                                                  atcTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTTTTTAACTGCCTGGCGGGTACGCCCGGCAGAAAACGCTTGGGAcc                                                               >  1:405813/1‑251 (MQ=255)
                                                                                                                                                                                  atcTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGGAAACCGCTTTGGAcc                                                               >  1:139700/1‑251 (MQ=255)
                                                                                                                                                                                      ttACGTTCGTGCTGGGGCTGCTTGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGGGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACGGCTATTCTTAATCCCGTTTACCGGCTGGCTTATTATCATCGTCTTGAACTGCCAGGCGGGTACGCCAGGCGAAAACCGCTTTGGAcc                                                               <  1:182125/247‑1 (MQ=255)
                                                                                                                                                                                       tACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGAt                                                           <  1:137764/250‑1 (MQ=255)
                                                                                                                                                                                              gTGCTGGGGCTGCTGGAAAAAATGGTAGGCTGGCAACGCGGTGGTGGCGAAGGCACCCTGACGACAATCTCCGGTATTCTGGTGTTTGGACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTAGCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGAAAGCCGCTTTTGACCATATCCAAAACt                                                   <  2:432001/251‑1 (MQ=255)
                                                                                                                                                                                                gcgggggCTGCTGGATAAAATGTTGGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTGCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGGGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTgg                                                 <  2:299263/248‑1 (MQ=255)
                                                                                                                                                                                                                aaaaTGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATc                                 >  1:101119/1‑251 (MQ=255)
                                                                                                                                                                                                                       tAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAAATGGAAACTTAAAAAAAATCTGccccc                          >  2:201760/1‑247 (MQ=255)
                                                                                                                                                                                                                        aGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGGAAACCGCTTTGGCACAGATCCAAAACCGGAACCTTAAAAAAAATTTGCTCGcc                         >  1:352153/1‑251 (MQ=255)
                                                                                                                                                                                                                              gCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGAAAAACGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATgg                   >  1:124490/1‑250 (MQ=255)
                                                                                                                                                                                                                                     gcTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTc            <  1:401093/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gcTGGTGGCGAAGGCACCCGGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTc            <  2:383849/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ggCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCg       <  1:174422/250‑1 (MQ=255)
                                                                                                                                                                                                                                              gAAGGCACCCTGACGACAATCTCCGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCtt   <  2:121919/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCtt   <  1:294064/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aaGGCAACCTGACGACAATCGACGGTGTTCTGGCGTGTTTAGCGGTGTGGGCGGTTCAGTTCCGCCTCCTGCACGACACCGACCGCTCGGCCTGTTGGGCACTGCTATTCTTAATCCCGTCTATCTGCTAGCTTAATATCATCGTCTTTAACTTTCAGGCGGGTACTCCAGGCGAAAACCGATTTGGGCCAGATCCAAAACGCGAACCTTAAAACAAATCTGCTCGCCGGATGTCTTACTCATCCGGCttc  <  2:15582/251‑2 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
GAACGAATATTTACTGCCATTTCATTCACGTCTATTCTTAATTTGCTGCTTTAAAAACAAGTGATTGAGCAAAATTGAGGCAAAAATGGACTGGTATCTGAAAGTACTAAAAAATTATGTCGGTTTCCGGGGTCGTGCGCGGCGCAAAGAGTACTGGATGTTTATTCTGGTCAACATCATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAGGCGGGTACGCCAGGCGAAAACCGCTTTGGACCAGATCCAAAACTGGAACCTTAAAACAAATCTGCTCGCCGGATGTCTGACTCATCCGGCTTT  >  NC_000913/3252219‑3252708

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