Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 536,173 T→C L155P (CTC→CCC)  hyi → hydroxypyruvate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913536,1730TC97.6% 141.0 / ‑6.1 42L155P (CTC→CCC) hyihydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (23/18);  total (24/18)
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.82e-01

ATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGT  >  NC_000913/535926‑536419
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
aTTGCCTGTATTCCTGGCCGTGCAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCAc                                                                                                                                                                                                                                                     >  1:206239/1‑251 (MQ=255)
aTTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCAc                                                                                                                                                                                                                                                     >  1:68723/1‑251 (MQ=255)
        tATTCATGGCCGTGAAGAAGAGGTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCAGGTGCGATGGGTAATAAAAAAAATAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGGAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTAGTACTGATTGAACCTATTAACAAGTATGATATTCCTGGTTTCCATACAACCGGAACTc                                                                                                                                                                                                                                             >  2:483272/1‑251 (MQ=255)
                   gTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTg                                                                                                                                                                                                                                  >  1:427462/1‑251 (MQ=255)
                                     agggCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAATCGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCTGCATGCGCTGAAATTGATTGATGATgtt                                                                                                                                                                                                                <  1:358830/249‑1 (MQ=255)
                                         cGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGCACTCGGCAGGCGCTGAAATTGCTTGATGATgttggtt                                                                                                                                                                                                            >  2:59308/1‑251 (MQ=255)
                                          gTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTACAGCAGTGAACAGATTCATGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgttggttg                                                                                                                                                                                                           <  1:221324/251‑1 (MQ=255)
                                            agcagcagcGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTtgct                                                                                                                                                                                                         >  1:402060/1‑251 (MQ=255)
                                                     gATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGa                                                                                                                                                                                                <  1:390718/251‑1 (MQ=255)
                                                         cGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAAt                                                                                                                                                                                            >  2:64542/1‑251 (MQ=255)
                                                         cGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGCCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAAt                                                                                                                                                                                            >  1:190422/1‑251 (MQ=255)
                                                               gCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGt                                                                                                                                                                                       <  1:231151/250‑1 (MQ=255)
                                                                   gTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGAc                                                                                                                                                                                  <  1:166049/251‑1 (MQ=255)
                                                                   gTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGAc                                                                                                                                                                                  <  2:489353/251‑1 (MQ=255)
                                                                      cgcTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACAtt                                                                                                                                                                               >  1:509461/1‑251 (MQ=255)
                                                                        cTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGATTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACAGGAAATCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAACTTCAGTATGACAttt                                                                                                                                                                              >  2:401843/1‑250 (MQ=255)
                                                                                taaAACAATTAACTTTCTGGTCGGTAAAACGCCGGCTGGTTTCATCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCGGATGCAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCGATCCCACCTGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTACAATTTGAAAATTCAGTATGACATTTAtgatttgc                                                                                                                                                                     >  2:12761/1‑244 (MQ=255)
                                                                                  aacaaaTTAACTGTATGGTCGGAAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGGAGACAACCTGCGTTATGCCGCGAATATGCTGATGAAAGCAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGCAATTGATTGGTGATGTTGGTTGCTGCAATTTGAACATTCCGTATGACATTTATCatattc                                                                                                                                                                     >  2:187230/4‑247 (MQ=255)
                                                                                                  ggTCTGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAAt                                                                                                                                                   <  1:411207/251‑1 (MQ=255)
                                                                                                                  gCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGAc                                                                                                                                   >  1:311271/1‑251 (MQ=255)
                                                                                                                                         tCACGCAACGCTTGTAGAAATCCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGt                                                                                                             <  2:485248/250‑1 (MQ=255)
                                                                                                                                                     tGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAAtt                                                                                                 <  1:355277/250‑1 (MQ=255)
                                                                                                                                                                 gCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAAGCCg                                                                                     <  2:359639/250‑1 (MQ=255)
                                                                                                                                                                   gTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCAt                                                                                  <  1:354585/251‑1 (MQ=255)
                                                                                                                                                                                 atGCTGATGAAAGAAGATATTTTATTACTTATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGTAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGAGGATGAAAGGTGAATTAACCAACAACATGACTCAGTGGGCGGATAAAATTGGTCACCTGCAAATAGCCGATAATCCGCATCGCGGAgaaccggg                                                                    >  2:355003/1‑251 (MQ=255)
                                                                                                                                                                                  tGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCgaaccggg                                                                    >  2:21917/1‑250 (MQ=255)
                                                                                                                                                                                       atgaAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCgaaccgggaaccgg                                                              <  2:442411/251‑1 (MQ=255)
                                                                                                                                                                                                 atatTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAaattaat                                                    >  1:140344/1‑251 (MQ=255)
                                                                                                                                                                                                 atatTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAaattaat                                                    <  1:174296/251‑1 (MQ=255)
                                                                                                                                                                                                  tatTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAaattaatt                                                   >  2:7666/1‑251 (MQ=255)
                                                                                                                                                                                                        attaCTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGAtt                                             >  1:107545/1‑251 (MQ=255)
                                                                                                                                                                                                        aataCTGATTGAACATATTAACCATTTTGATATTCCTGGTTTCAATCCCACCGGAACTATGCAGGCGCTGAAATTGATTGATGAAGTTGGTTGCTGCACTTTGAAAGTTCAGAATGACATTTATCATATGCAGCGGATGGAAGGAGAATTAAACAACAACATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGACCCGGAGAAATTCATCATGAtt                                             >  2:378203/3‑251 (MQ=255)
                                                                                                                                                                                                                  gAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTa                                    <  1:147552/250‑1 (MQ=255)
                                                                                                                                                                                                                       tATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTa                               >  2:131119/1‑250 (MQ=255)
                                                                                                                                                                                                                                  ttGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTc                    <  1:116814/250‑1 (MQ=255)
                                                                                                                                                                                                                                       atTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACTTTTATCATATGCAGCGGATGGAAGGTGACTTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACt               >  1:235402/1‑250 (MQ=255)
                                                                                                                                                                                                                                       atTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAAGTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGCTTATCTCTTTAAGGTAATCGAAAATTCTGACt               <  1:404525/250‑1 (MQ=255)
                                                                                                                                                                                                                                             ggTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAAc          <  1:209307/249‑1 (MQ=255)
                                                                                                                                                                                                                                             ggTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACgg        >  2:341063/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ttCCATCTCACAGGAAATCGCCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGACTTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCAGATAATCCGCATTGCGGCGCACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACCACagttg     >  2:149067/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ttCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACggtt      <  1:7666/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   cATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACggttgggt  >  1:455004/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
ATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGT  >  NC_000913/535926‑536419

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