Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 501,710 T→C 16.0% E510G (GAG→GGG)  ybaL ← inner membrane putative NAD(P)‑binding transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913501,7100TC16.0% 190.3 / 5.5 81E510G (GAG→GGG) ybaLinner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base T (47/20);  new base C (2/12);  total (49/32)
Fisher's exact test for biased strand distribution p-value = 1.55e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCAGAGGCGAGCAATT  >  NC_000913/501472‑501935
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 
tCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGTCGCGGGCAATAAtct                                                                                                                                                                                                                                   >  1:427896/1‑239 (MQ=255)
               tcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGAtttt                                                                                                                                                                                                                    >  1:752506/1‑239 (MQ=255)
               tcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGGGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAACTCCACAATCCGGAtttt                                                                                                                                                                                                                    >  2:627716/1‑239 (MQ=255)
                         aTGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGc                                                                                                                                                                                                          >  1:616363/1‑239 (MQ=255)
                                    gcgcGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCCCAATATCCGGATTTTTCGCGCGGGCAGATGCccca                                                                                                                                                                                                >  2:240770/1‑238 (MQ=255)
                                         tCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGGTTTTTCGCGCGGGCAGATGCCACAAtctc                                                                                                                                                                                           >  1:196135/1‑238 (MQ=255)
                                             tttctcccTATATTTTTTTATTTTTTTACTCCCGCCCACCCCCCTTTCCCCCCCCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCCCCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCAccc                                                                                                                                                                                       <  2:616363/233‑1 (MQ=255)
                                                    cctaattttttttATATTTTTCCCTCCGACACCCCCCTTCTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCCCGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCAt                                                                                                                                                                                <  1:209159/232‑1 (MQ=255)
                                                       cATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCCGATGCCACAATCTCACCCCCTTCATAAc                                                                                                                                                                             >  1:739462/1‑238 (MQ=255)
                                                       cATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAAcc                                                                                                                                                                            >  1:745391/1‑239 (MQ=255)
                                                            ttttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTgg                                                                                                                                                                       >  1:383900/1‑239 (MQ=255)
                                                              ttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCCCCCGGGCAGATGCAACAATCTAACCCGCTTCATAACCGTTGGGa                                                                                                                                                                     >  2:899230/1‑239 (MQ=255)
                                                                       ttACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCCCAATCTCACCCGCTTCATaa                                                                                                                                                                              >  1:720322/1‑222 (MQ=255)
                                                                               ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAAcc                                                                                                                                                                            >  1:837870/1‑216 (MQ=255)
                                                                                        ccttttccccctgcccccccccccccGCCGGTGGCGTTTCCAGCAGTTCCCGCATGGCTCGGGCGATTTCACGCTCCCCCATCACTACCGTATTCGCACCACGTTCGGTCATATACGCCATTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGca                                                                                                                                             <  1:899230/216‑1 (MQ=255)
                                                                                         gTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAAtc                                                                                                                                                                                             >  1:225068/1‑189 (MQ=255)
                                                                                          tttACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACAt                                                                                                                                          >  1:769593/1‑239 (MQ=255)
                                                                                               cccGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCCATAATCTCAATATCCGGATTTTTCGCTCGGGCAGATGCCACAATCTCACCCGCTTCATAACCCTTGGGAATCGTCAGGTACACCCATTTTGCACATTCCAg                                                                                                                                     >  2:229335/1‑239 (MQ=255)
                                                                                                  gTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCCGact                                                                                                                                  >  1:236851/1‑237 (MQ=255)
                                                                                                  gTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAAATg                                                                                                                                  >  2:93062/1‑239 (MQ=255)
                                                                                                                  gccgGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTCGGAATCGTCAGGACCCGCCATTTTGCACATTCCAAATGCCCAAGttaccaaat                                                                                                                   >  2:628270/1‑232 (MQ=255)
                                                                                                                              tCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACACGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCCTTTTGCACATTCCAGATTCGCCAGTTGCATAATTTCttcgttcc                                                                                                        >  1:207816/1‑237 (MQ=255)
                                                                                                                                  ttttaaagaccagggTCCGGGCGATTTCCCGCCCCCCCATCACCACCTGATTGCCCCCACGTTCGGTGATATCCGCCACTTCATCGTCATAATGGGCGCGGGCAATAACCTCAAATTCCGGATTTTTCGCGCGGGAAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCa                                                                                                   <  2:928684/226‑1 (MQ=255)
                                                                                                                                  tcattctccccagggtgggggCGTTTTCACCGTTGCCCATCAATACCCTATACGCCCCACGTTCGGTGATATACGCCTATTCACCGCCATAATGGGCGCGGGCAATAATCTCAATTCCCGGTTTTTTCGCGCGGGCAGAGGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCa                                                                                                   <  2:126555/221‑1 (MQ=255)
                                                                                                                                  gctcttCCAGCCTGGCACGGGCAGTTGCCCGCTCCGCCCACTCTTCCTGATTCGCCCCAAGTTCGGGCATATACGCCCCTTCATCGTAATAATGGGCGCGGGAAATAATCCCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCa                                                                                                   <  2:1013359/235‑1 (MQ=255)
                                                                                                                                  gcatttcctggctGGTTCGGGCGATTTCCCGCCCCCCCCTCCCTCCCTTATTCGCCCCACGTTCGGTGATATACGCCCCTTTATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCa                                                                                                   <  2:156737/227‑1 (MQ=255)
                                                                                                                                  gcattgccaccaTGGGCCGGGGGATTCCCCGCCCGCCCATAACTACCTGATTCGCACCACTTTCGGTGATATACGCCACTTCATCGTCATAAGGGGCGCGGGCAATAATCCCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCa                                                                                                   <  2:616092/233‑1 (MQ=255)
                                                                                                                                  gcagTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCa                                                                                                   <  1:93062/239‑1 (MQ=255)
                                                                                                                                  cggtttCCCGCTTGGTCGGGGCGTTTTCCCGCCCGCCCTTCACTACCGGATTCGCCCCACGTTGGGTAATAACCGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCa                                                                                                   <  2:530040/235‑1 (MQ=255)
                                                                                                                                  cgctttccccctgggTACGGGCAGTTCTCACGCCGCCCATCTTCACCTGTTCCCCCACCCGTTGGGGCAAATAAGCCCTTTCATCTTCAAACTGGGCGCGGCCAATAATATCAATATCCGGTTTTTTGGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCa                                                                                                   <  2:771272/226‑1 (MQ=255)
                                                                                                                                  cctgttccccatggggcggggTGATTTCACTCCCCGCCATCCCTACCGGGTTCGCCCAGCGTTGGGGGATTTACGCCCTTTCATCGGCATAATGGGGGCGGGCAATAATCCCAATACCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCa                                                                                                   <  2:958791/221‑1 (MQ=255)
                                                                                                                                  ccctttCCACCCTGGTCGGGGCGTTTTCACCCCCGCCCATCCCTCCCTGATTCCCCCCCCGTTCGGTTATATACCCCCCTTCATCGTCATAATGGGCGCGGGCAATAATCCCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCa                                                                                                   <  2:1080074/235‑1 (MQ=255)
                                                                                                                                  acattccccgctgggTCCGGGCTTTTTCACCGCCCCCCATCCCTTCCGTATTCCCACCCCGTTCGGTATAATACGCCCCTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCa                                                                                                   <  2:196135/226‑1 (MQ=255)
                                                                                                                                       tCCACCCTGGTACGGGCGTTTTCCCGCTCGCCCATCACTCCCTGATTCGCACCACGTTCGGTGATATACGCCACTTCCTCGGCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGcc                                                                                              <  2:752506/239‑1 (MQ=255)
                                                                                                                                        ccccagctgtttcggcggttttCCCGCTCGCCCTCTCCTACCTTATTCCCACCACCGTCGGTAATATACGCCTCTTCACCGTCATAATGGGCGCGGGCAATAATCCCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGcc                                                                                              <  2:917260/220‑1 (MQ=255)
                                                                                                                                             aTGGCACGGGCGTTTTCACGCCGGCCCATCACTTCCTGATTCTCTCCACGTTCGGTGATATACGCCACTTCAGCGTCATAATGGGCGCGGGCAATAATCTCAATAGCCGGATTTGTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATa                                                                                         <  1:445282/238‑1 (MQ=255)
                                                                                                                                                gTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTgc                                                                                     >  1:653522/1‑239 (MQ=255)
                                                                                                                                                gTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAAAACTg                                                                                      >  1:751786/1‑238 (MQ=255)
                                                                                                                                                      gCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATCCTCCACggg                                                                                >  1:1059801/1‑237 (MQ=255)
                                                                                                                                                      gCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCCCggc                                                                                >  1:220914/1‑237 (MQ=255)
                                                                                                                                                          aTTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGAcac                                                                             >  1:485689/1‑236 (MQ=255)
                                                                                                                                                           tttCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATa                                                                                         >  1:928684/1‑225 (MQ=255)
                                                                                                                                                            ttCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATTGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAAAACTCCGCGGCcccccc                                                                          >  2:110460/1‑237 (MQ=255)
                                                                                                                                                                  ccccccccTTCCCTCCCTGGTTCGCCCCCCGTTCGGTGATAAACGCCATTTCACCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGGGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGctct                                                                     <  2:751786/234‑1 (MQ=255)
                                                                                                                                                                      gCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGc                                                                  >  1:88952/1‑238 (MQ=255)
                                                                                                                                                                       cccATCCCCACCGGTTTCGCACCCCGTTCGGGGATATACGGCCTTTCATCGGAAAAATGGGGGCGGGAAAGAATCTCAAGATCCGGATTTTTCGCGCGGGAAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAg                                                                <  2:915770/239‑1 (MQ=255)
                                                                                                                                                                          aTCCCTACCTGATTCGCACCACGTTCGGTGATATACCCCCCTTCATTGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCtc                                                             <  1:110460/239‑1 (MQ=255)
                                                                                                                                                                           tCACTACCTGATTCGCACCCCGTTCGGGGATATACGCCACTTCCTCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCa                                                                                                   <  2:599863/200‑1 (MQ=255)
                                                                                                                                                                                aCCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCCCGCCCCCGCAGCCCACCaaa                                                       >  2:131832/1‑238 (MQ=255)
                                                                                                                                                                                 ccgtaTTCCCCCCACGTTCGGTTATATCCCCCACTTCATCCTCATAAGGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGGAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAaca                                                      <  1:904142/235‑1 (MQ=255)
                                                                                                                                                                                 ccTGTTTCGCCCCCCTTGCGGTGATTTACGCCACTTCCCCGCCAAAATGGGGCAGGGCAATAACCCCAAAACCCGGATTTTGCGGGCGGGAAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCCCATCAaca                                                      <  2:769593/239‑1 (MQ=255)
                                                                                                                                                                                 ccTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAac                                                       >  1:788461/1‑238 (MQ=255)
                                                                                                                                                                                    gATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATCCTGCGCGGCCCCCGCGCTCTCGCAGCCCATCAACACgg                                                   >  2:904142/1‑239 (MQ=255)
                                                                                                                                                                                     aTTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATa                                                                                         >  1:599863/1‑200 (MQ=255)
                                                                                                                                                                                         gCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCgcgc                                                                        <  1:1001196/213‑1 (MQ=255)
                                                                                                                                                                                         gCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCgcgc                                                                        >  2:1001196/1‑213 (MQ=255)
                                                                                                                                                                                           accacGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACAGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTATTCGTTCGCCGCATGGCCCAATACTGCGCGGACCCCGCGCCCTCCCCGCTCACCAAAACGGGTTCgg                                             >  2:405075/1‑237 (MQ=255)
                                                                                                                                                                                           accacGTTCGGTGATATACGCCACTTCATCGTCATAATAGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGTGGGCAGATGCCACAATTTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCAAGATGCCCCAGTTGCATAATTTCCTCGTTCGACGCATTGCCCCAACCTTCCCGGACCCCGGGCTTCCCCACCTCCTCAACACGGGtcctga                                            >  2:904794/1‑234 (MQ=255)
                                                                                                                                                                                               cGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGACCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTGGCCGCATTGCCCAATACTCCCCGGGCCCCGCCCTCTCGCCGCCCCTCAACACGGGTTCGTGCCg                                         >  2:24236/1‑238 (MQ=255)
                                                                                                                                                                                                  tCGGTGAAATACGCCACTTCCACGTCATAATGGGCTCGGGCAATAATCCCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGtct                                      <  1:1014342/238‑1 (MQ=255)
                                                                                                                                                                                                     gTGATATACGCCACTTCATCGTCATAATGTGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTTGCCGCATTGCCCAATACTGCGCGGACCCCGCGCCCTTGCACCTCACCACAACGGGTTCGTGACGCTTCaa                                   >  1:279166/1‑238 (MQ=255)
                                                                                                                                                                                                              ccctcttttcggtaaaaaTGGGCGCGGGCAAATATCCCCATTCCCGGATTTTTCGGGCGGGCAAAAGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAg                         <  2:1019204/223‑1 (MQ=255)
                                                                                                                                                                                                                cACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCCGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTCCTTGACGTCTCAATCCcccccacgg                       >  2:596587/1‑236 (MQ=255)
                                                                                                                                                                                                                cACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCAccaccacc                        >  1:1019204/1‑236 (MQ=255)
                                                                                                                                                                                                                cACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCACCg                       >  1:188120/1‑239 (MQ=255)
                                                                                                                                                                                                                 ttttCATCGGCAAAAGGGGCGGGGGAAAATACCCAAATATCCGGTTTTTTCGCGCGGGGAGATTCCAAAATCTCCCCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCgg                      <  2:1059801/237‑1 (MQ=255)
                                                                                                                                                                                                                 attttctgttCATAATGGGGGCGGGAAATATACCCAATATCCGGATTTTCCAGGCGGGCAGATGCCACAATCGCAAACTGTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTGCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCgg                      <  2:207816/230‑1 (MQ=255)
                                                                                                                                                                                                                  cTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGca                                                                                                                                             >  1:242182/1‑119 (MQ=255)
                                                                                                                                                                                                                  cTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCa                                                                                                   >  1:115861/1‑161 (MQ=255)
                                                                                                                                                                                                                  cTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAAAACGGGTTCGTGACGTCTCAATCACCACCAGCGGa                     >  1:635471/1‑239 (MQ=255)
                                                                                                                                                                                                                  cTTCATCGTCATAATGGGCGCGGGCAATAATCCCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCa                                                                                                   <  2:115861/161‑1 (MQ=255)
                                                                                                                                                                                                                  cTTCATCGTCAAAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGca                                                                                                                                             <  2:242182/119‑1 (MQ=255)
                                                                                                                                                                                                                   ttCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCCCCGCTCATCAACACGGGTTCCTGACGCCCCAACCCCCACCACCGGaa                    >  2:561503/1‑239 (MQ=255)
                                                                                                                                                                                                                       tCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGTTCATCAACACGGGTTCTGGACGTACCAATCCCCCCCCGCGGaaaacc                >  2:549334/1‑235 (MQ=255)
                                                                                                                                                                                                                        cGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCa               >  1:47932/1‑239 (MQ=255)
                                                                                                                                                                                                                         gTCATAAGGGGCGCGGGCAATAAGCTCAATACCCGGATTTTTCGCGCGGGCAGATGCCACAAGCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCa               <  1:740021/238‑1 (MQ=255)
                                                                                                                                                                                                                                   gcgcGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTAGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCCCCGCTCATCAACACGGGTTCGGGACGGCTCAATCCCCCCCAGCGGACTATCCGCGGCGAGcca    >  2:740021/1‑237 (MQ=255)
                                                                                                                                                                                                                                    cgcgGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCAtt                                                                                                 >  1:510289/1‑145 (MQ=255)
                                                                                                                                                                                                                                    cgcgGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGCCGTCTCAATCACCACCACCGGAATATCAGAGGCGAGCAAt   >  1:800107/1‑239 (MQ=255)
                                                                                                                                                                                                                                    cgcgGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCAAAGGCGGGCAAt   >  1:192664/1‑239 (MQ=255)
                                                                                                                                                                                                                                    cgcgGGCAATAATCCCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCAtt                                                                                                 <  2:510289/145‑1 (MQ=255)
                                                                                                                                                                                                                                      cgGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCAGAGGCGAGCAAtt  >  1:30369/1‑238 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 
TCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCAGAGGCGAGCAATT  >  NC_000913/501472‑501935

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