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,2280TC98.2% 194.1 / NA 55intergenic (+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);  major base C (24/30);  minor base A (1/0);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGATGTGTTATGTCTGATACCAAACGTAATACAATCGGCAAATTCGGCTTGCTCTCGCTGACTTTTGCCGCCGTTTACAGCTTTAACAACGTTATCAATAATAATATTGAGCTTGGACTGGCCTCGGCACCGATGTTTTTCCTCGCGACGATTTTTTATTTTATTCCCTTCTGTCTGATCATCGCAGAATTTGTTTCGTTAAATAAAAACTCAGAAGCCGGTGTCTACGCGTGGGTAAAAAGTTCGCT  >  NC_000913/3225998‑3226472
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
cGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGTAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagc                                                                                                                                                                                                                                   <  1:83725‑M1/250‑15 (MQ=255)
  tgatgaAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTACACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAACAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATNTTCGGAGCCGGGTTaataaatccttagcagg                                                                                                                                                                                                                                >  2:297698‑M1/1‑234 (MQ=255)
   gatgaAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggt                                                                                                                                                                                                                               <  2:575397‑M1/251‑19 (MQ=255)
   gatgaAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCCACTCATTTTCGGAGACGGGTTgataaatccttagcaggt                                                                                                                                                                                                                               >  2:99677‑M1/1‑233 (MQ=255)
            gACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATAGCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaa                                                                                                                                                                                                                       <  1:276698‑M1/250‑27 (MQ=255)
                 tGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatc                                                                                                                                                                                                                                         <  2:332089‑M1/227‑9 (MQ=255)
                 tGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatc                                                                                                                                                                                                                                         >  1:332089‑M1/1‑219 (MQ=255)
                         tAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgctt                                                                                                                                                                                                          <  1:287018‑M1/250‑40 (MQ=255)
                          aaaaGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttca                                                                                                                                                                                                        <  1:304594‑M1/251‑42 (MQ=255)
                              ggCTGCGGAGAAACCGTTGAGGGCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacatt                                                                                                                                                                                                    <  2:573173‑M1/251‑46 (MQ=255)
                                cTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCAATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattat                                                                                                                                                                                                  >  1:94431‑M1/1‑204 (MQ=255)
                                  gCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatg                                                                                                                                                                                                 <  1:387906‑M1/250‑49 (MQ=255)
                                   cGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTTCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcccattatga                                                                                                                                                                                                >  2:90798‑M1/1‑201 (MQ=255)
                                     gagaAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaac                                                                                                                                                                                             <  1:297698‑M1/251‑53 (MQ=255)
                                                 aGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAACAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatg                                                                                                                                                                                                                            >  2:197685‑M1/1‑187 (MQ=255)
                                                 aGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatg                                                                                                                                                                                                                            <  1:197685‑M1/208‑22 (MQ=255)
                                                     ccACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaa                                                                                                                                                                              >  1:171295‑M1/1‑183 (MQ=255)
                                                                  aTCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgac                                                                                                                                                                <  2:109905‑M1/251‑82 (MQ=255)
                                                                    cGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGCCGGGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaccgccgttttttattaaccactcatttcgactt                                                                                                                                                              >  2:573619‑M1/1‑168 (MQ=255)
                                                                               atatCCATGAAGCGTATCGGTTTAGCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttac                                                                                                                                                   <  1:84776‑M1/251‑95 (MQ=255)
                                                                                   ccATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGATCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataactccttagcaggtatggaaaagcaaacggccattcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttccgccg                                                                                                                                               >  2:14763‑M1/1‑153 (MQ=255)
                                                                                       gAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGGTTGCCGCCACCGCTACCCTAATCATTTTCGGAGACGTGTTgatcaatcactagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagggttaagccgcttt                                                                                                                                           >  2:388504‑M1/1‑149 (MQ=255)
                                                                                        aaGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgctttt                                                                                                                                          <  1:337669‑M1/251‑104 (MQ=255)
                                                                                        aaGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgctttt                                                                                                                                          <  2:452788‑M1/251‑104 (MQ=255)
                                                                                          gCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCTCCATCGAAGATGGTTATGTCCATAACAAATATCAACTTCTGCGGCAAAAGGAGTGTGTCGCCACCGCTACCCTTCTCATTTTCGTGGACGTGTTgataaatcctttgcaggtatggaaaagcaaactgagcttcacattattaaacgccgttttttattaacaactcatttctactttatagcgttacgccgcttttga                                                                                                                                        <  2:205130‑M1/251‑106 (MQ=255)
                                                                                          gCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgcctttttttattaacaactcatttcgactttatagcgttacgccgcttttga                                                                                                                                        <  1:322614‑M1/251‑106 (MQ=255)
                                                                                                  gTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgcc                                                                                                                                <  1:413309‑M1/251‑114 (MQ=255)
                                                                                                          gCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcg                                                                                                                          <  2:38044‑M1/249‑120 (MQ=255)
                                                                                                          gCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCAACGCTACCCTACTCATTTTCGGAGAAGTGTTtataaatccttagcaggtatggcaaagcaaccggcgcttcacattatgaaacgccgttttttatttacaactcatttcgactttatagcgttacgccgcttttgaagttcgccagttcgcg                                                                                                                        >  2:310011‑M1/1‑130 (MQ=255)
                                                                                                               aaCGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCCTTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgcagcttttgaagatcgccagttcgcggaagt                                                                                                                   >  1:546094‑M1/1‑125 (MQ=255)
                                                                                                                                ggTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatcgcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttaccg                                                                                                  <  2:534480‑M1/251‑144 (MQ=255)
                                                                                                                                      cAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgc                                                                                            <  1:461568‑M1/251‑150 (MQ=255)
                                                                                                                                         aGTCACCATCGAAGATGGTTATTTCAATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttagactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgctt                                                                                          >  1:215430‑M1/1‑99 (MQ=255)
                                                                                                                                             aCCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTAATTTTCGGAGACGTGTTgctaaatccttagcaggactggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcatttgaagatcgccagttcgcggaagtcgttacgcacgttacaggattgcttaccgt                                                                                     >  1:579580‑M1/1‑95 (MQ=255)
                                                                                                                                                   gAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaa                                                                               <  1:579787‑M1/251‑163 (MQ=255)
                                                                                                                                                               ttCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatgg                                                                   >  1:91767‑M1/1‑77 (MQ=255)
                                                                                                                                                                     aacaaataacaaCTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcg                                                              <  1:427590‑M1/250‑180 (MQ=255)
                                                                                                                                                                      acaaataacaaCTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcgat                                                            <  1:575811‑M1/251‑182 (MQ=255)
                                                                                                                                                                       caaataacaaCTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcgat                                                            >  1:72288‑M1/1‑69 (MQ=255)
                                                                                                                                                                       caaataacaaCTACGGCGGCAAAAAGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggggtcgat                                                            >  1:460526‑M1/1‑69 (MQ=255)
                                                                                                                                                                        aactaacaactaCGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTggtaaatccttagcaggtatggaaaagcaaacggagcttcaaattatgaaacgccggtttttattaacaactcctttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgtaacaggtttgcttaccgttggcaaactcaaagatggggtcgata                                                           >  2:557266‑M1/4‑68 (MQ=255)
                                                                                                                                                                        aaataacaaCTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcgataa                                                          <  1:61211‑M1/251‑184 (MQ=255)
                                                                                                                                                                         aataacaaCTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcgataa                                                          >  1:401074‑M1/1‑67 (MQ=255)
                                                                                                                                                                                       cggcAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcgataaattcttccagcaact                                           <  1:241022‑M1/251‑199 (MQ=255)
                                                                                                                                                                                             aaGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatc                                                                                                                                   >  2:213913‑M1/1‑47 (MQ=255)
                                                                                                                                                                                             aaGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatc                                                                                                                                   <  1:213913‑M1/157‑111 (MQ=255)
                                                                                                                                                                                              aGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcgataaattcttccagcaactgcggcat                                    >  1:489252‑M1/1‑46 (MQ=255)
                                                                                                                                                                                                       ccgccACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcgataaattcttccagcaactgcggcatcgctttacc                           >  1:577352‑M1/1‑37 (MQ=255)
                                                                                                                                                                                                        cgccACCGCTACCCTACTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcgataaattcttccagcaactgcggcatcgctttaccg                          <  1:225684‑M1/251‑216 (MQ=255)
                                                                                                                                                                                                          ccACCGCTACCCTTCTCATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgctgaagtcgttacgctcgttacatgttggcttaccgttggcaaactcaacgatggtgtcgataaattcttccagcaactgcggcatcgctttaccgtg                        <  1:559034‑M1/251‑218 (MQ=255)
                                                                                                                                                                                                                  accatactcATTTTCGGAGACGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgctggtgtcgataaattcttccagcaactgcggcatcgctttaccgtggatcaact                >  2:156170‑M1/5‑26 (MQ=255)
                                                                                                                                                                                                                    actACTCATTTTCGGAGACATGTTgatcaatccttagcaggtatggaatagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatagccagttcgcggcagtcgttacgctcgttacaggtttgcttcccgttggcaaactcaacgatggtgtagataacttcttccagcaactgcggcatcgctttaccgtggatcacctg               >  2:91179‑M1/2‑24 (MQ=255)
                                                                                                                                                                                                                                cGGAGACGTTTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacatctcatttcgactttatagcgttacgccgcgtttgaagatctccagttcgcggaagtcgttacgctcgtttcaggtttgcttaccgttggcaaactcaacgatggtgtcgataaattcttccagcaacggcggcatcgctttaccgtggatcaactgacccgcgtcaaag  <  2:41208‑M1/251‑240 (MQ=255)
                                                                                                                                                                                                                                  gagaCGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcgataaattcttccagcaactgcggcatcgctttaccgtggat                     <  1:95980‑M1/230‑221 (MQ=255)
                                                                                                                                                                                                                                  gagaCGTGTTgataaatccttagcaggtatggaaaagcaaacggcgcttcacattatgaaacgccgttttttattaacaactcatttcgactttatagcgttacgccgcttttgaagatcgccagttcgcggaagtcgttacgctcgttacaggtttgcttaccgttggcaaactcaacgatggtgtcgataaattcttccagcaactgcggcatcgctttaccgtggat                     >  2:95980‑M1/1‑10 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
CGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGATGTGTTATGTCTGATACCAAACGTAATACAATCGGCAAATTCGGCTTGCTCTCGCTGACTTTTGCCGCCGTTTACAGCTTTAACAACGTTATCAATAATAATATTGAGCTTGGACTGGCCTCGGCACCGATGTTTTTCCTCGCGACGATTTTTTATTTTATTCCCTTCTGTCTGATCATCGCAGAATTTGTTTCGTTAAATAAAAACTCAGAAGCCGGTGTCTACGCGTGGGTAAAAAGTTCGCT  >  NC_000913/3225998‑3226472

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