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,1730TC98.1% 189.0 / NA 53L155P (CTC→CCC) hyihydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (29/23);  minor base A (1/0);  total (30/23)
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

ATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGG  >  NC_000913/535935‑536413
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
aTTCCTGGTCGTGAAGAGGAGGTCCAGGATGGCGTAGCAGCAGCGATTCGTTTTGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCg                                                                                                                                                                                                                                        >  2:206000/1‑251 (MQ=255)
aTTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCg                                                                                                                                                                                                                                        <  1:379109/251‑1 (MQ=255)
          gTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTg                                                                                                                                                                                                                              >  1:547004/1‑251 (MQ=255)
               gaagaGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCGCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGATCAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCAGTTTGATATTCCAGGTTTCCAACCCACCGGAACTCGGCAGGCGCTGAAatt                                                                                                                                                                                                                         <  2:120182/251‑1 (MQ=255)
                            aTGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCGATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgtt                                                                                                                                                                                                            <  1:295464/251‑1 (MQ=255)
                                cGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGGAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgttggtt                                                                                                                                                                                                        >  2:324849/1‑251 (MQ=255)
                                             aTTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAGAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGaa                                                                                                                                                                                           <  1:208012/251‑1 (MQ=255)
                                                 gTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAAtt                                                                                                                                                                                       <  1:513735/251‑1 (MQ=255)
                                                           tGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGAc                                                                                                                                                                              >  1:401844/1‑250 (MQ=255)
                                                              gcTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACAttt                                                                                                                                                                          <  2:184693/251‑1 (MQ=255)
                                                              gcTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACAttt                                                                                                                                                                          >  1:442084/1‑251 (MQ=255)
                                                                 ggggaataaAACAATTAACTGTCTGGCCGGCAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGTTTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGCAATTGATTAATGATGTTGGTTGCTGCAATTTGAAAATTCAGTCTGACATTTATc                                                                                                                                                                       >  2:474572/5‑251 (MQ=255)
                                                                   gTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGAc                                                                                                                                                                              >  2:286467/1‑242 (MQ=255)
                                                                   gTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGAc                                                                                                                                                                              <  1:286467/242‑1 (MQ=255)
                                                                     aataaAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCatat                                                                                                                                                                   <  1:259915/251‑1 (MQ=255)
                                                                                 atcggTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCGTGTAGAAAACCTGCGTTATGCAGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTAGTAACCATTTTGATATTCCTGGGTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAAGTCAGTATGACATTTATCATATGAAGCGGATgg                                                                                                                                                        <  2:618508/246‑1 (MQ=255)
                                                                                        gggTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGaa                                                                                                                                                <  1:327096/250‑1 (MQ=255)
                                                                                                aaaaCGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑Ca                                                                                                                                        <  2:442084/250‑1 (MQ=255)
                                                                                                         gCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGAc                                                                                                                              >  1:366565/1‑251 (MQ=255)
                                                                                                                      gTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGAt                                                                                                                 <  2:530683/251‑1 (MQ=255)
                                                                                                                                 cACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑c                                                                                                      >  2:96395/1‑250 (MQ=255)
                                                                                                                                  aCGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CAc                                                                                                    >  1:289747/1‑251 (MQ=255)
                                                                                                                                            tgaaGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTAAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTg                                                                                          >  1:292783/4‑251 (MQ=255)
                                                                                                                                                gAAAACCGGCGTTACGCCGCGGATATGATGACGCAAGAAGATCTTTTATTACTGATTGAACCTATTCATCCTTTTAATATTCCTGGGTTCCATCCCACCGGCACTCAGCAGGCGCTGCAATTGATTACTGCTGTTGGCTGCTGCAACTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGAGGGCTGATAAAATTGGA‑CACATGCAAATTGCCgc                                                                                      >  2:502954/1‑250 (MQ=255)
                                                                                                                                                   caaatgcgcaaTGCCGCACATATGCTGATGAAAGAAGATATTTTAGTACTGATTGATCCTATTAACTATTTTAATAGGCCCGATTTCCATCACACCGGAACTCGGCAGTCGTTGATATTGATAGTTGATAATGTTTGCTGTAATTGGAACATGTTGTCTGCCATTTATCATATGCAGCGGATGGAAGGTGAATTAACACAA‑ACCTTGACTACGTGGGCTGCTAAAATTGGT‑CACATGGAAAttcgcgcga                                                                                    >  2:6530/11‑243 (MQ=255)
                                                                                                                                                   aaCCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATa                                                                                    >  1:629138/1‑250 (MQ=255)
                                                                                                                                                          gTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCAt                                                                            <  1:531839/251‑1 (MQ=255)
                                                                                                                                                          gTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCAt                                                                            >  1:53736/1‑251 (MQ=255)
                                                                                                                                                          gTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCAt                                                                            <  1:7480/251‑1 (MQ=255)
                                                                                                                                                            tATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATc                                                                           <  2:513647/250‑1 (MQ=255)
                                                                                                                                                             aTGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATcgc                                                                         >  2:235586/1‑251 (MQ=255)
                                                                                                                                                                ccGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCgcggc                                                                      >  1:501034/1‑251 (MQ=255)
                                                                                                                                                                    gAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAat                                                                                                                                                                                                                          <  1:270491/101‑1 (MQ=255)
                                                                                                                                                                    gAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAat                                                                                                                                                                                                                          >  2:270491/1‑101 (MQ=255)
                                                                                                                                                                       tatGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGGTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCgaaccgg                                                               >  2:354407/1‑251 (MQ=255)
                                                                                                                                                                             gatgaAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCgaaccgggaaccg                                                         >  1:170505/1‑251 (MQ=255)
                                                                                                                                                                                gaAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCgaaccgggaaccgg                                                        >  1:396344/1‑249 (MQ=255)
                                                                                                                                                                                gaAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGgag                                                      <  2:53736/251‑1 (MQ=255)
                                                                                                                                                                                    gaagaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAaat                                                  >  1:363273/1‑251 (MQ=255)
                                                                                                                                                                                       gaTATTTTATGATTTATTGAACCAATTAACTATTTTGATATTCCTGGTTTCCTTCCCACCGGAAGACGGCAGGCTCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGCCATTTATCATATGCAGCGGATGGAAGGTGAATGAAC‑CCACACCATGACTCAGTGGGCCGATAAA‑CTGTTACACCTGCAAATTGCCGATAATCCGCATAGCGGCGAACGGGGAACCGCAGAattta                                                >  2:391390/1‑250 (MQ=255)
                                                                                                                                                                                        atatTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAaattaat                                              >  1:22925/1‑251 (MQ=255)
                                                                                                                                                                                          atTTTATTACTGATTGCACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTTAGTATGACCTTTATCATATGCAGTGAATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGACCCGGGAACCGGAGAaattaatt                                             >  1:290528/1‑250 (MQ=255)
                                                                                                                                                                                             ttagtaCTGATTGAACCTATTAACCATTGTGAGATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATg                                          <  2:6505/246‑1 (MQ=255)
                                                                                                                                                                                                  aCTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTAtc                                    >  1:255988/1‑251 (MQ=255)
                                                                                                                                                                                                   cTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTAtct                                   >  1:28048/1‑251 (MQ=255)
                                                                                                                                                                                                    tGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTAtctc                                  >  2:30521/1‑251 (MQ=255)
                                                                                                                                                                                                            ccTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGt                          <  2:255988/251‑1 (MQ=255)
                                                                                                                                                                                                                         ttGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAAGTAATTATGATTATCTCTTTAAGGTAATCGAAAATTc              <  2:452857/250‑1 (MQ=255)
                                                                                                                                                                                                                           gATATTCCTGGTTTCCATCCTACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACAGGGAACCGGAGAAATTAATTAATATTATCTCTTTAAGGTCATCGAAAATTCTg            >  2:431743/1‑250 (MQ=255)
                                                                                                                                                                                                                                tCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTAc       >  1:555602/1‑250 (MQ=255)
                                                                                                                                                                                                                                tCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTAc       <  1:206000/250‑1 (MQ=255)
                                                                                                                                                                                                                                    ggTTTCCATCCCACCGGAACTCTGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACgg  <  1:466232/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ggTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATGCTGACTACAAct   <  2:501034/250‑2 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
ATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAAC‑CAACACCATGACTCAGTGGGCTGATAAAATTGGT‑CACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGG  >  NC_000913/535935‑536413

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: ‑