Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 79,040 A→G A138A (GCT→GCC)  leuD ← 3‑isopropylmalate dehydratase subunit LeuD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091379,0400AG97.5% 123.7 / ‑6.2 40A138A (GCT→GCC) leuD3‑isopropylmalate dehydratase subunit LeuD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (19/20);  minor base C (1/0);  total (20/20)
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.39e-01

AAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTTTCTCGTGCCAGCAAAATGGAAGCGCCCTGATACTGCGGGAAGTTCAGCACGAAGTCC  >  NC_000913/78819‑79277
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
aaGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCAc                                                                                                                                                                                                                  >  1:219175/1‑251 (MQ=255)
                          ggATTAATTCATAAACGCAGGTTGTTTTGATTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGACGGCGGAAGGCATCGCTGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AAAGGATACCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGAATCGCTTAATTTctgatcgg                                                                                                                                                                                        >  1:557746/1‑244 (MQ=255)
                               aaTTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGcagcag                                                                                                                                                                                    <  1:239085/250‑1 (MQ=255)
                                aTTCATAAACGCAGGTTGTTTTGCTTCATAAGAGGCAATGGCGTCGTCGCGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCCTCATGCAGTGGCGGCGGAAGGCATCGATGGAAAAGCGATAAGTTTTTGCTCCCGCTTTCACCTCTTGCACTTAAAGATCCACGGCGACAATA‑AGCCCCGGATTGGCTTTCACCCGCGCAAAATGTTCGTCCACTTATGccgagcttcatgttaacgacagcagc                                                                                                                                                                                   >  2:48743/1‑225 (MQ=255)
                                         cGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGGTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTaa                                                                                                                                                                         <  2:138190/251‑1 (MQ=255)
                                                 gatttGCTTCCTAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGCCGA‑AATGGATCCCCGGATTGGCTTTCCACAGCGCCACCAGTTCGTCCACTTCTGCATCGCTACATTTCACCGGCAGCAGCTTGTTGTTAAAGCTATTg                                                                                                                                                                 >  2:111749/3‑251 (MQ=255)
                                                   tttGCTTCATAAGCGGCAATGGCGGCGTCGTGCTGCAAGGTAAGCCCAATCCTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGc                                                                                                                                                                <  2:329523/250‑1 (MQ=255)
                                                          cATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTagaag                                                                                                                                                        <  1:304926/251‑1 (MQ=255)
                                                          cATAAGCGGCAATGGAGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGAATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCCCCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTAACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTaga                                                                                                                                                          >  2:81201/1‑249 (MQ=255)
                                                           ctAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCCCCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCCGCAGCTGGGTTTTTGAGCTATTGCCGTaga                                                                                                                                                          >  2:163370/2‑248 (MQ=255)
                                                                   cAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGc                                                                                                                                                >  1:366348/1‑250 (MQ=255)
                                                                    aaTGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCaa                                                                                                                                              <  1:145026/251‑1 (MQ=255)
                                                                       ggcttcgtcgtTCTGCAAGGTATGCACAATAAGGTGCAGACCGTTCATCTTGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAAc                                                                                                                                           <  2:455920/247‑1 (MQ=255)
                                                                         cgtcgtcgtGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTc                                                                                                                                         <  1:620327/251‑1 (MQ=255)
                                                                          gtcgtcgtGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCTGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGGCAGCAAAACTCg                                                                                                                                        <  1:473427/251‑1 (MQ=255)
                                                                             gtcgtGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGgcg                                                                                                                                     >  1:61375/1‑251 (MQ=255)
                                                                                     gcAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATca                                                                                                                               <  1:605645/249‑1 (MQ=255)
                                                                                                     actttCCAGACCGTTCAGCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGt                                                                                                             <  2:708783/247‑1 (MQ=255)
                                                                                                          ccAGACCGTTCAGCAGGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTc                                                                                                        <  2:405768/251‑1 (MQ=255)
                                                                                                          ccAGACCGTTCAGCAAGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTTTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAACCCGTAGTCGGTc                                                                                                        <  2:348939/251‑1 (MQ=255)
                                                                                                            aGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCa                                                                                                       >  1:608564/1‑250 (MQ=255)
                                                                                                             gACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAAt                                                                                                     <  1:641211/251‑1 (MQ=255)
                                                                                                                    catcatGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAgg                                                                                              >  1:85371/1‑251 (MQ=255)
                                                                                                                     atcatGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGAGCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAggg                                                                                             <  2:595785/251‑1 (MQ=255)
                                                                                                                                  ggcggAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCAcgc                                                                                >  1:51795/1‑251 (MQ=255)
                                                                                                                                     ggAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAACCAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCACAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGa                                                                              >  1:663861/1‑250 (MQ=255)
                                                                                                                                        aGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCAc                                                                                                                                                                                             >  2:486239/1‑136 (MQ=255)
                                                                                                                                        aGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCAc                                                                                                                                                                                             <  1:486239/136‑1 (MQ=255)
                                                                                                                                           cATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGAAGAGcc                                                                        >  1:98086/1‑250 (MQ=255)
                                                                                                                                                                              ttCACCTCGTTCGCCTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTTTCTCGTGCCAGCAAAATGGaa                                      <  2:608564/249‑1 (MQ=255)
                                                                                                                                                                                cACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTTTCTCGTGCCAGCAAAATGGAAgcgc                                  >  1:29633/1‑251 (MQ=255)
                                                                                                                                                                                        gcgcTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTTTCTCGTGCCAGCAAAATGGAAGCGCCCTGATAc                          <  2:219175/251‑1 (MQ=255)
                                                                                                                                                                                            ttCCAGATCCACGTCGA‑CATGGATCCCCGGATTCGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTACATCGCTTACTTTCACCGGCAGAAGCTGGTTGTTCAAGCTATTGCCGTAGAAGCTGTCAGCAAAACTCGGCGCAATCAACACTTTAAAACCGTCGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTTTCTCGTGCCAGCAAAATGGCAGCGCCCTGAtcccacg                      >  2:89230/1‑245 (MQ=255)
                                                                                                                                                                                                       cggcGA‑AATGTATCCCCGTATTGTCTTTCACCAGCTCAAAAAGTTCGGCCCCTTCGTCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCAGTAGAAGATGTCAGCACAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGTTCAATGCCCAGGCCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTGTCTCGTGCCAGCAAAATGGAAGCGCCCTGATACTGCGGGAAGTTCAGc           <  2:440470/248‑1 (MQ=255)
                                                                                                                                                                                                           ta‑aaTGGATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGGCAGCATAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGccc                                                                                                 <  1:39869/160‑1 (MQ=255)
                                                                                                                                                                                                           gA‑AATGGAACCACGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATAGCTTAATTTCACCGGCAGCAGCTGGTTGGTAACGCTATTGCCGTAGCAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACAGTAGTCGGTCAATGccc                                                                                                 >  2:39869/1‑161 (MQ=255)
                                                                                                                                                                                                                 ggctCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAAATCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCAAGGGCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTTTCACGTGCCAGCAAAATGGAAGCGCCCTGATACTGCGGGAAGTTCAGCACGAAGTcc  >  2:354279/4‑251 (MQ=255)
                                                                                                                                                                                                                 ggATCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTTTCTCGTGCCAGCAAAATGGAAGCGCCCTGATACTGCGGGAAGTTCAGCACGAAGTcc  >  1:502549/1‑251 (MQ=255)
                                                                                                                                                                                                                   aTCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTa                                                                                                                                                            >  2:405575/1‑95 (MQ=255)
                                                                                                                                                                                                                   aTCCCCGGATTGGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTa                                                                                                                                                            <  1:405575/95‑1 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
AAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGA‑AATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAGCAGCTGGTTGTTAAAGCTATTGCCGTAGAAGATGTCAGCAAAACTCGGCGCAATCACCACTTTAAAACCGTAGTCGGTCAATGCCCAGGGCGCGTGCTCACGCGAAGAGCCACAGCCGAAGTTTTCTCGTGCCAGCAAAATGGAAGCGCCCTGATACTGCGGGAAGTTCAGCACGAAGTCC  >  NC_000913/78819‑79277

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