Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,407 C→T S357S (AGC→AGT)  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,4070CT96.7% 168.4 / ‑2.2 64S357S (AGC→AGT) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (28/31);  total (31/33)
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.37e-01

CAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTG  >  NC_000913/1165‑1642
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             
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‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCTTTCCGCAAAGCGACtgtg‑tg                                                                                                                                                                                                             <  2:227230/251‑1 (MQ=255)
                         ccGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGGGTTTGCAGCGGTGTCACGCGCCCCGATTTTCGT‑GGTGCTGATTACGCGATCATCTTCCGAAAACAGTATCAGTTTCTGCGTTTCGCAAAGCGACagtg‑tg                                                                                                                                                                                                             >  2:195379/1‑251 (MQ=255)
                          cGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGc                                                                                                                                                                                                            <  1:171764/251‑1 (MQ=255)
                               aTCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCt                                                                                                                                                                                                       <  1:445854/251‑1 (MQ=255)
                               aTCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑GGCGAGCt                                                                                                                                                                                                       <  2:219384/251‑1 (MQ=255)
                                          ccgtgtacgcgcagcggttccatccGTGATGAAGACTAATTACCGGGCAAGTGCATGTCCATTCTGAATAAGATGGCAATGTGCCGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTTCAGCGATGTCACTCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCTGAATTCAGTATAAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAAt                                                                                                                                                                                            <  2:437118/228‑1 (MQ=255)
                                          ccAGGTACGCTCATTGGTGCCAGCCGTGATGAATACGAATGACCGGTCAAGGGCATTTCCAATCTGAATAAGAGGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGTCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAAt                                                                                                                                                                                            <  2:462611/251‑1 (MQ=255)
                                          ccAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTCTTTCCGT‑GGTGCTGATTACGCCATCATCTTCCGAATACAGAATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAAt                                                                                                                                                                                            >  2:170842/1‑251 (MQ=255)
                                            aGGTACGCTCATTGGAGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTATGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGc                                                                                                                                                                                          >  1:303970/1‑251 (MQ=255)
                                              gTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAg                                                                                                                                                                                        <  1:453856/251‑1 (MQ=255)
                                                aCGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGGTGTCACGCGCCCGTATTTCCGT‑GGGGCTGATTACGCAATCATCTTCCGGATACAGGATCAGTTTCTGc                                                                                                                                                                                                                                   >  2:117780/1‑207 (MQ=255)
                                                aCGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGc                                                                                                                                                                                                                                   <  1:117780/207‑1 (MQ=255)
                                                            gCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCt                                                                                                                                                                          <  1:244995/251‑1 (MQ=255)
                                                             ccAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCCGGAAGAGTTCTACCTg                                                                                                                                                                         >  2:183381/1‑251 (MQ=255)
                                                                   gTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGATCTACCTGGAACTg                                                                                                                                                                   >  1:441056/1‑251 (MQ=255)
                                                                      aggaAGCCGAATTAACGGTCAAGGGCATTTCCAATCTGAATAATATGGGAATGTTTAGCGTTTACGGCACGGGGATGAAAGGGATGGTGGGCATGGGGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTGTCGGCGTTCCGCAAAGCGACGGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCGACCTGGAACTGaaa                                                                                                                                                                <  2:110819/249‑1 (MQ=255)
                                                                       tgaAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGCTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCCACTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCCACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGaa                                                                                                                                                                 >  1:66902/1‑249 (MQ=255)
                                                                              gAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGGTTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGGGCTGAACGGGAAATGACGGAAGGGTTCTACCTGGAACTGAAAGACGGCtt                                                                                                                                                        >  2:258914/1‑251 (MQ=255)
                                                                                         cAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGGGcc                                                                                                                                               >  2:426550/1‑249 (MQ=255)
                                                                                         cAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGcc                                                                                                                                               <  1:426550/249‑1 (MQ=255)
                                                                                          aaGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGGGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGCATGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTCCTGGAGCCGCt                                                                                                                                            >  1:304401/1‑251 (MQ=255)
                                                                                          aaGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGGATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGGGCTGAACGGGCAATGCAGGAAGATGTCTACCTGGGACTGAAAGAAGGCTTACATGAGCCGCt                                                                                                                                            >  2:481842/1‑251 (MQ=255)
                                                                                             ggCATTTCCAATCTGGATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGc                                                                                                                                         <  2:274454/251‑1 (MQ=255)
                                                                                                     cAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGG‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGGTTATGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGGGTTCTACCTGGAACTGACAGAAGACTTACTGGAGCCGCTGGcggttaccg                                                                                                                                 >  1:260793/1‑243 (MQ=255)
                                                                                                              taataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGTTTACGCAATCATCTTCCGCATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGaacgaaccgccgg                                                                                                                        >  1:209783/5‑239 (MQ=255)
                                                                                                                 taTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTggac                                                                                                                      >  1:400077/2‑248 (MQ=255)
                                                                                                                  aTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGTATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTGCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCAt                                                                                                                    <  1:305385/251‑1 (MQ=255)
                                                                                                                       aaTGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTGGCCCTTAtct                                                                                                               >  1:178941/1‑251 (MQ=255)
                                                                                                                        atggtCAGCGTGTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCAGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTAGCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCGGGCGGTGACGGTTCTCCTGGCCATTAtctc                                                                                                              <  1:197310/247‑1 (MQ=255)
                                                                                                                        aggTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTAtctc                                                                                                              <  1:234003/249‑1 (MQ=255)
                                                                                                                          gTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACCCGCCCGTCTTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGAGACTGTG‑TGCGAGCTGAACGGGCAAGGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGTCGGAAAGGCGGGCCACTTTCTCgg                                                                                                            >  2:499872/1‑251 (MQ=255)
                                                                                                                                    tccggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTGGCCCTTATCTCggggggaggggg                                                                                                  >  1:211457/7‑241 (MQ=255)
                                                                                                                                     ccggcccGGGGATGAAAGGGATGGTTGGCATGGGGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGGGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGa                                                                                                 <  2:178941/245‑1 (MQ=255)
                                                                                                                                         cctggggATGAAAGGGATGATTGGCACGGGGGCTCGCGTGTTTGCAGGGATTCCACGCGCCTGTATTTCCGTTGGT‑CTGATTTCGCAATCATCCTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGGG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGATCTGAAAGAAGTCTTACTGGAGCCGCTGGCGGGGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTa                                                                                            <  2:363392/248‑1 (MQ=255)
                                                                                                                                              ggATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGGGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTGGCCCTTATCCCGGGGG‑GAGGTGATGGGATgggg                                                                                       >  1:114108/1‑248 (MQ=255)
                                                                                                                                                     aGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCa                                                                                                                                                                                         >  1:143058/1‑147 (MQ=255)
                                                                                                                                                     aGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCa                                                                                                                                                                                         <  2:143058/147‑1 (MQ=255)
                                                                                                                                                     aGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGCGCACCTTGc                                                                                <  1:263026/251‑1 (MQ=255)
                                                                                                                                                          tGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGCTGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACCCAATCCTCTTCCGAATACCGTATCAGTTTCTGCGTTCCGCAAAGCGCCTGTG‑GGCGAGCTGAACGGGCAATGCTGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGAAGGGAAGGCCGGCCATTATCTCGGGGG‑GAGGGGAGGGGTGGCACACCTTGCGTggg                                                                           >  2:155773/1‑251 (MQ=255)
                                                                                                                                                               ggCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTGCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCAGTATCTCGGTGG‑TAGGTGATGGTATGCGCACCTTGCGTGGGAtctc                                                                      <  1:23288/251‑1 (MQ=255)
                                                                                                                                                                             ggcttttttaGCGATGTCACGCGTCCGTATTTCCGG‑GGGGCTGATTACGCAATCCTCGTTCTACTTCAGTATGAGTTTCTGCGTTCCGCAAAGCGACTGGG‑TGCGAGCTGAACGGGCAATGCGGTAAGAGTGCTACCTGTAACTGAAAGAAGGCTTATTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGCGCACCTTGCGTGGGATCCCGGCGAAACTCTTTg                                                        <  2:107071/242‑1 (MQ=255)
                                                                                                                                                                                    gCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACt                                                 <  1:183381/251‑1 (MQ=255)
                                                                                                                                                                                        cGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCCTTATCTCGGGGG‑TAGGTGATGGGATGCGCACCTTGCGGGGGATCTCGGCGGAATTCTTTTCCGcgcctgcc                                             >  2:371728/1‑244 (MQ=255)
                                                                                                                                                                                         gATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGAAGGAACGGCTGGCCCTTATCTCGGTGG‑TAGGTGCTGGTATTCGCACCATGCTTGGGATTTCGGCGGAATTTTTTTCCGcggctgccc                                            >  2:265921/1‑243 (MQ=255)
                                                                                                                                                                                         gATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCGGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGGTCTACCTGGAACTGAAAGAAGGATTGCTGGAGCCGCTGGAGGTGTCGGGACGGCTGGCCATTATATCGGTGGGTGGGGGAGGGGTTGGACACCT‑GCGTGGGATCTCGGGGAAATTATTTGCCGcggatgccc                                            >  2:242115/1‑243 (MQ=255)
                                                                                                                                                                                              cACGCGCCCGTATTTCCGT‑GGGGCTGGTTACGCATTCTTCTTCCGATGTCAGGATCAGTTTCTGCGGTCCGCAAATCGACTGTG‑TGCGAGCTGAACGGGGAATGCAGGAAGAGTTCTACCGGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGGCATTATCGCGGTGG‑TAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCTCGCTGGCCCGCGcc                                       <  2:319268/251‑1 (MQ=255)
                                                                                                                                                                                                   gcCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGGGATGGGATGGGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAatat                                  >  1:283040/1‑251 (MQ=255)
                                                                                                                                                                                                    cccGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGGGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCtt                                                          >  1:39199/1‑226 (MQ=255)
                                                                                                                                                                                                    cccGTATTCCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTCCCTGGAACTGAAAGAAGGCGTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCtt                                                          <  2:39199/226‑1 (MQ=255)
                                                                                                                                                                                                      cGTAGTTCCGA‑GGTGCTGATTTCTCAATCATCTCCCGAATACAGTAACAGTTTCTGCGGTCCGCAAAGCTACCGTG‑TGCGAGCTGAGCGGGCAAGGCAGGGAGAGGTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGCGCACCTTGCGTGGGATCGCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCaa                               <  2:5130/251‑1 (MQ=255)
                                                                                                                                                                                                                     gCTGATTACGCATTCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGCGCACCTTGCGTGTGATCTCGGCGAAATTTTTTGCCGCGCTGGGCCGCGTCAATATCAACATTGTCGCCAtt                  <  1:189230/250‑1 (MQ=255)
                                                                                                                                                                                                                     gCTGATTACGCAATCAACTTCCGAAGGCAGTATCAGTTTCTGCGTTCCGCAAGGCGAC‑GTGGTGCGAGCGGAACGGGCAATGCAGGAAGAGTTCTCCCTGGATCTGAAAGAAGGCTCACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCGTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTg                 <  2:50197/251‑1 (MQ=255)
                                                                                                                                                                                                                        gATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAACGCGACTGTG‑TGCGAGCTGAACGGGCACTGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAAGGAACGGCTGGCCATTATCTCGGGGG‑GAGGGGATGGTTAACGCACCATGCGGGGGAACACGGCGAAATTCTTTGCCGGGCTGGCCCGCGCCAATATCCACAATGTACCCATTGCTc              >  2:224212/1‑251 (MQ=255)
                                                                                                                                                                                                                         aTTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCACCCTTGTCGCCATTtctcc             >  1:216865/1‑250 (MQ=255)
                                                                                                                                                                                                                                    atcTTCCGCATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGGAGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGGATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTTCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGAtcttatg  >  1:232817/1‑248 (MQ=255)
                                                                                                                                                                                                                                               aCAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGTAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGCGCACCTTg                                                                                 <  2:490530/161‑1 (MQ=255)
                                                                                                                                                                                                                                               aCAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGCGCACCTTg                                                                                 >  1:490530/1‑161 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             
CAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTG  >  NC_000913/1165‑1642

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: ‑