Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 80,529 C→G G112G (GGG→GGC)  leuC ← 3‑isopropylmalate dehydratase subunit LeuC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091380,5290CG95.5% 205.1 / 3.0 68G112G (GGG→GGC) leuC3‑isopropylmalate dehydratase subunit LeuC
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base G (28/36);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

CCGGTTTTACCGATAATTGCCAGCACGATATCTTTTGCGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGAGGTCACTTCATGCACCAGG  >  NC_000913/80290‑80775
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
ccGGTGTTACCGATAATTGCCAGCACGATATCGTTTGCGGTCATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTGGCCAGTACGTGTTCAACTGCGGAAGTGCCGATACCAAAGGCCAGGGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCTCTGTTCCGGGCCCATTACGTg                                                                                                                                                                                                                                             <  2:678390/251‑1 (MQ=255)
   gTTTTACCGATAATTGCCAGCACGATATCTTTTGCGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGAc                                                                                                                                                                                                                                          >  2:232674/1‑251 (MQ=255)
         ccGATAATTGCCAGCACGATATCTTTTGCGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTAGCCGCAGACAATGGTCATCCCCGGCAAGGGGACGCCCTGTTCCGGGCCCATTACGTGGACGATccc                                                                                                                                                                                                                                    >  2:438663/1‑251 (MQ=255)
          cGATAATTTCCAGCACGATATCTTTTGCGGTAATGCCCGGCGCGGCTTTGGCCTGGACTGCAATTTTCATGGTTTTTGCGCGGCCCTGTTGCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATcctc                                                                                                                                                                                                                                   <  2:614527/251‑3 (MQ=255)
             tAATTGCCAGCACGATATCTTTTGCGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGCAGATCCCCTg                                                                                                                                                                                                                                 >  1:381087/1‑250 (MQ=255)
                  gCCAGCACGATATCTTTTGCGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGGTCCGGGCCCATTACGTGGACGATCCCCTGATACgg                                                                                                                                                                                                                           <  2:568606/251‑1 (MQ=255)
                  gCCAGCACGATATCTTTTGCGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGATCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCTAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATGCCCTGATACgg                                                                                                                                                                                                                           >  1:1895/1‑251 (MQ=255)
                                ttttGCGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCCAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGCGGATTCAGGTCAt                                                                                                                                                                                                             >  1:35924/1‑251 (MQ=255)
                                  ttGCGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATa                                                                                                                                                                                                            <  1:368515/250‑1 (MQ=255)
                                  ttGCGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATAc                                                                                                                                                                                                           >  1:330270/1‑251 (MQ=255)
                                     cGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGt                                                                                                                                                                                                        >  1:187494/1‑251 (MQ=255)
                                     cGGTAAGGCCCGGCGCGGCTGTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTTCGTTGCCAGTACGTGGTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGt                                                                                                                                                                                                        <  2:56990/251‑1 (MQ=255)
                                      ggTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGCCGCTCCCCTGTTACGGGTGATTCAGGTCATACAGtt                                                                                                                                                                                                       >  1:637817/1‑251 (MQ=255)
                                      ggTAATGCCCGGCGCGGCTTTGCCCTGGACGTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGtt                                                                                                                                                                                                       <  1:338402/251‑1 (MQ=255)
                                       gatattcgcGGCGCGTCTTGTCCCTGGAATTCAAGTTTCTTGTGTTGTGCGCTGTCCTGTTTCAGGGTTTGCGTGGCCAGTACGTGTTCAACTTCGGAAGTGACGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGGCGCCGCAGACAATGGGCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGAGCCCCTGATACGGGTGATTCAGGTCATACAGTTc                                                                                                                                                                                                      <  2:391052/243‑1 (MQ=255)
                                        gaaTGCCCGGCGCGTCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTTCGTGGTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGGTCCGTGCCCATGACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCg                                                                                                                                                                                                     <  1:37251/250‑1 (MQ=255)
                                              ccGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGcc                                                                                                                                                                                                >  1:479775/1‑250 (MQ=255)
                                               cGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCa                                                                                                                                                                                               >  1:94637/1‑250 (MQ=255)
                                                   gcgGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTc                                                                                                                                                                                          <  1:235597/251‑1 (MQ=255)
                                                          tGCCCTGGACTTCAATTTTCATGGTTTTTGCGGGGCCCGGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAAGTCTTTGCa                                                                                                                                                                                    <  2:347963/250‑1 (MQ=255)
                                                           gCCCTGGATTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGGACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGt                                                                                                                                                                                  <  2:374878/251‑1 (MQ=255)
                                                                      cAATTATCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGGCAAACGCGCCGTGGGTGGCGGTATGCGAGGCGTCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAg                                                                                                                                                                        <  1:689594/250‑1 (MQ=255)
                                                                           ttCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGGGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATGCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCt                                                                                                                                                                   <  2:330270/250‑1 (MQ=255)
                                                                            tCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCAGCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTg                                                                                                                                                                  <  2:174044/250‑1 (MQ=255)
                                                                                         gcgGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGCGTTCAACTTCGGAAGTGCCGCTACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCACGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACCGTTCGCCGCCAAATTCTTTGCAGTTTTTTATCATTTCCTGCATCTGGATAcgcg                                                                                                                                                    >  2:168996/1‑251 (MQ=255)
                                                                                          cgGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCACCTTCGGAAGTGCCGCTACCAAAAGCCCGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCACCGCAGACACTGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTAAGTGGACGGTCCCCTGATACGGGTGATTAGGGTCCTACAATTCGACGCCCAAATCTTTGCGGTTTTTGATCACTTCCTTCATccaggtcctcgc                                                                                                                                                   >  2:156265/1‑240 (MQ=255)
                                                                                                tGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGAc                                                                                                                                                                                                   >  1:304358/1‑197 (MQ=255)
                                                                                                tGTTTCAGGGTTTGCGTGGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGAc                                                                                                                                                                                                   <  2:304358/197‑1 (MQ=255)
                                                                                                         gTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCGTCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGc                                                                                                                                    >  1:461173/1‑251 (MQ=255)
                                                                                                           ttgggTTGCCAGTACGTGTTCAACTTAGGAAGTGCCGATACCAAAGGCCAGTGCGGCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCa                                                                                                                                   <  2:479775/246‑1 (MQ=255)
                                                                                                           ttGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCa                                                                                                                                   <  1:492477/250‑1 (MQ=255)
                                                                                                            tGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGAGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGAGATTCAGGTCATACCGTTCGACGCAAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCAtt                                                                                                                                 >  1:599459/1‑251 (MQ=255)
                                                                                                            tGCGTTGCCAGTACGTGTTCAACGTCGTAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCAt                                                                                                                                  <  2:576983/250‑1 (MQ=255)
                                                                                                               gTTGCCAGTACGTGGTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTTCATCTGGATACGCGCCATTTCACCGCAGGCATTAAt                                                                                                                              <  2:264902/251‑1 (MQ=255)
                                                                                                                 tGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGTACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATg                                                                                                                             <  1:256447/250‑1 (MQ=255)
                                                                                                                          gtgttcaactcaggcggTGCCGATACCAAACACCACTGCACCAAACGCGCCGTGGGTGGCGGTATGCGAGTAGCCGCAGACAATGGTCATCCCCGGCAAGGAGACGCCCTGTTCCAGGCCCATTACGTGCACGCTCCACTGTTACGGGGGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCCGCATCTGGATCCGCCCCATTTACCCGAAGGCATTAATGTCTTTGGtat                                                                                                                   >  2:653719/17‑249 (MQ=255)
                                                                                                                               cAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTa                                                                                                              <  1:285307/251‑1 (MQ=255)
                                                                                                                                aaCTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGt                                                                                                                                                                                                        <  1:445857/160‑1 (MQ=255)
                                                                                                                                aaCTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGATCATACAGt                                                                                                                                                                                                        >  2:445857/1‑160 (MQ=255)
                                                                                                                                     cGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGAGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACg                                                                                                        >  2:591142/1‑251 (MQ=255)
                                                                                                                                       gAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGtt                                                                                                      >  1:508312/1‑251 (MQ=255)
                                                                                                                                         aGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTtgt                                                                                                    <  1:288955/251‑1 (MQ=255)
                                                                                                                                          gTGCCTATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGGATGCGAGTCTCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTTCGTGTACGATCCCCTGATACTGGTGATTCAGGTCATACAGTTCGACTCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTtgtg                                                                                                   <  1:123036/251‑1 (MQ=255)
                                                                                                                                               gATACCAAAGGCCAGTGCGCCAATCGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCAACCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCa                                                                                              <  2:596130/251‑1 (MQ=255)
                                                                                                                                                aTACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGCGGACGATCCCCTGATACGGGTGATTCAGGTCTTACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGCTACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCAt                                                                                             >  1:655976/1‑251 (MQ=255)
                                                                                                                                                    cAAATGCCCGTGCGCCAAACTCGCCGTGGGTGGCGGTATGCTGGTCGCCGCAGACAATGGTCATTCCCGGCAAGGTTCCGCCCTGTGCCGGGCCCATTACGTGGACGATCCCCTTATACGGGTGATTCAGGTCATACATTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGAGACGCGCCATTTCACCGCAGGAATTAATGTCTTGGGTCTGGGTAGAGACGTTGTGATCCATgtta                                                                                         <  2:73811/251‑4 (MQ=255)
                                                                                                                                                    cAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTa                                                                                         <  1:19037/251‑1 (MQ=255)
                                                                                                                                                      aaGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTAACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGAGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGc                                                                                       >  1:437707/1‑251 (MQ=255)
                                                                                                                                                                       cgcgCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCt                                                                      <  1:231121/251‑1 (MQ=255)
                                                                                                                                                                            cGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCAGCCCAGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGAGGACGATCCCCTGATACGGGTGATTAAGGACATACAGTTCGACGACAAATTCTTTGAAGTTTTTGATAATTTACTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGAAGTTGTGATCCATGGAAGCGAAGTATTTGCCCGGCTGACGt                                                                 >  1:468344/1‑251 (MQ=255)
                                                                                                                                                                                 gtggCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCgg                                                            <  1:577818/251‑1 (MQ=255)
                                                                                                                                                                                         aTGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTGGCCCGGCTGCCGTACCGGGCGACCGt                                                    >  2:56299/1‑251 (MQ=255)
                                                                                                                                                                                             gAGTCGCCGCAGACAATTGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGATGACGTACCGGGCGACCGTGGgc                                                >  2:17863/1‑251 (MQ=255)
                                                                                                                                                                                              aGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCGTGGTAGCGAAGGTTTTGCCGGGCTGACGTACCGGGCGACCGTGGgcg                                               >  2:193565/1‑251 (MQ=255)
                                                                                                                                                                                                ccgccgcAGACAATGGTCACCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTGTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGgcgcg                                             <  2:598399/250‑1 (MQ=255)
                                                                                                                                                                                                         ccAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCTACCGTGGGCGCGCAGACCATc                                    <  2:315658/250‑1 (MQ=255)
                                                                                                                                                                                                                tCATCCCCGGCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGCTACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGcc                             >  1:617173/1‑251 (MQ=255)
                                                                                                                                                                                                                        ggCAAGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCGGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGTCCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGa                     <  1:503306/251‑1 (MQ=255)
                                                                                                                                                                                                                           aaGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGAGGt                  <  1:168996/251‑1 (MQ=255)
                                                                                                                                                                                                                            aGGTGACGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGCCGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGACCGCCTGCGGTGAggcc                 >  2:527125/1‑249 (MQ=255)
                                                                                                                                                                                                                                gACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGAGGTCACtt             <  1:193565/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aCGCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGAGGTCACTTc            <  1:438663/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cGCCCTGTTCCGGGACCATTACGTTGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTGTTAAGTTTTTGATCAGTTCCGGCATCTGGAGACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTGTGCCCGGCGGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGAGGTCACTTCa           <  2:599459/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gCCCTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGAGGTCACTTCAt          <  1:92835/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cTGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGATCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGGGGTCACTTCATGcc       >  1:480640/1‑250 (MQ=255)
                                                                                                                                                                                                                                       tGTTCCGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGAGGTCACTTCATGCAc      >  1:164465/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ggTTCCGGTCCCATTACGTGGACGATCCCCTGATACGTGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGAGCAGTTCCTGCATCTGGATACGCGCCATTTCACCTCAGGCATTAATGTCTGTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGGACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGAGGTCACTTCATGCAc      <  2:634037/250‑1 (MQ=255)
                                                                                                                                                                                                                                           ccGGGCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGAGGTCACTTCATGCACCAgg  <  2:94637/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
CCGGTTTTACCGATAATTGCCAGCACGATATCTTTTGCGGTAATGCCCGGCGCGGCTTTGCCCTGGACTTCAATTTTCATGGTTTTTGCGCGGCCCTGTTTCAGGGTTTGCGTTGCCAGTACGTGTTCAACTTCGGAAGTGCCGATACCAAAGGCCAGTGCGCCAAACGCGCCGTGGGTGGCGGTATGCGAGTCGCCGCAGACAATGGTCATCCCCGGCAAGGTGACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATGGTAGCGAAGGTTTTGCCCGGCTGACGTACCGGGCGACCGTGGGCGCGCAGACCATCGAACGCCTGCGGTGAGGTCACTTCATGCACCAGG  >  NC_000913/80290‑80775

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