Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,509 A→G A391A (GCA→GCG thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,5090AG97.3% 224.0 / ‑3.3 73A391A (GCA→GCGthrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (35/36);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

TTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGA  >  NC_000913/1263‑1747
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
ttCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTgg                                                                                                                                                                                                                                                     >  1:552754/1‑250 (MQ=255)
 tCCAATCTGAATAATATGGCAAGGTTCATCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTTCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGa                                                                                                                                                                                                                                                     <  1:394236/250‑1 (MQ=255)
  ccAATCTTAATAATATGTCAATGTGCAGCGTTTCCGGCCCGGGGATTAAAGGGATGGTTTGCATGGGGGCGCGCGTGATTTCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCACTCATCATCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGAAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTTGCGGTGAc                                                                                                                                                                                                                                                    <  2:594500/250‑1 (MQ=255)
    aaTCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGACggg                                                                                                                                                                                                                                                >  1:340465/1‑250 (MQ=255)
           ataataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTg                                                                                                                                                                                                                                         <  1:35666/251‑1 (MQ=255)
           ataataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTg                                                                                                                                                                                                                                         <  2:266950/251‑1 (MQ=255)
              ataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGGATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGcc                                                                                                                                                                                                                                      >  1:251315/3‑251 (MQ=255)
               taTGGCAATGTTCAGCGTTTCCGGCCCGGGGCTGCACGGGCTGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGCTGTCCCGCGCCCGTATTTCCGTGGTGCTGGTTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGGACTTAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACTGGAACgcttgccc                                                                                                                                                                                                                                      >  1:254445/2‑244 (MQ=255)
                  ggCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCAtt                                                                                                                                                                                                                                   <  1:604702/250‑1 (MQ=255)
                    cAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGGCGGAAAGGcggggcactttt                                                                                                                                                                                                                                >  2:237356/1‑240 (MQ=255)
                       tGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGTTTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGGGTGCGAGCTGAACGGGCAATGAAGGAAGAGTTCTACCAGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGGCGGGACGGCTgggcacttatt                                                                                                                                                                                                                               >  2:570280/1‑240 (MQ=255)
                           cAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAGATGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGGACGGCAGGCCATTATCTCgaggg                                                                                                                                                                                                                         >  1:24906/1‑247 (MQ=255)
                            agcgtgtccggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCggtggt                                                                                                                                                                                                                        <  2:302104/239‑1 (MQ=255)
                                   ccggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGAt                                                                                                                                                                                                               <  1:123992/246‑1 (MQ=255)
                                     ggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGGCATTATTTCGGTGGGAGGGGATgg                                                                                                                                                                                                            >  1:374563/4‑251 (MQ=255)
                                       cccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGAGACGGAACGGCTGGCCATTCACTCGGTGGTAGGGGATggg                                                                                                                                                                                                           >  1:198460/2‑249 (MQ=255)
                                        ccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTAt                                                                                                                                                                                                         <  1:3999/251‑1 (MQ=255)
                                        ccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTAt                                                                                                                                                                                                         <  1:323048/251‑1 (MQ=255)
                                        ccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGACTGGACGGCTGGCCCTTATCTCGGTGGTAGGGGATGGTAt                                                                                                                                                                                                         >  1:319608/1‑251 (MQ=255)
                                                gAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTTTCTCGGGGGTAGGTGATGGTATGCGCACCt                                                                                                                                                                                                 >  1:606164/1‑251 (MQ=255)
                                                         gggtggCAGGGCGGCGCGCGTGGTTGCTGAGATATGACGCGTCCGTATTGCCGTGGTGCTGATTGCGCAATCATCTTCCGAATACAGTATCAGGTTCCGCGTTCCGCAAAGCGACTGGGTGCTAGCTGAACGGGCCATGCAGGAGGAGTTCTACTGGGAACTGAAAGAAGGCTTACTGGAGCCGCGGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGa                                                                                                                                                                                        <  2:162740/247‑1 (MQ=255)
                                                             ggCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGAtct                                                                                                                                                                                     <  1:570280/250‑1 (MQ=255)
                                                             ggCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCACGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGAtctc                                                                                                                                                                                    <  1:343346/251‑1 (MQ=255)
                                                                    cggcgCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGaaa                                                                                                                                                                            <  1:196652/251‑1 (MQ=255)
                                                                          gcgTGTTTGCAGCGATGTCACGCGCTCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGTATCAGTTTCTGCGTTCCGCAAAGCGACCGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCttt                                                                                                                                                                      <  1:643932/251‑1 (MQ=255)
                                                                              gtttGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGGACGGATGGCCATTAACACGGGGTGAGGGTGTTGTATGCGCACCATGGGTGGGATTTCGGCGCAAATTTTTGccc                                                                                                                                                                  >  2:268561/2‑250 (MQ=255)
                                                                               tttGCAGCGATGTCACGCGCCCGTATGTCCTTGGTGGTGATTACGCAAGCGTCTACCGAATGCAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCGGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCg                                                                                                                                                                  <  2:164807/250‑1 (MQ=255)
                                                                                             aCGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGGGATGGTATGCGGAACTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGccc                                                                                                                                                   >  1:10537/1‑250 (MQ=255)
                                                                                                     gTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAAc                                                                                                                                          <  1:150850/251‑1 (MQ=255)
                                                                                                         ttCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGGATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTg                                                                                                                                      >  1:93488/1‑251 (MQ=255)
                                                                                                           ccGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGGGATGTGTTGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCACCATTGTc                                                                                                                                    >  2:110947/1‑251 (MQ=255)
                                                                                                             gtggtgCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCAACTTGCGTGGGGTCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTTTCGc                                                                                                                                 >  1:655808/1‑251 (MQ=255)
                                                                                                                         aCGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAggg                                                                                                                    <  1:8296/251‑1 (MQ=255)
                                                                                                                          cGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAggg                                                                                                                    >  1:288402/1‑250 (MQ=255)
                                                                                                                          cGCAATCATCTGCCGAATACGGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGa                                                                                                                   <  2:340465/251‑1 (MQ=255)
                                                                                                                            cAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGGAGGGGATGGTCAGGGCACCCTTCGGGGGATCTCGGCGAAAATCTTTGCCGCCCTGGCCCCCCCCAAGAACCACATTGTCCGCCTTGCTCAggggt                                                                                                                  >  2:290569/1‑248 (MQ=255)
                                                                                                                                 atcTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGGGATGGGATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCCTTGCTCAGGGATCTTCTg                                                                                                            >  1:400836/1‑251 (MQ=255)
                                                                                                                                    ctCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCCACATTGTCGCCATTGCTCAGGGATCTTCTGAAc                                                                                                         >  1:193722/2‑251 (MQ=255)
                                                                                                                                     tCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGa                                                                                                           <  1:352537/248‑1 (MQ=255)
                                                                                                                                      ccGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACg                                                                                                        <  1:274611/250‑1 (MQ=255)
                                                                                                                                      cagaatgcttgatctgTTTCTGCGTTCCGCAAAGCGACTGTGTGCGTGCTGAACGGGCACTGCAGTAAGTGTTCTTCCTGTATCTGACAGAATGCTTCCTGGTGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCGTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATACCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGc                                                                                                       <  2:514732/236‑1 (MQ=255)
                                                                                                                                         aataaagtaTCAGGTTCTGCGTTCGGCAAAGCGACGGTGTGCGAGCGGAACGGGCAATGCAGGAAGAGTTCTACCAGGAACTGAAAGAAGGCTTGCTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCa                                                                                                    <  2:251315/243‑1 (MQ=255)
                                                                                                                                         aaTACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGATCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGGAGGGGAGGGGAAGAGAACCCTGAGTGGGCTATCGGCGAAAATCTTTGCCGCGCTGGCCCGCGACAAAATAAACAATGTCGCACTTGCTCAGGAGACTTCTGAAcccacc                                                                                                    >  2:439439/1‑246 (MQ=255)
                                                                                                                                            tcagtaTCAGATTCTGCGTTCCGCAAAGCGACTGTGTGAGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACAGCTGGCAATTATCTCGGTGGTAGGCGGTAGTAAGCGCACCTTGCGTGGGATATCGAGGAAATTCTTTGCAGGTCTGGCCCGCCAACATTTCAACCTTTTCGCCCTAGCACAGGGATCTTCTTAACGcccaca                                                                                                  >  2:510935/6‑245 (MQ=255)
                                                                                                                                            aCAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAAtc                                                                                                 <  1:237356/251‑1 (MQ=255)
                                                                                                                                             cagtatccgTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGGTCTTCTGAACGCTCAAtct                                                                                                >  1:240711/9‑251 (MQ=255)
                                                                                                                                              agtaTCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGGGGGAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCCATATCCACATTGTCGCCATTGCTCCGGGATCTTCTGAACGCTCAAtctc                                                                                               >  1:375278/4‑251 (MQ=255)
                                                                                                                                               gtaTCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGACGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGGACGGCTGGCCATTATCTCGGTGGTAGGCGATGGTATGCGCACCTTTCGGGGGGTCCCCGCGGAAATCTTTGGCGCGCTGGCCCCCGCCAATATCAACGTTTTCGCCCTTGCTCCGGGGTCTTCTggctgaccaatctct                                                                                              >  2:202015/3‑237 (MQ=255)
                                                                                                                                                      tttCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGGGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTgg                                                                                       >  1:420413/1‑251 (MQ=255)
                                                                                                                                                       ttCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGACGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGCGATGGGATGCGCACACTGCGTAGGATCTCCGCGCAATTTTTTGCCCCGCAGGCCCGCGCCCATATCAACATTTCCGCCCTTGCCCAGGGACCTTCTGAACGCACACTCTCCGTCGTggg                                                                                      >  2:253589/1‑250 (MQ=255)
                                                                                                                                                           gCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTaaa                                                                                  <  1:202015/250‑1 (MQ=255)
                                                                                                                                                                 cgcaAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAATAGTTCTACCTGGAACTGAAAGAAGGCTTACGGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                                           <  1:233340/249‑1 (MQ=255)
                                                                                                                                                                        cGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGGCGGAAAGGCTGGCCATTATCTCGGTGGTAGGTGATGGGATGCGGACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTTTCGCCATTTCTCCGGGATCTTCTGAACGCTCAATATCTGTCGTGGTAAAAAACGATGATGCg                                                                     >  1:356202/1‑250 (MQ=255)
                                                                                                                                                                         gACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCATGGATCTTCTGAACGCTCAATCTCTGTCGTGGGTAAAAACGATGATGCGAc                                                                   >  1:109989/1‑251 (MQ=255)
                                                                                                                                                                         gACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGGGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGAc                                                                   >  1:348264/1‑251 (MQ=255)
                                                                                                                                                                          aCTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGAcc                                                                  <  1:357840/251‑1 (MQ=255)
                                                                                                                                                                                ggCGAGCTGAACGGGCAATGCAGGAAGATTTCTACCTGGAACGGAAAGAAGTCTGACTGGAGCCGCGGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCGCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGc                                                            <  2:422321/250‑1 (MQ=255)
                                                                                                                                                                                 gCGAGCTGAACGGGCAATGCAGGAAGAGTGCTACCTGGAACTGAAAGAAGGCTTACGGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGAGGATGCGACCACTGGc                                                            <  2:288402/250‑1 (MQ=255)
                                                                                                                                                                                   gAGCTGAACTGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGGAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCTATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTg                                                         <  1:96304/251‑1 (MQ=255)
                                                                                                                                                                                             ggCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGGTGGTATGCGCACCGTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCCTTGCTCAGGGATCTTCTGAACGCCCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTCGCGTGCGCGTTACtt                                               >  1:616517/1‑250 (MQ=255)
                                                                                                                                                                                                  ggCAGGAAGAGTTCTACCTTGAACGGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCACTATCTCGGTGGTAGGTGATGGTATGCGCTCCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACtcatca                                           <  2:10537/249‑1 (MQ=255)
                                                                                                                                                                                                           agTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGGGATGGGATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTggtc                                 >  1:510813/1‑248 (MQ=255)
                                                                                                                                                                                                               atcccTGGAACTGACAGAAGGCTTATTGGAGGCGGTGGTGGGGACGGAACGGGTGGGCTTTATCTCGGTGGTAGGGGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAAGTCTTTGCCGCGCTGGCCCGCGCCAATATCAACAGTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTTATCAGATGTTGTTCAAt                              <  2:606164/247‑1 (MQ=255)
                                                                                                                                                                                                                tACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGGGGTAGGGGGTGGTATGCGCAACTTGCGGGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGGTCTTTTGAACCCACCATCCCTGTCGGGGTAAATAACGCTAATGCGACCACTGGCGTGCGCGGTTCtcatccaggtcttgtcaatac                            >  2:156721/1‑235 (MQ=255)
                                                                                                                                                                                                                 aCCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCGCGGTGGTAGGTGATGGTATTCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATAc                            <  1:428889/250‑1 (MQ=255)
                                                                                                                                                                                                                      gAACTGAAAGAAGGATTACTGGAGCCGCTGACAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGAAGGTATGCGCACCCTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCTCCAATATCAACATTGTAGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTAATCATCAGATGCTGTTCAATACCGATc                       >  1:345450/1‑250 (MQ=255)
                                                                                                                                                                                                                              agaagGCTTCTTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCg              <  2:93488/251‑1 (MQ=255)
                                                                                                                                                                                                                              agaagGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCg              <  1:13668/251‑1 (MQ=255)
                                                                                                                                                                                                                                 agGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGGAGGGGATGGCATGCGCACCTTGGGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGGTCTTCTGGACGCCCAATCTCTGTCGTGGGAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGATGGTCAATACCGATAGGGTTATCGa             >  1:289873/1‑249 (MQ=255)
                                                                                                                                                                                                                                       aCTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGttt      <  1:439439/250‑1 (MQ=255)
                                                                                                                                                                                                                                       aCTGGAGCCGCTGGCGGGGACGGAACGGCTGGCCATTATCTCGGTGGTCGGTGCTGGGATGCGCCCCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCCATATCAACATTTTCGCCCTTGCTCAGGGTTCTTCTGAACGCCCACTCTCTGTCGGTGTAAATAACGCTGATGCGACCAATGGCGGGCGCGTTACTCATCAGGTGCTGTTCAATTCCGGTCCGGTTATCGGAGTttttt     >  2:343543/1‑250 (MQ=255)
                                                                                                                                                                                                                                         gggtgCCGTTGGCGGTGACGGAACGGCGGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTTGGATCTCGTCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGAGCAGGTTATCGAAGTGTTtgtg   <  2:109989/247‑1 (MQ=255)
                                                                                                                                                                                                                                          ggAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGa  <  1:87618/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
TTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGA  >  NC_000913/1263‑1747

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