Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 83,579 (A)7→8 intergenic (‑50/+43) leuA ← / ← leuL 2‑isopropylmalate synthase/leu operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091383,5791.A91.9% 112.2 / 1.7 37intergenic (‑50/+43)leuA/leuL2‑isopropylmalate synthase/leu operon leader peptide
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base A (15/19);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

GCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAA‑CCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGC‑TCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCGTAAA  >  NC_000913/83341‑83820
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
gcgggcgATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCTCGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCGTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATg                                                                                                                                                                                                                                         <  2:138582/251‑1 (MQ=255)
                       cAAAATCGCTCGGCGCAGAGCCGGGCAAACCGTATTCTATCACGTCAACACCCTTACGGGCAAGGGCCGGCCCAATTTGCAGTTTTTCTTTCACCATCAAGCTTGCCTGAAAAGCCTGTTCACCGTAGCGCAATGAGGTATCGCCAATATTGTCTTGTTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCTCCTCGCTTCGGGCAT‑A‑AAAAAAACCAGCGCAATGGCGCGGGTTTTTTGTTTGACTGc                                                                                                                                                                                                                  >  2:210571/1‑251 (MQ=255)
                       cAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTg                                                                                                                                                                                                                   <  1:6018/250‑1 (MQ=255)
                         aaaTCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGt                                                                                                                                                                                                                <  1:12535/251‑1 (MQ=255)
                          aaTCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATCA‑AAAAAA‑CCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTg                                                                                                                                                                                                               >  2:97945/1‑251 (MQ=255)
                                ccGGCGAAGAGACGGGGAAACCGACTTCCATCCCGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCACGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCt                                                                                                                                                                                                         >  2:405679/1‑251 (MQ=255)
                                           ccGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCGATTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCACAATGGCTCGGGTTTTTGGGTTGACTGCGTGCTGGCTTAATGCTGGAt                                                                                                                                                                                              <  2:87157/250‑1 (MQ=255)
                                           aCGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGTTATCGAAAATAATGACTTGCTGCCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCTCAATGGCGCGTGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGAt                                                                                                                                                                                              <  1:407346/251‑1 (MQ=255)
                                           aCGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGAt                                                                                                                                                                                              >  2:157561/1‑251 (MQ=255)
                                            cGGGTAAACCGACTTCCATCACGTCAACACCCATACGCTCAATGGCCAGCGCAATTTGCAGTATTTCTTTCACACTCAAGCGTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCGTTTAGAGCGCCGCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAAGGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATg                                                                                                                                                                                             <  2:226513/251‑1 (MQ=255)
                                               ggAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCg                                                                                                                                                                                          >  1:339083/1‑251 (MQ=255)
                                                gAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGc                                                                                                                                                                                         >  1:127170/1‑251 (MQ=255)
                                                  aaCCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTc                                                                                                                                                                                       <  1:131674/251‑1 (MQ=255)
                                                       aCTTCCATCACGTCAACACCCATACGCTCAAGGTCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCg                                                                                                                                                                                  <  1:128944/251‑1 (MQ=255)
                                                                  gTCAACACCCATACGCTCAAGTGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCa                                                                                                                                                                       <  1:424732/251‑1 (MQ=255)
                                                                        acCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaagaag                                                                                                                                                                 <  1:330966/251‑1 (MQ=255)
                                                                         cccATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACGCAATCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGATCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCGTAATGCTGGAGGCCGCTCATTCGGCTTCCGCGCAAagaaga                                                                                                                                                                <  2:382205/251‑1 (MQ=255)
                                                                                   cAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagt                                                                                                                                                      <  2:127170/251‑1 (MQ=255)
                                                                                         ccAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtagt                                                                                                                                                <  1:157561/251‑1 (MQ=255)
                                                                                           aGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTCATTGAAACCCCTACTAACCGTCGAGAAATCTGGCACCGAAAATACTGAATTACTGGCTCACAGATTGGCTCCATGTCCCTTTTAGAGCGCCTCGCTTCGGGCAT‑ACAACAAA‑CCCCCGCAATGGCGCGGGTTTTTTGTTTGCCTGCGTGCTCGCTTAATGCTGGATGCCGCCCACTCGTCTACCGCGCCAAGAAGATGCGTTtagtagtcgtcg                                                                                                                                              >  1:396984/1‑246 (MQ=255)
                                                                                                      cAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCCAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGc                                                                                                                                   >  1:271607/1‑251 (MQ=255)
                                                                                                                ttcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTGAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtg                                                                                                                           <  1:405679/249‑1 (MQ=255)
                                                                                                                tgcacaCTCAAGCTTGCGTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTGGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTGTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtgt                                                                                                                          <  1:20856/248‑1 (MQ=255)
                                                                                                                   acacTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGt                                                                                                                      >  1:80995/1‑251 (MQ=255)
                                                                                                                       tCAAGCTTGCCTGTAACGCATGTTAACCGTCGCGCCATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGCTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGATTCGGGCAT‑A‑AAAAAAACCCGCGCCATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAGGTGAGTCAtt                                                                                                                  >  2:281032/1‑251 (MQ=255)
                                                                                                                                                aCCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGC‑TCCAGATGAATGCGat                                                                                         >  2:231152/1‑250 (MQ=255)
                                                                                                                                                  cGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGC‑TCCAGATGAATGCGATATg                                                                                      >  1:320135/1‑251 (MQ=255)
                                                                                                                                                          aaTGTGGTATCGAAAATACTGACTTGCCGGCTCATGGTTTGGGACCTTGTATCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCACTGGCGCGGGATTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGC‑TCCAGATGAATGCGATATGCttttcg                                                                               >  1:368224/1‑248 (MQ=255)
                                                                                                                                                          aaTGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGATGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGGGAGTCATTAAATCAGCCTCC‑GATGAATGCGATATGCtttt                                                                                 >  2:431316/1‑248 (MQ=255)
                                                                                                                                                                 tATCGAAAATACTGACTTGCTGGCTCATCGTTTGGGTCCTTGTCTCTTTTAGTCCGCCTCGCTACGGGCATCA‑CAAAAA‑CGCGCGCCATGGCGCGGGTTTTTTGTTTGATTGCATGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAATGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGCTGTGAGTCATTAAATCAGC‑ACCAGATGCCTGCGATCTGCTTTTAGAGGTACTg                                                                       >  2:215761/1‑251 (MQ=255)
                                                                                                                                                                            atGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGC‑TCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAAc                                                            <  1:100348/251‑1 (MQ=255)
                                                                                                                                                                             tGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTCGCGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCg                                                                                                                                                                                                                 >  1:410531/1‑102 (MQ=255)
                                                                                                                                                                             tGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTATGTTTGACTGCg                                                                                                                                                                                                                 <  2:410531/102‑1 (MQ=255)
                                                                                                                                                                              gACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGC‑TCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACgg                                                          >  2:305182/1‑251 (MQ=255)
                                                                                                                                                                                     tGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCCCTCGTCTAACGCGCAAAGAAGATGCGTTTAGTAGTAGTAGCCCGATAAAGCGAACGATGTGAGTCATTAAATAAGC‑TCCAGATGAATg                                                                                             >  1:207199/1‑209 (MQ=255)
                                                                                                                                                                                     gggCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGC‑TCCAGATGAATg                                                                                             <  2:207199/208‑1 (MQ=255)
                                                                                                                                                                                                                            cGCTTCGGGCAT‑A‑AAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGC‑TCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGGCCCGCCAGCGTGAGTTc              <  1:305182/249‑1 (MQ=255)
                                                                                                                                                                                                                                    ccatca‑aaaaag‑CCCGCGCAATGGCGCGTGTTTTTTTTTTGGCTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGTTGCGTTTAGGAGTAGTAGACCGATAAAGCGAACGATGTGGGTCATTAAATCAGC‑TCCAGATGAATGCGATTTGCTTTTAGAGTTACTGGATACAAAAACAGATGTCACCCCTGACGCAATAAAAACGTCCCGCCAGCGTGATTTCTGCATCCGTa    <  1:74609/246‑1 (MQ=255)
                                                                                                                                                                                                                                      ataa‑aaaaaa‑CCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTGGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGC‑TCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATTTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCGTaaa  <  1:431316/248‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
GCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑A‑AAAAAA‑CCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGC‑TCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCGTAAA  >  NC_000913/83341‑83820

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