Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 537,024 C→T P147S (CCG→TCG)  glxR → tartronate semialdehyde reductase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913537,0240CT97.6% 130.3 / NA 42P147S (CCG→TCG) glxRtartronate semialdehyde reductase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/21);  minor base G (0/1);  total (20/22)
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.99e-01

GACACACCTCAGGTTGAAGAAGTTCTGTTCGGTGAAAATGGTTGTACCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAACCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCGGTACGTGTGCGCCAGGCGCTGATGGGCGGCTTTGCTTCCTCACGTATTCTGGAAGTTCATGGCGAGCGTA  >  NC_000913/536778‑537246
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               
gACACCCCTCAGGGTGAAGAAGTTCTGTTCGGTGAAAATGGTTGTACCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATGGAAACTAAGCGGTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGTGGTGAAATCGGTGCGCGTTAAGTGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCt                                                                                                                                                                                                                             <  2:242967/251‑1 (MQ=255)
    cacCTCAGGTTGAAGAAGTTCTGTTCGGTGAAAATGGTTGTACCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGttt                                                                                                                                                                                                                        >  1:423562/1‑251 (MQ=255)
                 agaagTACTGTCCGGTGCAGATGATTGTACCAACGCCTCGCTGAAGGGCAAAACCATTGTTGATATCAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGCTGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGCCGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGt                                                                                                                                                                                                            >  2:89058/1‑250 (MQ=255)
                         tGTTCGGTGAAAATGGTTGTACCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGACGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATc                                                                                                                                                                                                   >  1:500562/1‑251 (MQ=255)
                          gTTCGGTGAAAATGGTTGTACCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATTGACACTAAGCGTTTCGCTAGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATc                                                                                                                                                                                                   >  1:360610/1‑250 (MQ=255)
                              ggTGAAAATGGGTGTACCAAAGCCTCGCTGAAGGGCAAAACCATGGTTGATATGAGCTCCATTTCCCCGATTGAAACGAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCt                                                                                                                                                                                              <  2:337011/251‑1 (MQ=255)
                                tGAAAATGGTTGTACCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCg                                                                                                                                                                                            >  1:99468/1‑251 (MQ=255)
                                   aaaTGGTTGTACCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGTTATGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGAGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATAAAGCGGTATTTGAATGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTgg                                                                                                                                                                                         >  1:7564/1‑251 (MQ=255)
                                       ggTTGTCCCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCTTTTCCCCGATTGAAACTAAGCGTTGCGCTCGTCAGGTTAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGACGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCTCCCTCGTGGGCg                                                                                                                                                                                      <  1:118081/250‑1 (MQ=255)
                                               cAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCg                                                                                                                                                                             >  1:212700/1‑251 (MQ=255)
                                                        gCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCaaa                                                                                                                                                                    <  1:516738/251‑1 (MQ=255)
                                                        gCTGAAGGGCAAAACCATTGTTGATATGAGCGCTATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCaaa                                                                                                                                                                    <  1:285331/251‑1 (MQ=255)
                                                                 cAAAACCATTGCTCATATGAGATCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGACGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGAGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATAACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCaa                                                                                                                                                             >  2:77269/1‑249 (MQ=255)
                                                                          tgttgATATGAGCTCCATGTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAAGCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAAt                                                                                                                                                   <  1:86225/250‑1 (MQ=255)
                                                                                 atGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTAt                                                                                                                                           >  1:434760/1‑251 (MQ=255)
                                                                                     gCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAATCCTGCAAAGTGGCAAATCAGATTATCGTg                                                                                                                                       <  1:122191/251‑1 (MQ=255)
                                                                                        ccATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGgc                                                                                                                                     <  1:488905/250‑1 (MQ=255)
                                                                                        ccATTTCCCCGATTGAAACGAAGCGTTTCGCGCGTCAGGTGAATGAACTGGGCGGCGATTAGCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGGTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACGGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTgg                                                                                                                                      <  1:534888/249‑1 (MQ=255)
                                                                                            ttCCCCGATTGCAACTAAGCGTTTCGCTCGTCCGGTGCATGAATTGGGCGGAGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGGCAAAACTGCAAAGTGGCAAATCAGATTATCGGGGCGCTcc                                                                                                                                >  2:316004/1‑250 (MQ=255)
                                                                                                   aTTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGa                                                                                                                         <  2:283480/251‑1 (MQ=255)
                                                                                                       aaaCTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCg                                                                                                                     <  1:395911/251‑1 (MQ=255)
                                                                                                       aaaCTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCg                                                                                                                     <  2:212700/251‑1 (MQ=255)
                                                                                                       aaaCTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACATGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAg                                                                                                                       >  2:437634/1‑249 (MQ=255)
                                                                                                               cGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGa                                                                                                             <  2:423562/251‑1 (MQ=255)
                                                                                                                     gCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCt                                                                                                       >  1:399851/1‑251 (MQ=255)
                                                                                                                      cTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTc                                                                                                                                                                                             >  1:225519/1‑164 (MQ=255)
                                                                                                                      cTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAAGCGCTGGTTGAACTGCTCGGTAAAAATATCGCCCTc                                                                                                                                                                                             <  2:225519/164‑1 (MQ=255)
                                                                                                                                    gaaCTGGGCGGCGATTCTCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAAGGACGTTGTCGATTATGGTTAGCGGAGATGAAGCGGCATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCaa                                                                                         >  2:253045/1‑250 (MQ=255)
                                                                                                                                       cTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGc                                                                                     >  2:436733/1‑251 (MQ=255)
                                                                                                                                         gggCGGCGATTATCTCGAGGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCg                                                                                   <  1:379506/251‑1 (MQ=255)
                                                                                                                                             ggcgATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGAGTTAAATAGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAATCGGTTTCTGAAGCCCTGCTATTTGGTTCAAAAGCCGGTGc                                                                               >  2:478358/1‑251 (MQ=255)
                                                                                                                                                    aTCTCGATGCGCCAGTCTACGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCAACCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCg                                                                        >  1:527355/1‑251 (MQ=255)
                                                                                                                                                          aTGCGCCCGTCTCCGCCGGTGAAATCGGTTCGCGTCAACGGCCGTTGTCGATTTTGGTTGGCGGTGCTGAAGCGGTATTTGGACGTGTTCAATCGCTGTTTGAACTGCTCGGTCAAAATATCACCCTCGTGGGCGGTCCCGGCGATGTTCAAACCTGCACAGTGGCACATCAGATTATCGTGGCGCTCACTATTGTAGCGGTTTCTGCAGCCCTGCTATTTGCTTCAAAAGCCGGTGCTGACCCGGTCCgt                                                                  >  2:241840/1‑251 (MQ=255)
                                                                                                                                                            gcgcCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAATCCGTTGCGGACCCGGTACgtgt                                                                <  2:99468/251‑1 (MQ=255)
                                                                                                                                                               ccAGTCTCCGGCGGTGAAATCGGTGCTCGTGAAGTGACGTTTTCGATTATGGGTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCGCGTGGGCGGTAACGGCGATGGTCAACCCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCGGTACGTGTgcg                                                             <  2:423949/251‑1 (MQ=255)
                                                                                                                                                                             ggAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCGGTACGTGTGCGCCAGGCGCTGATgg                                               <  1:322983/250‑1 (MQ=255)
                                                                                                                                                                              gAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGGAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAATACCGGTGCGGACACGGGACGTGTGCGCCAGGCGCGGATgtg                                              >  2:308355/1‑249 (MQ=255)
                                                                                                                                                                                aaTCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCGGTACGTGTGCGCCAGGCGCTGATGggcg                                            >  2:516257/1‑251 (MQ=255)
                                                                                                                                                                                                          cGATGATGGTTGGCGGTGATGAAGCGGTATGTGAACGTGTTAAATCGCGTGTTGAACTGCTCGGTAAAAATATCACCCACGGGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCGCAATAGTGAAGCGGTTTCTGAAGCCCTGCTAGTTGTTTCAAAAGCCGGTGCGGACCCGGTACGTGGGCGACAGGCGCTGATGGGCGGCTTGGCTTCCTCACGTATTCTGGaa                  <  2:94135/251‑1 (MQ=255)
                                                                                                                                                                                                                     ggCGGTGATGAAGCGGTATTTGAACGGGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCTCCCTCGTGGGCGTTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTTAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCGGTACGTGTGCGTCAGGCGCTGATGGTCGGATTTGCTTCCTCACGTATTCTGGAAGTTCATGGCGa       <  2:66598/251‑1 (MQ=255)
                                                                                                                                                                                                                         gTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCGGTACGTGTGCGCCAGGCGCTGATGGGCGGCTTTGCTTCCTCACGTATTCTGGAAGTTCATGGCGAGCGt   <  1:14736/251‑1 (MQ=255)
                                                                                                                                                                                                                           gatgaAGCGGTATTTGAACGTGTTAAATCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCGGTACGTGTGCGCCAGGCGCTGATGGGCGGCTTTGCTTCCTCACGTATTCTGGAAGTTCATGGCGAGCGTa  <  2:115304/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               
GACACACCTCAGGTTGAAGAAGTTCTGTTCGGTGAAAATGGTTGTACCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAACCGCTGTTTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCGGTACGTGTGCGCCAGGCGCTGATGGGCGGCTTTGCTTCCTCACGTATTCTGGAAGTTCATGGCGAGCGTA  >  NC_000913/536778‑537246

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: