Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,725,176 T→G E48A (GAG→GCG)  glyQ ← glycine‑‑tRNA ligase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,725,1760TG97.5% 108.6 / NA 40E48A (GAG→GCG) glyQglycine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (21/18);  minor base A (1/0);  total (22/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 = 1.00e+00

CCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAG‑CTCTTTCAGAGAACCG‑AGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCAT  >  NC_000913/3724955‑3725421
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
cacAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAG‑CTCTTTCAGAGAACCG‑AGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATATTAGGGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGt                                                                                                                                                                                                                           <  1:25278/248‑1 (MQ=255)
      tATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAG‑CTCTTTCAGAGAACCG‑AGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGt                                                                                                                                                                                                                           <  1:261183/244‑1 (MQ=255)
      tATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAG‑CTCTTTCAGAGAACCG‑AGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGGTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGGCCCCGCGCGCGCAGGCCGGTCATTGGGGgagggggcagtcgg                                                                                                                                                                                                                    >  2:261183/1‑240 (MQ=255)
          ttCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAG‑CTCTTTCAGAGAACCG‑AGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGGGAGAGGTTCCCGCGCCgc                                                                                                                                                                                                                >  1:187195/1‑250 (MQ=255)
            ccACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAG‑CTCTTTCACAGAAACG‑AGGTAAAGATACTGAATATTGTCCGGCGATGGATTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTTGGGTTTTCGCAGTAAAGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCTTCGGTTCTTGCCCCAGAGCGCGACGACAGGTCATTGGGTGAGAAGTTCCCGACCCGACt                                                                                                                                                                                                              >  2:408889/1‑251 (MQ=255)
                gAAACGGATGTCGTGAATAGTCTGGTCCATGCCCAG‑CTCTTTCAGAGAACCG‑AGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACAATCGGTCGGGCGACGAGAAGGATGAACATAAGCAGCCGCCATCGGTTCAGGCCCCAGCACACGCAGACAGGAACTTGGGTGAGAGGTTCCGGCGCCGACTTcc                                                                                                                                                                                                           >  2:374365/1‑250 (MQ=255)
                 aaaCGGATGTCGTGAATAGTCGGGTCCATGCCCAG‑CTCTTTCAGAGAACCG‑AGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCAt                                                                                                                                                                                                         >  1:277514/1‑251 (MQ=255)
                         gTCGTGAATAGTCGGGTCCATGCCCAG‑CTCTTTCAGAGAACCG‑AGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAAt                                                                                                                                                                                                  <  1:101830/250‑1 (MQ=255)
                              gAATAGTCGGGTCCATGCCCAG‑CTCTTTCAGAGAACCG‑AGGTACAGCTCCTGAATATTGTCCGGCGATGTCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTGCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGa                                                                                                                                                                                            <  1:210141/251‑1 (MQ=255)
                              gAATAGTCGGGGACATGCCCCG‑CTCTTTCCGTGAACCG‑AGGTACAGCTCCTGAATATTGTCCGGCGGTGGCTTAATGACCACCTGGAACTGCTAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGCCGAGAAGGCTGAACCTAAGCAACCGCCCTCGGTTCTGTCCCCAGCGCGGGCAGACCGGGCCTTGGGTGCGAGGTTCCCGCGCCGCCTTGCATGTCCAATGGttac                                                                                                                                                                                            >  2:406793/1‑249 (MQ=255)
                              gAATAGTCGGGAACATGCCCAG‑CTCTTTCAGAGAACCG‑AGGTACAGCTCCTGAATAATGTCCGGCGATGGCTTAAAGACCACCTGGAACTGATAGTAGTGATGTAACCGGTTGGGGATTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACCTAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGGCATGTCATTGGGGGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGttt                                                                                                                                                                                             >  2:283150/1‑249 (MQ=255)
                                  aGTCGGGTCCAGGCCCAGGCT‑TGTCAGAGAACCG‑AGGTACAGCTCCGGAATATTGTCCGGCGAGGGTGTAATGACCACCTGGAACTGATAGTAGTGCGGGAAACGGTTGGGGTTTTCGCCTTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCTGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGGGCCGACTTCCATGTCCAATGGTTGAACaa                                                                                                                                                                                        <  2:149400/251‑1 (MQ=255)
                                         tCCATGCCCAG‑CTCTTTCAGAGAACCG‑AGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGc                                                                                                                                                                                  >  1:201297/1‑250 (MQ=255)
                                               cccAG‑CTCTTTCAGAGAACCGCGGGTTC‑GCTCCTGAATATTGTCCGGCGATGTCTTAATGACCTCCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTGTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCTGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGTCGACTTCCATGTCCAATGGTTGAATAATGGTGCAGCCCTg                                                                                                                                                                           <  2:328286/251‑1 (MQ=255)
                                                  aG‑CTCTTTCAGAGAACCG‑AGGTACAGCTCCTGCATATTGTCCGGCGATGGCTTAATGACCACCTGGACCTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGAAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGCGAGAGGTTCCCGCGCCGCCTTCCATGTCCAATGGTTGACCAATGGTGAAGCCCTGGCg                                                                                                                                                                        >  2:235942/1‑251 (MQ=255)
                                                                          aCAGCTCCTTAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTc                                                                                                                                                  <  1:303653/251‑1 (MQ=255)
                                                                            aGCTCCGGAATATTGGCCGTCGATGGCGTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGATTTCGCCGTAGCGACCATCGGTCGGGAGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGACCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAg                                                                                                                                                <  1:45602/251‑1 (MQ=255)
                                                                                    aaTATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGGAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCAGCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGc                                                                                                                                        <  2:408867/251‑1 (MQ=255)
                                                                                     atatTGTCCGGCGATGGCTTCATAACCACCTGGAACTGCTAGTAGTGCTGTAAACGGATGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGCCGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGAGCACGCAGACAGGTAATTGGGTGAGAGGTTCCCGCGCCGACTTCCCTGTCCAATGGTTGAACAATGGTGCAGCACTGGCGAGCCCAGTCATCCTGTAAGGACAgagtcaagac                                                                                                                                       >  2:451718/1‑242 (MQ=255)
                                                                                          gTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTgg                                                                                                                                   >  1:455669/1‑250 (MQ=255)
                                                                                           tCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGGGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGACCTGGa                                                                                                                                  >  2:150763/1‑250 (MQ=255)
                                                                                             cGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGCCCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGTTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAgg                                                                                                                               >  1:65856/1‑251 (MQ=255)
                                                                                                cgaaggctaaaagataaccTGGAAATGATAGTAATGATGTAAACGGTTGGGGTTTTCGCCGCAGCGCCCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCCGCCGACATCGGTTCTGGCCACAGCGCGCGCAGACAGGTCATTGGGAGAGAGGTTCCCGCACAGACTTCCATGACCAATGGTTGAACCATGCAGCAGCCCTGGCGAGCCCAGACAACATGTATGGTCAGGATCAAGCCCTGGACGGTcc                                                                                                                            >  2:470234/17‑251 (MQ=255)
                                                                                                       ttAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATc                                                                                                                     >  1:241089/1‑251 (MQ=255)
                                                                                                       ttAATGACCACCTGGAACTGATAGTAGGGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGAGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTCATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATc                                                                                                                     >  2:416408/1‑251 (MQ=255)
                                                                                                        tAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATc                                                                                                                     >  1:332581/1‑250 (MQ=255)
                                                                                                               caccTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTGCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGTCCCCAGCGCGCGCAGACAGGTCATTGGGGGAGAGGTGCCCGCGCCGACTTCCATGTTCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACtttt                                                                                                             <  1:451718/251‑1 (MQ=255)
                                                                                                                   gtgatccgATAGTAGGGCTGTAAACGGTGGGGGTTTTCGCCGTAGCGACCATCGTGCGGGCGCCGATACGGGTGAACATGATCTGCTGCCATCGGTTCTGGCCCTAGCGCGCGCAGACAGGTCATTGGGTGAGAGGGTCCCGCGCCGTCTTCCATGTCCAGTGGTTGTTCAATGGTGCAGCGCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGAGCAAGCCCTGGAAGGTCCTGGTATCCAACTTTTGCat                                                                                                         <  2:65856/244‑1 (MQ=255)
                                                                                                                             agtagtGCTGTAAACGGTTGGGGTTTTCGCCGTGGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCg                                                                                                <  2:183672/250‑1 (MQ=255)
                                                                                                                                      ggaaaCGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATAc                                                                                      <  1:374365/249‑1 (MQ=255)
                                                                                                                                             gTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTgg                                                                                >  1:133491/1‑250 (MQ=255)
                                                                                                                                                         gCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGTTTGAACAATGGTACAGCCCTGGCGAGCCCAGGAATCCTGTAATGTCAGGATCAAGACCTGACGGGTCCTGGTATCACACTTTTGCATATTATTTCGATCTGGATACGTTTGGATTTCaaggaagg                                                                   >  2:124158/1‑251 (MQ=255)
                                                                                                                                                            gTAGCGACCATCGGTCGGGCGACGAGGAGGCTGACCATAAGCAGCCGCCATCGGTTCTAGCCCCAGCGCGCGCAGACAGGTCATTGGGTGCGAGGTTCCCGCGCCGCCTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTCAGGTCAGGATCACGCCCTGGAAGGGCCTGGTATCAAACTTTTGACTCTTATTACGTGCTGGTTAAGTGTGGATTTAaaggaaggaat                                                                >  2:280282/1‑248 (MQ=255)
                                                                                                                                                              agctaCCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTAGCAAACTGTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCa                                                              <  2:277514/247‑1 (MQ=255)
                                                                                                                                                              aGCGACCAGCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGGAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCa                                                              <  1:6341/251‑1 (MQ=255)
                                                                                                                                                                       cggtcggGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCg                                                     <  2:36851/251‑1 (MQ=255)
                                                                                                                                                                            gggCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCtgg                                                 >  1:6475/1‑250 (MQ=255)
                                                                                                                                                                                                                ggTTCTGGACCCAGCGCGCGCAGACTGGTCATTGGGTGAGAGGTTCCCGCGCCGACTGCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGACCTGGTATCAAACTTTTGAATATTATTTCGTGCTGGATACGTGAGGGTTTAAAGGAAGTGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAggg            >  2:67316/1‑251 (MQ=255)
                                                                                                                                                                                                                         cccAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCa   <  1:448694/251‑1 (MQ=255)
                                                                                                                                                                                                                          ccAGCGCGCGCAGACGGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATGTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCGTAACGCAt  <  2:241089/251‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
CCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAG‑CTCTTTCAGAGAACCG‑AGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCAT  >  NC_000913/3724955‑3725421

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