Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,278 C→T N314N (AAC→AAT 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,2780CT98.3% 179.6 / ‑4.4 60N314N (AAC→AATthrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (34/25);  total (35/25)
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.72e-01

ATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTG  >  NC_000913/1034‑1524
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
ataCCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATgg                                                                                                                                                                                                                                                    <  1:642952/249‑1 (MQ=255)
   ccTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATg                                                                                                                                                                                                                                               >  1:302104/1‑251 (MQ=255)
                  aGGTGCCCGATGAGCGGTTTTTGAAGACGCTGTCCTACCAGGAAGAGCTGGAGCTTTCCTACTTCGGCGCTAAAGTTCTACACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTCCAACTACCGGAAATCCTCACGCTACAGGCACGCTCATTGGTGCCCGCCGTGATGAAGACGACTTACCGGTCAGTGACATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTccggc                                                                                                                                                                                                                                >  2:96304/1‑247 (MQ=255)
                    gTGCCCGTTGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTACGGCGCTAAAGTTCTTCCCCCCCGCACCATTACCCCCCTCGCCCCGTTCCAGATCACTTGCCTGCTTAAAACTACCAGAAATCCTCACGCCCCCGGCACGCTCCTTGGTGCAAGCCGCGATGATGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCCGCGtatccgaccc                                                                                                                                                                                                                              >  2:611862/1‑242 (MQ=255)
                              cGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGaa                                                                                                                                                                                                                     >  1:434527/1‑250 (MQ=255)
                              cGAGGGTGTTGAAGTTGATGTCCTACCAGGAATCGATGGAGCTTTCCTACTTCGGCTCTAAAGTTCTTCGCCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAATCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGaaa                                                                                                                                                                                                                    <  1:268415/251‑1 (MQ=255)
                                     gttgAAGTCGATGTCCTACCATGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCCCCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAAGCCTCAATCTCCAGGTACGCTCATTGTTGCCAGCCGTGATGAAGACGAATTACCGGGCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCATCGTTTCCGGCCCGGGGATGAAAGGGATgg                                                                                                                                                                                                             <  1:656303/251‑1 (MQ=255)
                                         aaGTCGCTGTCCTACCAGGAAGCGATGGCGCTTTCCTACTTCGGCGCTAACGATCTTCACCCCCGCACCATTACACCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAACATACCGGCAATCCTCAAGCTCCAGGAACGCTCATTGTTGCCAGCCGAGATGTAGAAGAATTACCGGTAACGGGCTTTTCCACTCTGAATCATCTAGCATTGTACAGCGTGTCCGGCCCGGGGAGGAACGGGATGGtggg                                                                                                                                                                                                         >  2:260233/1‑248 (MQ=255)
                                                    cTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCCTTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCACGGGGAGGAAAGGGATGGTTGGCATggcggg                                                                                                                                                                                                >  2:3999/1‑248 (MQ=255)
                                                     tACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCGATGTTCAGCGTTTCCGGCCCGGGGATGata                                                                                                                                                                                                                    >  2:165335/1‑226 (MQ=255)
                                                     tACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGaaa                                                                                                                                                                                                                    <  1:165335/228‑1 (MQ=255)
                                                            aaGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGtggtt                                                                                                                                                                                      >  1:24031/1‑247 (MQ=255)
                                                               cGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGtgtt                                                                                                                                                                                       >  1:674828/1‑244 (MQ=255)
                                                               cGATGGAGCTTTCCTACTTCGCCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAAGACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTTAATAATAGGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGtgtt                                                                                                                                                                                       <  2:674828/247‑4 (MQ=255)
                                                                 aTGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGc                                                                                                                                                                                 >  1:266950/1‑251 (MQ=255)
                                                                           ccTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCAcgc                                                                                                                                                                       <  1:357358/251‑1 (MQ=255)
                                                                                      gcTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTc                                                                                                                                                            >  1:164807/1‑251 (MQ=255)
                                                                                               gttaccccccGCACCATTGCCCCCATCGCCCAGTTCCAGATCCCTGGCCTGATTAAAAATTCCGGAAATCCTCAAGTGCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCTATCTGAATAATATGGTAATGTTCAGCGGTTCCGGCCCGGGGATGAAAGGGATGTTTGTCATTGCGGCGCGCGTGTTTGCAGCGATGTCACGCGTCCGTATTTCCGTGGTGc                                                                                                                                                    <  2:679959/245‑1 (MQ=255)
                                                                                                ttcACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGGGCTg                                                                                                                                                  >  1:551203/1‑251 (MQ=255)
                                                                                                  cACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTTGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAAGCTGAATAATATGGCAATGTTCAGCGTTTccggcccg                                                                                                                                                                                                                             <  2:511539/174‑8 (MQ=255)
                                                                                                  cACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTccggcccg                                                                                                                                                                                                                             >  1:511539/1‑167 (MQ=255)
                                                                                                  cACCCCCGCACCATTACCACCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGAt                                                                                                                                                >  1:267495/1‑251 (MQ=255)
                                                                                                     ccccGCCCCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTAc                                                                                                                                             <  2:143342/251‑1 (MQ=255)
                                                                                                     ccccGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGTTGCTGATTAc                                                                                                                                             <  1:601722/251‑1 (MQ=255)
                                                                                                     ccccGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGAGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGGGCTGATTAc                                                                                                                                             >  1:61154/1‑251 (MQ=255)
                                                                                                      cccGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGGGCTGATTAc                                                                                                                                             >  1:143342/1‑250 (MQ=255)
                                                                                                              aTTACCCCCATCGCCCAGTTCCAGACCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGGGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCTCGCGCCCGTATTTCCGTGGTGCTGATTACGCAatcatc                                                                                                                                    <  2:434527/251‑1 (MQ=255)
                                                                                                                   ccccATCGCCCAGTTCCAGATCCCTTTCCTGATTAAAAATACCGAAAACCCTCCAGCTCCAGGTACGCTCATTGGTGCCAGCCGTTATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGAGGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTGCCGTGGTTCTGATTACGCAATCATCTTCCg                                                                                                                               <  2:55842/251‑1 (MQ=255)
                                                                                                                      cctCGCCCAATTCCAGATTCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGTGGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAAt                                                                                                                            <  2:508107/249‑1 (MQ=255)
                                                                                                                         cGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGGCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATAc                                                                                                                          <  2:63041/250‑1 (MQ=255)
                                                                                                                           cccAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGCTGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCCATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGGGCGCGGGTTTGCAGGGATGTCACGCGCCCGTTTTTCAGTGGGGCTGATTACGCatt                                                                                                                                        >  2:4361/1‑232 (MQ=255)
                                                                                                                           cccAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAat                                                                                                                                        <  1:4361/234‑1 (MQ=255)
                                                                                                                                ttCCAGATCCATTGCCTGCATAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTAACGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGAAGGCCCGCGCCCGTGTTTCCGGGGGGCTGATTACGCAATCATCTTTCGAATACAgactca                                                                                                                   >  2:428889/1‑245 (MQ=255)
                                                                                                                                       tCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATTCAGCATCAGTTTCTGc                                                                                                           >  1:422321/1‑251 (MQ=255)
                                                                                                                                          cTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGGGTGTTTGCAGGGATGTCACGGGCCCGTATTTCCGGGGGGGTGATTTAGAGATAATCTTCCGAATACAGTATCAGTTTCTGCggt                                                                                                        >  2:268555/1‑249 (MQ=255)
                                                                                                                                             gCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGGACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCACGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTTAAGCGATGTCACGCGCCCGTATTTCAGTGGTGCTGATTACACAATCATCTTCCTAAACCAGTATCCGTTTCTGCTTTccg                                                                                                     >  1:117153/1‑250 (MQ=255)
                                                                                                                                                    tAAAAATACCGGAAATCCTCAAGCTCCAGGTACGATCATTGGCGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGAGCGCGTGTTTGCAGCGAGGTCACGAGCCCGTATTTCCGTGGAGCTGATTACGCCATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCg                                                                                              >  1:523969/1‑251 (MQ=255)
                                                                                                                                                     aaaaaTACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCg                                                                                              >  1:543366/1‑250 (MQ=255)
                                                                                                                                                        aaTACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGCCtg                                                                                          >  1:241550/1‑251 (MQ=255)
                                                                                                                                                                    ccTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGAGCGCGTGCTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAAACCTCTTCCGAATACAGTATCAGTTTCTACGGTTCGCAAAGCGACTGTGTGCgggccg                                                                               >  2:501096/1‑245 (MQ=255)
                                                                                                                                                                              cAGGTACGCTCATGGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTGCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATg                                                                    <  2:284096/251‑1 (MQ=255)
                                                                                                                                                                               aGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGc                                                                   >  1:321352/1‑251 (MQ=255)
                                                                                                                                                                                            ggtggctgcctgtattaaGACGAACTACCGGGCAAGGGCATTTCCAAGCTGAATAATATGCAAATGTTCAGCGTGTCCTGTCCGGGGATGAAAGGGATGGTTGGCTGGGCGGCGCGCTGGGTTACATCAAGGTCACGCTCCCGTATTACCGTGGTGCTGATTACGCAATCATCTTCCGACTACAGGATCAGTTTCTGCTTTCCGCAAAGCGACTGGCTGCGAGCTGAACGGGCAATGCAGCAagtggtcta                                                      <  2:250280/235‑8 (MQ=255)
                                                                                                                                                                                                    ccgggatgaatgcgcatgacctgTCAAGGGCTTTGCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCTTGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCTCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGGGCGAGCTGAACGGGCAAGGCAGGAAGAGTTCTTCCTGGAAc                                              <  2:551203/229‑1 (MQ=255)
                                                                                                                                                                                                       tgatgaAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGAGTGTTTTCAGCGGTGTCCCGCGCCCGGATTTCCGTGGGGCTGGTTACGCAATCATCTTCCGGATACAGGATATGTTTTTGCGTTCCGCAAAGCGACTGTGTGGGCGCTAGACGGGGAATGCCGGAagggatctcaccgggactg                                            >  2:13668/1‑233 (MQ=255)
                                                                                                                                                                                                                   aTTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCTGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTAc                               <  1:245272/251‑1 (MQ=255)
                                                                                                                                                                                                                      acctgTCAAGGGCATTTCCAAGCAGAATAATATGGCAATGTGCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCGTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTg                             <  2:239904/246‑1 (MQ=255)
                                                                                                                                                                                                                      aCCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGACAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGGGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGACGAGTTCTACCTGGAACTGCaagacagcttcatgg                            >  2:35666/1‑240 (MQ=255)
                                                                                                                                                                                                                      aCCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCCATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTgg                            >  1:324984/1‑251 (MQ=255)
                                                                                                                                                                                                                             aaGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTGCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCt                     <  1:368464/251‑1 (MQ=255)
                                                                                                                                                                                                                               gggCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACggg                                                                         >  1:644821/1‑197 (MQ=255)
                                                                                                                                                                                                                               gggCAGTTCCAAGCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACggg                                                                         <  2:644821/197‑1 (MQ=255)
                                                                                                                                                                                                                                     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)
                                                                                                                                                                                                                                         aaTCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGGTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCa                                                                       >  1:246934/1‑189 (MQ=255)
                                                                                                                                                                                                                                          aTCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGCTGACACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCCATGCAGGAAGAGGTTTAACTGGAACTGAAAGAAGGCTTACAGGGGCCGCTGGcgggggggg          >  2:87618/1‑241 (MQ=255)
                                                                                                                                                                                                                                          aTCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTgg                                                 <  1:643848/210‑1 (MQ=255)
                                                                                                                                                                                                                                          aTCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGGATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAACGACTGTGTGCGAGCTGAACGGGCACTGCAGGAAAAGTTCTACCTgg                                                 >  2:643848/1‑210 (MQ=255)
                                                                                                                                                                                                                                                ataataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTg  <  1:35666/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ataataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTg  <  2:266950/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
ATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTG  >  NC_000913/1034‑1524

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