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

GACACACCTCAGGTTGAAGAAGTTCTGTTCGGTGAAAATGGTTGTACCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAACCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCGGTACGTGTGCGCCAGGCGCTGATGGGCGGCTTTGCTTCCTCACGTATTCTGGAAGTTCATGGCGAGCGTA  >  NC_000913/536778‑537246
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               
gACACCCCTCAGGGTGAAGAAGTTCTGTTCGGTGAAAATGGTTGTACCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATGGAAACTAAGCGGTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGTGGTGAAATCGGTGCGCGTTAAGTGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCt                                                                                                                                                                                                                             <  2:242967/251‑1 (MQ=255)
    cacCTCAGGTTGAAGAAGTTCTGTTCGGTGAAAATGGTTGTACCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGt‑tt                                                                                                                                                                                                                        >  1:423562/1‑251 (MQ=255)
                 agaagTACTGTCCGGTGCAGATGATTGTACCAACGCCTCGCTGAAGGGCAAAACCATTGTTGATATCAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGCTGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGCCGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGt                                                                                                                                                                                                            >  2:89058/1‑250 (MQ=255)
                         tGTTCGGTGAAAATGGTTGTACCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGACGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATc                                                                                                                                                                                                   >  1:500562/1‑251 (MQ=255)
                          gTTCGGTGAAAATGGTTGTACCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATTGACACTAAGCGTTTCGCTAGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATc                                                                                                                                                                                                   >  1:360610/1‑250 (MQ=255)
                              ggTGAAAATGGGTGTACCAAAGCCTCGCTGAAGGGCAAAACCATGGTTGATATGAGCTCCATTTCCCCGATTGAAACGAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCt                                                                                                                                                                                              <  2:337011/251‑1 (MQ=255)
                                tGAAAATGGTTGTACCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCg                                                                                                                                                                                            >  1:99468/1‑251 (MQ=255)
                                   aaaTGGTTGTACCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGTTATGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGAGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATAAAGCGGTATTTGAATGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTgg                                                                                                                                                                                         >  1:7564/1‑251 (MQ=255)
                                       ggTTGTCCCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCTTTTCCCCGATTGAAACTAAGCGTTGCGCTCGTCAGGTTAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGACGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCTCCCTCGTGGGCg                                                                                                                                                                                      <  1:118081/250‑1 (MQ=255)
                                               cAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCg                                                                                                                                                                             >  1:212700/1‑251 (MQ=255)
                                                        gCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCaaa                                                                                                                                                                    <  1:516738/251‑1 (MQ=255)
                                                        gCTGAAGGGCAAAACCATTGTTGATATGAGCGCTATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCaaa                                                                                                                                                                    <  1:285331/251‑1 (MQ=255)
                                                                 cAAAACCATTGCTCATATGAGATCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGACGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGAGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATAACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCaa                                                                                                                                                             >  2:77269/1‑249 (MQ=255)
                                                                          tgttgATATGAGCTCCATGTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAAGCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAAt                                                                                                                                                   <  1:86225/250‑1 (MQ=255)
                                                                                 atGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTAt                                                                                                                                           >  1:434760/1‑251 (MQ=255)
                                                                                     gCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAATCCTGCAAAGTGGCAAATCAGATTATCGTg                                                                                                                                       <  1:122191/251‑1 (MQ=255)
                                                                                        ccATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGgc                                                                                                                                     <  1:488905/250‑1 (MQ=255)
                                                                                        ccATTTCCCCGATTGAAACGAAGCGTTTCGCGCGTCAGGTGAATGAACTGGGCGGCGATTAGCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGGTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACGGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTgg                                                                                                                                      <  1:534888/249‑1 (MQ=255)
                                                                                            ttCCCCGATTGCAACTAAGCGTTTCGCTCGTCCGGTGCATGAATTGGGCGGAGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGGCAAAACTGCAAAGTGGCAAATCAGATTATCGGGGCGCTcc                                                                                                                                >  2:316004/1‑250 (MQ=255)
                                                                                                   aTTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGa                                                                                                                         <  2:283480/251‑1 (MQ=255)
                                                                                                       aaaCTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCg                                                                                                                     <  1:395911/251‑1 (MQ=255)
                                                                                                       aaaCTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCg                                                                                                                     <  2:212700/251‑1 (MQ=255)
                                                                                                       aaaCTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACATGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAg                                                                                                                       >  2:437634/1‑249 (MQ=255)
                                                                                                               cGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGa                                                                                                             <  2:423562/251‑1 (MQ=255)
                                                                                                                     gCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCt                                                                                                       >  1:399851/1‑251 (MQ=255)
                                                                                                                      cTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTc                                                                                                                                                                                             >  1:225519/1‑164 (MQ=255)
                                                                                                                      cTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAAGCGCTGG‑TTGAACTGCTCGGTAAAAATATCGCCCTc                                                                                                                                                                                             <  2:225519/164‑1 (MQ=255)
                                                                                                                                    gaaCTGGGCGGCGATTCTCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAAGGACGTTGTCGATTATGGTTAGCGGAGATGAAGCGGCATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCaa                                                                                         >  2:253045/1‑250 (MQ=255)
                                                                                                                                       cTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGc                                                                                     >  2:436733/1‑251 (MQ=255)
                                                                                                                                         gggCGGCGATTATCTCGAGGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCg                                                                                   <  1:379506/251‑1 (MQ=255)
                                                                                                                                             ggcgATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGAGTTAAATAGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAATCGGTTTCTGAAGCCCTGCTATTTGGTTCAAAAGCCGGTGc                                                                               >  2:478358/1‑251 (MQ=255)
                                                                                                                                                    aTCTCGATGCGCCAGTCTACGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCAACCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCg                                                                        >  1:527355/1‑251 (MQ=255)
                                                                                                                                                          aTGCGCCCGTCTCCGCCGGTGAAATCGGTTCGCGTCAACGGCCGTTGTCGATTTTGGTTGGCGGTGCTGAAGCGGTATTTGGACGTGTTCAATCGCTGT‑TTGAACTGCTCGGTCAAAATATCACCCTCGTGGGCGGTCCCGGCGATGTTCAAACCTGCACAGTGGCACATCAGATTATCGTGGCGCTCACTATTGTAGCGGTTTCTGCAGCCCTGCTATTTGCTTCAAAAGCCGGTGCTGACCCGGTCCgt                                                                  >  2:241840/1‑251 (MQ=255)
                                                                                                                                                            gcgcCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAATCCGTTGCGGACCCGGTACgtgt                                                                <  2:99468/251‑1 (MQ=255)
                                                                                                                                                               ccAGTCTCCGGCGGTGAAATCGGTGCTCGTGAAGTGACGTTTTCGATTATGGGTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCGCGTGGGCGGTAACGGCGATGGTCAACCCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCGGTACGTGTgcg                                                             <  2:423949/251‑1 (MQ=255)
                                                                                                                                                                             ggAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCGGTACGTGTGCGCCAGGCGCTGATgg                                               <  1:322983/250‑1 (MQ=255)
                                                                                                                                                                              gAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGGAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAATACCGGTGCGGACACGGGACGTGTGCGCCAGGCGCGGATgtg                                              >  2:308355/1‑249 (MQ=255)
                                                                                                                                                                                aaTCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCGGTACGTGTGCGCCAGGCGCTGATGggcg                                            >  2:516257/1‑251 (MQ=255)
                                                                                                                                                                                                          cGATGATGGTTGGCGGTGATGAAGCGGTATGTGAACGTGTTAAATCGC‑GTGTTGAACTGCTCGGTAAAAATATCACCCACGGGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCGCAATAGTGAAGCGGTTTCTGAAGCCCTGCTAGTTGTTTCAAAAGCCGGTGCGGACCCGGTACGTGGGCGACAGGCGCTGATGGGCGGCTTGGCTTCCTCACGTATTCTGGaa                  <  2:94135/251‑1 (MQ=255)
                                                                                                                                                                                                                     ggCGGTGATGAAGCGGTATTTGAACGGGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCTCCCTCGTGGGCGTTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTTAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCGGTACGTGTGCGTCAGGCGCTGATGGTCGGATTTGCTTCCTCACGTATTCTGGAAGTTCATGGCGa       <  2:66598/251‑1 (MQ=255)
                                                                                                                                                                                                                         gTGATGAAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCGGTACGTGTGCGCCAGGCGCTGATGGGCGGCTTTGCTTCCTCACGTATTCTGGAAGTTCATGGCGAGCGt   <  1:14736/251‑1 (MQ=255)
                                                                                                                                                                                                                           gatgaAGCGGTATTTGAACGTGTTAAATCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCGGTACGTGTGCGCCAGGCGCTGATGGGCGGCTTTGCTTCCTCACGTATTCTGGAAGTTCATGGCGAGCGTa  <  2:115304/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               
GACACACCTCAGGTTGAAGAAGTTCTGTTCGGTGAAAATGGTTGTACCAAAGCCTCGCTGAAGGGCAAAACCATTGTTGATATGAGCTCCATTTCCCCGATTGAAACTAAGCGTTTCGCTCGTCAGGTGAATGAACTGGGCGGCGATTATCTCGATGCGCCAGTCTCCGGCGGTGAAATCGGTGCGCGTGAAGGGACGTTGTCGATTATGGTTGGCGGTGATGAAGCGGTATTTGAACGTGTTAAACCGCTGT‑TTGAACTGCTCGGTAAAAATATCACCCTCGTGGGCGGTAACGGCGATGGTCAAACCTGCAAAGTGGCAAATCAGATTATCGTGGCGCTCAATATTGAAGCGGTTTCTGAAGCCCTGCTATTTGCTTCAAAAGCCGGTGCGGACCCGGTACGTGTGCGCCAGGCGCTGATGGGCGGCTTTGCTTCCTCACGTATTCTGGAAGTTCATGGCGAGCGTA  >  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: ‑