Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,224,628 C→T L666L (CTG→TTG)  ebgA → evolved beta‑D‑galactosidase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,224,6280CT96.7% 195.4 / ‑1.4 63L666L (CTG→TTG) ebgAevolved beta‑D‑galactosidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (35/24);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

ATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGCTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTG  >  NC_000913/3224386‑3224872
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
aTCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGcc                                                                                                                                                                                                                                              >  2:144770/1‑251 (MQ=255)
   aGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCcgc                                                                                                                                                                                                                                           >  1:382895/1‑251 (MQ=255)
         cGGGACCGGGCCTGAAAGAGTACAAACAGGTTAGCGCTCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCg                                                                                                                                                                                                                                     <  2:283443/251‑1 (MQ=255)
                gggCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTc                                                                                                                                                                                                                               >  1:289798/1‑250 (MQ=255)
                   ccTGAAAGAATACATACATGTTTTCTAGCCGGTAAAAATCCACGCGCGGAATCTGATTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCGCGCAGTGGTGCGCGCCGAAGGTGAAACGCTCGCTACGCAGCAGATTAAACTGCGCGACGATGCGCCGAACAGCGAAGCCCCTTTGCAGATCACGCTGCCGCAGGTGGACGCCCGCGAAGCGTTCCTCatca                                                                                                                                                                                                                           <  2:414382/251‑4 (MQ=255)
                      gAAAGAGTACAAACAGGGTATCGCGCCGGTAAAAATCCGCGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGAGGACTACACCCTGGACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCTCGACGCCGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTGGCAGAGCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTGCCTCAACAt                                                                                                                                                                                                                          <  2:187813/249‑1 (MQ=255)
                       aaaGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTAc                                                                                                                                                                                                                       >  2:404858/1‑251 (MQ=255)
                         agagTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACgg                                                                                                                                                                                                                     >  1:238935/1‑251 (MQ=255)
                           aggaCAAACAGGTTAGCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGGTCCTCAACATTACGGTg                                                                                                                                                                                                                   <  2:58242/248‑1 (MQ=255)
                            gTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGa                                                                                                                                                                                                                  <  2:238935/251‑1 (MQ=255)
                               cAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGGGAccc                                                                                                                                                                                                               >  2:256794/1‑250 (MQ=255)
                                 aaCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACGACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCaaa                                                                                                                                                                                                             <  2:235816/251‑1 (MQ=255)
                                         aTCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCg                                                                                                                                                                                                     >  1:226016/1‑251 (MQ=255)
                                           cgcgCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGCCCAAAGATTCCCGcc                                                                                                                                                                                                   >  2:12445/1‑250 (MQ=255)
                                               cctgTAAAAATCCACGCGAGTGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCGGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGACCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCg                                                                                                                                                                                               <  2:131198/248‑1 (MQ=255)
                                                 ggaaaaaaTCCTGGCGCGGGATCTGACTCGCGGCGTGTTGAAAGACGAAAATAAACTGTCGTTTACCACGCTTGATGACTACACCCTTCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACCGCGAAGCCCCCTCGCAGCTCACGCTGCCCCAGTCGGATGCCCGCGAAGCGGTCCTCAACATTACGGTGACCAAAGTTTCCCTCACCCGCt                                                                                                                                                                                             >  2:28034/4‑251 (MQ=255)
                                                        aTCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGCAAATAAACTGTGGTTTACCAAGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCACCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCg                                                                                                                                                                                       >  1:28312/1‑250 (MQ=255)
                                                           cACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGc                                                                                                                                                                                   >  1:126858/1‑251 (MQ=255)
                                                               cgcgGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCgg                                                                                                                                                                                >  1:234862/1‑250 (MQ=255)
                                                                            ccgcgtctaGTTGAAAGTCGAAAATAAACTGTGGTTTTCCTCGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCa                                                                                                                                                                  <  2:219243/243‑1 (MQ=255)
                                                                                 gcgAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTAt                                                                                                                                                             >  2:330564/1‑251 (MQ=255)
                                                                                   gAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCAAGCAGAGGTGCGCGCAGAAGGTGAAACGCTCGCGACGCAGCAGATTAAAATGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACATATc                                                                                                                                                            >  1:27634/1‑250 (MQ=255)
                                                                                      ttGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGtt                                                                                                                                                        >  1:315590/1‑251 (MQ=255)
                                                                                      ttGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGtt                                                                                                                                                        >  2:80653/1‑251 (MQ=255)
                                                                                           aGTCGAAAATAAACTGTTGGTTACCACGCTTGATGACTACTCCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTGCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGc                                                                                                                                                   <  1:58556/251‑1 (MQ=255)
                                                                                           aGTCGAAAATAAACTGTGGTTTACCACACTTGATGACTAGACCCTGCACGCAGAGGGGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCACGCGAAGCGTTCCTCAACATTACGGTGACCAAAGCTTCCCGCACCCGCTACAACGAAGCCGGACTCCCTATCGCCACTTATCAGTTCCCGc                                                                                                                                                   >  2:241883/1‑251 (MQ=255)
                                                                                              cGAAAATAAACTGTGGTTTACCACGCTCGATGACTACACCCTGCACACAGCGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTg                                                                                                                                                 >  2:166483/1‑250 (MQ=255)
                                                                                                 aaaTAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAgg                                                                                                                                             <  1:404858/251‑1 (MQ=255)
                                                                                                          gtgGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACcg                                                                                                                                    >  1:108424/1‑251 (MQ=255)
                                                                                                                tACCACGCTTGATGACTACACCATGAACGCAGAGGTGCGCGCAGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACACGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGACGTAAAACACCGCGCAGc                                                                                                                              >  2:8310/1‑251 (MQ=255)
                                                                                                                   cACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGcc                                                                                                                             >  1:228669/1‑249 (MQ=255)
                                                                                                                    aCGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGCAGGAAAACACCGCGCAGCCAGt                                                                                                                          >  1:384483/1‑251 (MQ=255)
                                                                                                                      gCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGc                                                                                                                        <  1:162859/251‑1 (MQ=255)
                                                                                                                         tgatgaCTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCtt                                                                                                                     >  2:306947/1‑251 (MQ=255)
                                                                                                                           atgaCTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTc                                                                                                                   >  2:83213/1‑251 (MQ=255)
                                                                                                                           atgaCTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGGTGTACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTc                                                                                                                   <  2:384483/251‑1 (MQ=255)
                                                                                                                              aCTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGc                                                                                                                 >  2:4794/1‑250 (MQ=255)
                                                                                                                                tACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCAcc                                                                                                              >  2:356590/1‑251 (MQ=255)
                                                                                                                                       tGCACGCAGAGGTGCGCGCCGAAGGGGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGCTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAaataat                                                                                                       >  2:404202/1‑251 (MQ=255)
                                                                                                                                           cGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATgcgc                                                                                                   <  1:275190/251‑1 (MQ=255)
                                                                                                                                                     cgcgcCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGAc                                                                                         <  2:346771/251‑1 (MQ=255)
                                                                                                                                                     cgcgcCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGAc                                                                                         >  2:332084/1‑251 (MQ=255)
                                                                                                                                                                   aaCGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTTTCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCCTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGCCGATc                                                                           >  2:5261/1‑251 (MQ=255)
                                                                                                                                                                      gCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGtt                                                                        <  2:126858/251‑1 (MQ=255)
                                                                                                                                                                                    agATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGtt                                                          <  1:120144/251‑1 (MQ=255)
                                                                                                                                                                                       ttAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACGCCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTcgc                                                       >  2:302209/1‑251 (MQ=255)
                                                                                                                                                                                       ttAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCACAGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGCAGACGATCGTTTGAGCTGCCCCGTTcgc                                                       >  2:268813/1‑251 (MQ=255)
                                                                                                                                                                                          aaCTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGc                                                    <  1:306947/251‑1 (MQ=255)
                                                                                                                                                                                             tGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTAc                                                 >  2:38310/1‑251 (MQ=255)
                                                                                                                                                                                                gcgACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCGCGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAAc                                              <  2:289798/251‑1 (MQ=255)
                                                                                                                                                                                                gcgACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAAc                                              <  1:330564/251‑1 (MQ=255)
                                                                                                                                                                                                   aCGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTc                                           <  1:345071/251‑1 (MQ=255)
                                                                                                                                                                                                   aCGTGGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTAGCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTc                                           <  1:256794/251‑1 (MQ=255)
                                                                                                                                                                                                        gcgcCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGAt                                      >  2:409147/1‑251 (MQ=255)
                                                                                                                                                                                                          gcCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCa                                    <  1:147003/251‑1 (MQ=255)
                                                                                                                                                                                                                cAGCGAAGCCCCCGTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATGCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTtct                              <  2:349343/251‑1 (MQ=255)
                                                                                                                                                                                                                  gCGAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGAAGCCCGCGAAGAGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGAGGCTACAACTTCGCGCTCACCTtctc                             >  1:381685/1‑250 (MQ=255)
                                                                                                                                                                                                                    gAAGCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTACTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGAGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCaa                           >  1:58106/1‑250 (MQ=255)
                                                                                                                                                                                                                       gCCCCCTTGCAGATCACGCTGCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAAt                       >  1:175066/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ccACGCTGCCTGAGTTGGACCCCCGCGAAGCTTTCCTCAACATTACGGTTACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACCCCCTATCGCCACTTATCATTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCTCCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCg          <  2:234862/250‑1 (MQ=255)
                                                                                                                                                                                                                                        gcggCCGCAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAATATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTGCGCACCAAAGAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTGCTCAAAAATGATTGGCAAACCGACAt      <  2:391923/248‑1 (MQ=255)
                                                                                                                                                                                                                                            ccGNAGTTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTg  <  2:371387/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ccGCAGCTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTg  <  1:404202/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
ATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGACTACACCCTGCACGCAGAGGTGCGCGCCGAAGGTGAAACGCTCGCGACGCAGCAGATTAAACTGCGCGACGTTGCGCCGAACAGCGAAGCCCCCTTGCAGATCACGCTGCCGCAGCTGGACGCCCGCGAAGCGTTCCTCAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTG  >  NC_000913/3224386‑3224872

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