Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,226,228 T→C intergenic (+57/‑6) ebgC → / → ygjI DUF386 domain‑containing evolved beta‑D‑galactosidase subunit beta/putative transporter YgjI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,226,2280TC100.0% 165.8 / NA 46intergenic (+57/‑6)ebgC/ygjIDUF386 domain‑containing evolved beta‑D‑galactosidase subunit beta/putative transporter YgjI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/20);  total (26/20)

GTAGTGGAATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGATGTGTTATGTCTGATACCAAACGTAATACAATCGGCAAATTCGGCTTGCTCTCGCTGACTTTTGCCGCCGTTTACAGCTTTAACAACGTTATCAATAATAATATTGAGCTTGGACTGGCCTCGGCACCGATGTTTTTCCTCGCGACGATTTTTTATTTTATTCCCTTCTGTCTGATCATCGCAGAATTTGTTTCGTTAAATAAAAACTCAGAAGCCGGTGTCTACGCGTGGGTAAAAAGTTCGCT  >  NC_000913/3225983‑3226472
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
gTAGTGGAATATTAGCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGtt                                                                                                                                                                                                                                                 <  1:258063/251‑1 (MQ=255)
          attatCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatcct                                                                                                                                                                                                                                       >  1:350634‑M1/1‑241 (MQ=255)
          attatCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatcc                                                                                                                                                                                                                                        <  1:313845‑M1/250‑10 (MQ=255)
                      aaaCTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatgg                                                                                                                                                                                                                           >  2:201409‑M1/1‑229 (MQ=255)
                           gACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaa                                                                                                                                                                                                                        >  1:16156‑M1/1‑224 (MQ=255)
                                         aaaaGGTTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATTTCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttc                                                                                                                                                                                                         <  1:74725‑M1/250‑41 (MQ=255)
                                           aaGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaaaggcgcttcac                                                                                                                                                                                                       >  1:465828‑M1/1‑208 (MQ=255)
                                           aaGGCTGCGGAGAAAACGATGAGGTCCACGCCGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAACGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGCGACGCGTTgattaatccttagccgggatggacaagcaaacggcgcttcac                                                                                                                                                                                                       >  2:189198‑M1/1‑208 (MQ=255)
                                                   ggcgaaacCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCCACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaa                                                                                                                                                                                               >  1:200942‑M1/4‑200 (MQ=255)
                                                   ggAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaa                                                                                                                                                                                               >  1:200950‑M1/1‑200 (MQ=255)
                                                      gaAACCGTTGATGTCCTCGAAGGGCAAATCGGTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAATTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTGTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcgcattatgaaacgc                                                                                                                                                                                           <  2:313707‑M1/251‑55 (MQ=255)
                                                                      aCGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaaca                                                                                                                                                                            >  1:34379‑M1/1‑181 (MQ=255)
                                                                              cAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacatctcatttc                                                                                                                                                                   <  1:193635‑M1/251‑79 (MQ=255)
                                                                              cAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttc                                                                                                                                                                   <  1:490281‑M1/251‑79 (MQ=255)
                                                                                                  ccATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCCAAAGGAGTTTGCCGCCACCGCTTCCCTACTCATTTTCGGAGCCGTGTTgataaatccttagcaggtatggacaagcaaacggcgcttcacattatgaaacgccgttttttattaacacctcatttcgactttatagcgttacgccg                                                                                                                                               <  1:479633‑M1/251‑99 (MQ=255)
                                                                                                  ccATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatgggaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccg                                                                                                                                               >  2:382006‑M1/1‑153 (MQ=255)
                                                                                                    aTGAATCGTATGGGTTTATCTGCAATAGCGCGGTAAAAAAAGGGGTTCTCAAAGTCACCAACGAAGATGGTTATTTCCATAAAAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggctcttcacattatgaaacgccgttttttataaacaactcatttcgactttatagcgttacgccgct                                                                                                                                             <  2:282651‑M1/251‑101 (MQ=255)
                                                                                                     tGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgctt                                                                                                                                            >  2:485406‑M1/1‑150 (MQ=255)
                                                                                                         gCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttga                                                                                                                                        <  2:311472‑M1/251‑106 (MQ=255)
                                                                                                             aTCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttga                                                                                                                                        >  1:304510‑M1/1‑142 (MQ=255)
                                                                                                                    tATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcagggatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatctccagt                                                                                                                             >  2:404986‑M1/1‑135 (MQ=255)
                                                                                                                     aTCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagccggtatggaaaagcaaacggcgcttcccattatgaaaagccgttttttattaccaacccatttcgactttatagcggtacgcctcttttgacgatctcccgtt                                                                                                                            >  2:32440‑M1/1‑134 (MQ=255)
                                                                                                                         gCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaaatcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcg                                                                                                                        >  1:392962‑M1/1‑130 (MQ=255)
                                                                                                                                cgcgGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcg                                                                                                                 >  1:234065‑M1/1‑123 (MQ=255)
                                                                                                                                        aaaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctc                                                                                                         <  1:312724‑M1/251‑137 (MQ=255)
                                                                                                                                             gTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttac                                                                                                    >  2:108783‑M1/1‑110 (MQ=255)
                                                                                                                                                     cAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggttt                                                                                              <  1:108783‑M1/249‑148 (MQ=255)
                                                                                                                                                               aTCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttgg                                                                                  <  1:480804‑M1/251‑160 (MQ=255)
                                                                                                                                                                 cGAAGTTGGTTATTTCCATAACAAATAACAANTTCGGCGGCAAAAGGAGTTTGCCGCCACCGCTTCCCTACTCAGTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaactgcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacagggtttcttaccgttggca                                                                                <  2:225910‑M1/251‑162 (MQ=255)
                                                                                                                                                                 cGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggc                                                                                 >  1:137319‑M1/1‑90 (MQ=255)
                                                                                                                                                                  gAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCCCCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcacacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagctcgccagttcgcggaagtcgttaagctcgttacagggttgcttaccgttggca                                                                                >  2:399997‑M1/1‑89 (MQ=255)
                                                                                                                                                                  gAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagccaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggca                                                                                >  1:147981‑M1/1‑89 (MQ=255)
                                                                                                                                                                  gAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggca                                                                                >  1:390590‑M1/1‑89 (MQ=255)
                                                                                                                                                                  gAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacatt                                                                                                                                                                                                    <  1:126936‑M1/134‑46 (MQ=255)
                                                                                                                                                                  gAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACACTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacctt                                                                                                                                                                                                    >  2:126936‑M1/1‑89 (MQ=255)
                                                                                                                                                                            aTTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTCTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgat                                                                     >  1:321293‑M1/1‑79 (MQ=255)
                                                                                                                                                                                    aacaaataacaaCTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgacgttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcg                                                              <  2:455368‑M1/250‑180 (MQ=255)
                                                                                                                                                                                       aaataacaaCTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcgata                                                           >  1:247174‑M1/1‑68 (MQ=255)
                                                                                                                                                                                              aacTACGGCGGCAAAAGGAGTTTGCCGCCAACGCTACCCTACTAATTTTCGGGGACGCGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgacacggcgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggacgtcgttacgctcgttaccggtttgcttaccgttggcgaactcaaagatgggggcgataaattcttc                                                   >  2:298577‑M1/1‑61 (MQ=255)
                                                                                                                                                                                                         cAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccg                                                                                      <  2:231880‑M1/205‑156 (MQ=255)
                                                                                                                                                                                                         cAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccg                                                                                      >  1:231880‑M1/1‑50 (MQ=255)
                                                                                                                                                                                                                   ttGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcgataaattcttccagcaactgcggcatcgctt                               <  1:125405‑M1/250‑211 (MQ=255)
                                                                                                                                                                                                                        gccACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcgataaattcttccagcaactgcggcatcgctttaccgt                         <  1:201409‑M1/251‑217 (MQ=255)
                                                                                                                                                                                                                                         cATTTTCGGAGACGTGTTgataaatccttancaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcgataaattcttccagcaactgcggcatcgctttaccgtggatcaactgacccgc         <  2:247174‑M1/250‑233 (MQ=255)
                                                                                                                                                                                                                                             ttCGGCGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcgataacttcttccagcaactgcggcatcgctttaccgtggatcaactgacccgcgtca     >  2:130219‑M1/1‑14 (MQ=255)
                                                                                                                                                                                                                                               cGGAGACGTGTTgataaatccttagcaggtatggaaaagcaatcggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcgataaattcttccagcaactgcggcatcgctttaccgtggatcaactgacccgcgtcaaag  <  2:352312‑M1/251‑240 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
GTAGTGGAATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGATGTGTTATGTCTGATACCAAACGTAATACAATCGGCAAATTCGGCTTGCTCTCGCTGACTTTTGCCGCCGTTTACAGCTTTAACAACGTTATCAATAATAATATTGAGCTTGGACTGGCCTCGGCACCGATGTTTTTCCTCGCGACGATTTTTTATTTTATTCCCTTCTGTCTGATCATCGCAGAATTTGTTTCGTTAAATAAAAACTCAGAAGCCGGTGTCTACGCGTGGGTAAAAAGTTCGCT  >  NC_000913/3225983‑3226472

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