Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,341 C→G V335V (GTC→GTG)  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,3410CG98.3% 160.9 / ‑3.4 59V335V (GTC→GTG) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (26/32);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.60e-01

AAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTC  >  NC_000913/1094‑1574
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          
aaGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCCCCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTTCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGttt                                                                                                                                                                                                                                         <  1:38423/251‑1 (MQ=255)
            tttCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGAATGGCGGCGCGCGGGTTTGCAGCGCTGTc                                                                                                                                                                                                                              >  2:360153/1‑250 (MQ=255)
                          gcTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTc                                                                                                                                                                                                               >  1:462611/1‑251 (MQ=255)
                             aaaGTTCTTCGCCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTcc                                                                                                                                                                                                              <  1:482177/249‑1 (MQ=255)
                                          cccGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACg                                                                                                                                                                                              <  1:37887/251‑1 (MQ=255)
                                             gCACCATTACCCCCATCGCCCAGTTCCAGAGCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAAGCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCaa                                                                                                                                                                                           <  1:244473/251‑1 (MQ=255)
                                               aCCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTCTTTCCGT‑GGTGCTGATTACGCAatc                                                                                                                                                                                         >  1:244979/1‑251 (MQ=255)
                                                      cccccATCGCCCAGTTCCCGATCCCTTGCCTGATTAAAAATACCGGGAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCATTGTTCAGCGTTTCCGGCCCGGGGGTGAAAGGGATGGTTGGCATGGCGGCACGCGCGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTAAGCAATAATCTTcc                                                                                                                                                                                  >  1:346563/1‑251 (MQ=255)
                                                      cccccAGCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTcc                                                                                                                                                                                  <  1:360153/251‑1 (MQ=255)
                                                             cGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACTCAATCATCTTCCGAATAc                                                                                                                                                                            <  1:116377/250‑1 (MQ=255)
                                                               cccAGGTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTAGTGCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAgt                                                                                                                                                                         <  2:196053/251‑2 (MQ=255)
                                                                ccAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCACGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGACAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAgt                                                                                                                                                                         >  1:458941/1‑249 (MQ=255)
                                                                       cAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGttt                                                                                                                                                                 <  1:482586/251‑1 (MQ=255)
                                                                        aGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTc                                                                                                                                                                >  1:430629/1‑251 (MQ=255)
                                                                          aTCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTg                                                                                                                                                              <  2:255160/251‑1 (MQ=255)
                                                                            cccTTGCCTGATTAAAAAGACCGGAAACCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCg                                                                                                                                                            <  2:3839/251‑1 (MQ=255)
                                                                                      aTTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAg                                                                                                                                                  <  2:430629/251‑1 (MQ=255)
                                                                                            actACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCCATCATCTTCCGAATACAGTATCCGTTTCTGCGTTCCGCaaaacggctg                                                                                                                                            >  2:171764/3‑244 (MQ=255)
                                                                                             agtccGGAAATCCTCAAGCTCCAGGTACGCTCATTGGGGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACtgt                                                                                                                                           <  2:25648/248‑1 (MQ=255)
                                                                                                ccGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTTAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCTTTCCGCAAAGCGACtgtgtg                                                                                                                                        <  2:227230/251‑1 (MQ=255)
                                                                                                ccGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGGGTTTGCAGCGGTGTCACGCGCCCCGATTTTCGT‑GGTGCTGATTACGCGATCATCTTCCGAAAACAGTATCAGTTTCTGCGTTTCGCAAAGCGACagtgtg                                                                                                                                        >  2:195379/1‑251 (MQ=255)
                                                                                                 cGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGc                                                                                                                                       <  1:171764/251‑1 (MQ=255)
                                                                                                      aTCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCt                                                                                                                                  <  1:445854/251‑1 (MQ=255)
                                                                                                      aTCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGGGCGAGCt                                                                                                                                  <  2:219384/251‑1 (MQ=255)
                                                                                                                 ccgtgtacgcgcagcggttccatccGTGATGAAGACTAATTACCGGGCAAGTGCATGTCCATTCTGAATAAGATGGCAATGTGCCGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTTCAGCGATGTCACTCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCTGAATTCAGTATAAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAAt                                                                                                                       <  2:437118/228‑1 (MQ=255)
                                                                                                                 ccAGGTACGCTCATTGGTGCCAGCCGTGATGAATACGAATGACCGGTCAAGGGCATTTCCAATCTGAATAAGAGGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGTCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAAt                                                                                                                       <  2:462611/251‑1 (MQ=255)
                                                                                                                 ccAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTCTTTCCGT‑GGTGCTGATTACGCCATCATCTTCCGAATACAGAATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAAt                                                                                                                       >  2:170842/1‑251 (MQ=255)
                                                                                                                   aGGTACGCTCATTGGAGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTATGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGc                                                                                                                     >  1:303970/1‑251 (MQ=255)
                                                                                                                     gTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAg                                                                                                                   <  1:453856/251‑1 (MQ=255)
                                                                                                                       aCGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGGTGTCACGCGCCCGTATTTCCGT‑GGGGCTGATTACGCAATCATCTTCCGGATACAGGATCAGTTTCTGc                                                                                                                                                             >  2:117780/1‑207 (MQ=255)
                                                                                                                       aCGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGc                                                                                                                                                             <  1:117780/207‑1 (MQ=255)
                                                                                                                                   gCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCt                                                                                                     <  1:244995/251‑1 (MQ=255)
                                                                                                                                    ccAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCCGGAAGAGTTCTACCTg                                                                                                    >  2:183381/1‑251 (MQ=255)
                                                                                                                                          gTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGATCTACCTGGAACTg                                                                                              >  1:441056/1‑251 (MQ=255)
                                                                                                                                             aggaAGCCGAATTAACGGTCAAGGGCATTTCCAATCTGAATAATATGGGAATGTTTAGCGTTTACGGCACGGGGATGAAAGGGATGGTGGGCATGGGGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTGTCGGCGTTCCGCAAAGCGACGGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCGACCTGGAACTGaaa                                                                                           <  2:110819/249‑1 (MQ=255)
                                                                                                                                              tgaAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGCTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCCACTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCCACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGaa                                                                                            >  1:66902/1‑249 (MQ=255)
                                                                                                                                                     gAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGGTTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGGGCTGAACGGGAAATGACGGAAGGGTTCTACCTGGAACTGAAAGACGGCtt                                                                                   >  2:258914/1‑251 (MQ=255)
                                                                                                                                                                cAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGGGcc                                                                          >  2:426550/1‑249 (MQ=255)
                                                                                                                                                                cAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGcc                                                                          <  1:426550/249‑1 (MQ=255)
                                                                                                                                                                 aaGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGGGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGCATGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTCCTGGAGCCGCt                                                                       >  1:304401/1‑251 (MQ=255)
                                                                                                                                                                 aaGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGGATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGGGCTGAACGGGCAATGCAGGAAGATGTCTACCTGGGACTGAAAGAAGGCTTACATGAGCCGCt                                                                       >  2:481842/1‑251 (MQ=255)
                                                                                                                                                                    ggCATTTCCAATCTGGATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGc                                                                    <  2:274454/251‑1 (MQ=255)
                                                                                                                                                                            cAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGG‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGGTTATGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGGGTTCTACCTGGAACTGACAGAAGACTTACTGGAGCCGCTGGcggttaccg                                                            >  1:260793/1‑243 (MQ=255)
                                                                                                                                                                                     taataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGTTTACGCAATCATCTTCCGCATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGaacgaaccgccgg                                                   >  1:209783/5‑239 (MQ=255)
                                                                                                                                                                                        taTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTggac                                                 >  1:400077/2‑248 (MQ=255)
                                                                                                                                                                                         aTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGTATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTGCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCAt                                               <  1:305385/251‑1 (MQ=255)
                                                                                                                                                                                              aaTGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTGGCCCTTAtct                                          >  1:178941/1‑251 (MQ=255)
                                                                                                                                                                                               atggtCAGCGTGTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCAGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTAGCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCGGGCGGTGACGGTTCTCCTGGCCATTAtctc                                         <  1:197310/247‑1 (MQ=255)
                                                                                                                                                                                               aggTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTAtctc                                         <  1:234003/249‑1 (MQ=255)
                                                                                                                                                                                                 gTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACCCGCCCGTCTTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGAGACTGTGTGCGAGCTGAACGGGCAAGGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGTCGGAAAGGCGGGCCACTTTCTCgg                                       >  2:499872/1‑251 (MQ=255)
                                                                                                                                                                                                           tccggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTGGCCCTTATCTCggggggaggggg                             >  1:211457/7‑241 (MQ=255)
                                                                                                                                                                                                            ccggcccGGGGATGAAAGGGATGGTTGGCATGGGGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGGGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGa                             <  2:178941/245‑1 (MQ=255)
                                                                                                                                                                                                                cctggggATGAAAGGGATGATTGGCACGGGGGCTCGCGTGTTTGCAGGGATTCCACGCGCCTGTATTTCCGTTGGT‑CTGATTTCGCAATCATCCTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGGGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGATCTGAAAGAAGTCTTACTGGAGCCGCTGGCGGGGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTa                        <  2:363392/248‑1 (MQ=255)
                                                                                                                                                                                                                     ggATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGGGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTGGCCCTTATCCCGGGGGGAGGTGATGGGATgggg                   >  1:114108/1‑248 (MQ=255)
                                                                                                                                                                                                                            aGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCa                                                                                                                    >  1:143058/1‑147 (MQ=255)
                                                                                                                                                                                                                            aGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCa                                                                                                                    <  2:143058/147‑1 (MQ=255)
                                                                                                                                                                                                                            aGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGc            <  1:263026/251‑1 (MQ=255)
                                                                                                                                                                                                                                 tGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGCTGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACCCAATCCTCTTCCGAATACCGTATCAGTTTCTGCGTTCCGCAAAGCGCCTGTGGGCGAGCTGAACGGGCAATGCTGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGAAGGGAAGGCCGGCCATTATCTCGGGGGGAGGGGAGGGGTGGCACACCTTGCGTggg       >  2:155773/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ggCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTGCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCAGTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGAtctc  <  1:23288/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          
AAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTC  >  NC_000913/1094‑1574

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