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

GCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCGTAAA  >  NC_000913/83341‑83820
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
gcgggcgATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCTCGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCGTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATg                                                                                                                                                                                                                                         <  2:138582/251‑1 (MQ=255)
                       cAAAATCGCTCGGCGCAGAGCCGGGCAAACCGTATTCTATCACGTCAACACCCTTACGGGCAAGGGCCGGCCCAATTTGCAGTTTTTCTTTCACCATCAAGCTTGCCTGAAAAGCCTGTTCACCGTAGCGCAATGAGGTATCGCCAATATTGTCTTGTTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCTCCTCGCTTCGGGCATAAAAAAAACCAGCGCAATGGCGCGGGTTTTTTGTTTGACTGc                                                                                                                                                                                                                  >  2:210571/1‑251 (MQ=255)
                       cAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTg                                                                                                                                                                                                                   <  1:6018/250‑1 (MQ=255)
                         aaaTCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGt                                                                                                                                                                                                                <  1:12535/251‑1 (MQ=255)
                          aaTCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATCAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTg                                                                                                                                                                                                               >  2:97945/1‑251 (MQ=255)
                                ccGGCGAAGAGACGGGGAAACCGACTTCCATCCCGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCACGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCt                                                                                                                                                                                                         >  2:405679/1‑251 (MQ=255)
                                           ccGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCGATTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCACAATGGCTCGGGTTTTTGGGTTGACTGCGTGCTGGCTTAATGCTGGAt                                                                                                                                                                                              <  2:87157/250‑1 (MQ=255)
                                           aCGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGTTATCGAAAATAATGACTTGCTGCCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCTCAATGGCGCGTGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGAt                                                                                                                                                                                              <  1:407346/251‑1 (MQ=255)
                                           aCGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGAt                                                                                                                                                                                              >  2:157561/1‑251 (MQ=255)
                                            cGGGTAAACCGACTTCCATCACGTCAACACCCATACGCTCAATGGCCAGCGCAATTTGCAGTATTTCTTTCACACTCAAGCGTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCGTTTAGAGCGCCGCGCTTCGGGCATAAAAAAAACCCGCGCAAGGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATg                                                                                                                                                                                             <  2:226513/251‑1 (MQ=255)
                                               ggAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCg                                                                                                                                                                                          >  1:339083/1‑251 (MQ=255)
                                                gAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGc                                                                                                                                                                                         >  1:127170/1‑251 (MQ=255)
                                                  aaCCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTc                                                                                                                                                                                       <  1:131674/251‑1 (MQ=255)
                                                       aCTTCCATCACGTCAACACCCATACGCTCAAGGTCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCg                                                                                                                                                                                  <  1:128944/251‑1 (MQ=255)
                                                                  gTCAACACCCATACGCTCAAGTGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCa                                                                                                                                                                       <  1:424732/251‑1 (MQ=255)
                                                                        acCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaagaag                                                                                                                                                                 <  1:330966/251‑1 (MQ=255)
                                                                         cccATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACGCAATCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGATCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCGTAATGCTGGAGGCCGCTCATTCGGCTTCCGCGCAAagaaga                                                                                                                                                                <  2:382205/251‑1 (MQ=255)
                                                                                   cAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagt                                                                                                                                                      <  2:127170/251‑1 (MQ=255)
                                                                                         ccAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtagt                                                                                                                                                <  1:157561/251‑1 (MQ=255)
                                                                                           aGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTCATTGAAACCCCTACTAACCGTCGAGAAATCTGGCACCGAAAATACTGAATTACTGGCTCACAGATTGGCTCCATGTCCCTTTTAGAGCGCCTCGCTTCGGGCATACAACAAACCCCCGCAATGGCGCGGGTTTTTTGTTTGCCTGCGTGCTCGCTTAATGCTGGATGCCGCCCACTCGTCTACCGCGCCAAGAAGATGCGTTtagtagtcgtcg                                                                                                                                              >  1:396984/1‑246 (MQ=255)
                                                                                                      cAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCCAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGc                                                                                                                                   >  1:271607/1‑251 (MQ=255)
                                                                                                                ttcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTGAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtg                                                                                                                           <  1:405679/249‑1 (MQ=255)
                                                                                                                tgcacaCTCAAGCTTGCGTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTGGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTGTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtgt                                                                                                                          <  1:20856/248‑1 (MQ=255)
                                                                                                                   acacTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGt                                                                                                                      >  1:80995/1‑251 (MQ=255)
                                                                                                                       tCAAGCTTGCCTGTAACGCATGTTAACCGTCGCGCCATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGCTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGATTCGGGCATAAAAAAAACCCGCGCCATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAGGTGAGTCAtt                                                                                                                  >  2:281032/1‑251 (MQ=255)
                                                                                                                                                aCCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGat                                                                                         >  2:231152/1‑250 (MQ=255)
                                                                                                                                                  cGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATg                                                                                      >  1:320135/1‑251 (MQ=255)
                                                                                                                                                          aaTGTGGTATCGAAAATACTGACTTGCCGGCTCATGGTTTGGGACCTTGTATCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCACTGGCGCGGGATTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCttttcg                                                                               >  1:368224/1‑248 (MQ=255)
                                                                                                                                                          aaTGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGATGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGGGAGTCATTAAATCAGCCTCCGATGAATGCGATATGCtttt                                                                                 >  2:431316/1‑248 (MQ=255)
                                                                                                                                                                 tATCGAAAATACTGACTTGCTGGCTCATCGTTTGGGTCCTTGTCTCTTTTAGTCCGCCTCGCTACGGGCATCACAAAAACGCGCGCCATGGCGCGGGTTTTTTGTTTGATTGCATGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAATGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGCTGTGAGTCATTAAATCAGCACCAGATGCCTGCGATCTGCTTTTAGAGGTACTg                                                                       >  2:215761/1‑251 (MQ=255)
                                                                                                                                                                            atGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAAc                                                            <  1:100348/251‑1 (MQ=255)
                                                                                                                                                                             tGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTCGCGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCg                                                                                                                                                                                                                 >  1:410531/1‑102 (MQ=255)
                                                                                                                                                                             tGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTATGTTTGACTGCg                                                                                                                                                                                                                 <  2:410531/102‑1 (MQ=255)
                                                                                                                                                                              gACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACgg                                                          >  2:305182/1‑251 (MQ=255)
                                                                                                                                                                                     tGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCCCTCGTCTAACGCGCAAAGAAGATGCGTTTAGTAGTAGTAGCCCGATAAAGCGAACGATGTGAGTCATTAAATAAGCTCCAGATGAATg                                                                                             >  1:207199/1‑209 (MQ=255)
                                                                                                                                                                                     gggCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATg                                                                                             <  2:207199/208‑1 (MQ=255)
                                                                                                                                                                                                                            cGCTTCGGGCATAAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGGCCCGCCAGCGTGAGTTc              <  1:305182/249‑1 (MQ=255)
                                                                                                                                                                                                                                    ccatcaaaaaagCCCGCGCAATGGCGCGTGTTTTTTTTTTGGCTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGTTGCGTTTAGGAGTAGTAGACCGATAAAGCGAACGATGTGGGTCATTAAATCAGCTCCAGATGAATGCGATTTGCTTTTAGAGTTACTGGATACAAAAACAGATGTCACCCCTGACGCAATAAAAACGTCCCGCCAGCGTGATTTCTGCATCCGTa    <  1:74609/246‑1 (MQ=255)
                                                                                                                                                                                                                                      ataaaaaaaaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTGGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATTTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCGTaaa  <  1:431316/248‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
GCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCGTAAA  >  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: