Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,425 A→G P363P (CCA→CCG 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,4250AG96.7% 169.7 / ‑2.5 60P363P (CCA→CCGthrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (29/29);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01

ATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTG  >  NC_000913/1180‑1662
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
aTTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAg                                                                                                                                                                                                                                            <  2:430629/251‑1 (MQ=255)
      actACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCCATCATCTTCCGAATACAGTATCCGTTTCTGCGTTCCGCaaaacggctg                                                                                                                                                                                                                                      >  2:171764/3‑244 (MQ=255)
       agtccGGAAATCCTCAAGCTCCAGGTACGCTCATTGGGGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACtgt                                                                                                                                                                                                                                     <  2:25648/248‑1 (MQ=255)
          ccGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTTAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCTTTCCGCAAAGCGACtgtgtg                                                                                                                                                                                                                                 <  2:227230/251‑1 (MQ=255)
          ccGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGGGTTTGCAGCGGTGTCACGCGCCCCGATTTTCGTGGTGCTGATTACGCGATCATCTTCCGAAAACAGTATCAGTTTCTGCGTTTCGCAAAGCGACagtgtg                                                                                                                                                                                                                                 >  2:195379/1‑251 (MQ=255)
           cGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGc                                                                                                                                                                                                                                <  1:171764/251‑1 (MQ=255)
                aTCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCt                                                                                                                                                                                                                           <  1:445854/251‑1 (MQ=255)
                aTCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGGGCGAGCt                                                                                                                                                                                                                           <  2:219384/251‑1 (MQ=255)
                           ccgtgtacgcgcagcggttccatccGTGATGAAGACTAATTACCGGGCAAGTGCATGTCCATTCTGAATAAGATGGCAATGTGCCGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTTCAGCGATGTCACTCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCTGAATTCAGTATAAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAAt                                                                                                                                                                                                                <  2:437118/228‑1 (MQ=255)
                           ccAGGTACGCTCATTGGTGCCAGCCGTGATGAATACGAATGACCGGTCAAGGGCATTTCCAATCTGAATAAGAGGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGTCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAAt                                                                                                                                                                                                                <  2:462611/251‑1 (MQ=255)
                           ccAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTCTTTCCGTGGTGCTGATTACGCCATCATCTTCCGAATACAGAATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAAt                                                                                                                                                                                                                >  2:170842/1‑251 (MQ=255)
                             aGGTACGCTCATTGGAGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTATGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGc                                                                                                                                                                                                              >  1:303970/1‑251 (MQ=255)
                               gTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAg                                                                                                                                                                                                            <  1:453856/251‑1 (MQ=255)
                                             gCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCt                                                                                                                                                                                              <  1:244995/251‑1 (MQ=255)
                                              ccAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCCGGAAGAGTTCTACCTg                                                                                                                                                                                             >  2:183381/1‑251 (MQ=255)
                                                    gTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGATCTACCTGGAACTg                                                                                                                                                                                       >  1:441056/1‑251 (MQ=255)
                                                       aggaAGCCGAATTAACGGTCAAGGGCATTTCCAATCTGAATAATATGGGAATGTTTAGCGTTTACGGCACGGGGATGAAAGGGATGGTGGGCATGGGGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTGTCGGCGTTCCGCAAAGCGACGGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCGACCTGGAACTGaaa                                                                                                                                                                                    <  2:110819/249‑1 (MQ=255)
                                                        tgaAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGCTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCCACTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCCACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGaa                                                                                                                                                                                     >  1:66902/1‑249 (MQ=255)
                                                               gAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGGTTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGGGCTGAACGGGAAATGACGGAAGGGTTCTACCTGGAACTGAAAGACGGCtt                                                                                                                                                                            >  2:258914/1‑251 (MQ=255)
                                                                          cAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGGGcc                                                                                                                                                                   >  2:426550/1‑249 (MQ=255)
                                                                          cAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGcc                                                                                                                                                                   <  1:426550/249‑1 (MQ=255)
                                                                           aaGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGGGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGCATGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTCCTGGAGCCGCt                                                                                                                                                                >  1:304401/1‑251 (MQ=255)
                                                                           aaGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGGATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGGGCTGAACGGGCAATGCAGGAAGATGTCTACCTGGGACTGAAAGAAGGCTTACATGAGCCGCt                                                                                                                                                                >  2:481842/1‑251 (MQ=255)
                                                                              ggCATTTCCAATCTGGATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGc                                                                                                                                                             <  2:274454/251‑1 (MQ=255)
                                                                                      cAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGGGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGGTTATGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGGGTTCTACCTGGAACTGACAGAAGACTTACTGGAGCCGCTGGcggttaccg                                                                                                                                                     >  1:260793/1‑243 (MQ=255)
                                                                                               taataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGTTTACGCAATCATCTTCCGCATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGaacgaaccgccgg                                                                                                                                            >  1:209783/5‑239 (MQ=255)
                                                                                                  taTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTggac                                                                                                                                          >  1:400077/2‑248 (MQ=255)
                                                                                                   aTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGTATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTGCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCAt                                                                                                                                        <  1:305385/251‑1 (MQ=255)
                                                                                                        aaTGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTGGCCCTTAtct                                                                                                                                   >  1:178941/1‑251 (MQ=255)
                                                                                                         atggtCAGCGTGTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCAGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTAGCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCGGGCGGTGACGGTTCTCCTGGCCATTAtctc                                                                                                                                  <  1:197310/247‑1 (MQ=255)
                                                                                                         aggTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTAtctc                                                                                                                                  <  1:234003/249‑1 (MQ=255)
                                                                                                           gTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACCCGCCCGTCTTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGAGACTGTGTGCGAGCTGAACGGGCAAGGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGTCGGAAAGGCGGGCCACTTTCTCgg                                                                                                                                >  2:499872/1‑251 (MQ=255)
                                                                                                                     tccggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTGGCCCTTATCTCggggggaggggg                                                                                                                      >  1:211457/7‑241 (MQ=255)
                                                                                                                      ccggcccGGGGATGAAAGGGATGGTTGGCATGGGGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGGGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGa                                                                                                                     <  2:178941/245‑1 (MQ=255)
                                                                                                                          cctggggATGAAAGGGATGATTGGCACGGGGGCTCGCGTGTTTGCAGGGATTCCACGCGCCTGTATTTCCGTTGGTCTGATTTCGCAATCATCCTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGGGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGATCTGAAAGAAGTCTTACTGGAGCCGCTGGCGGGGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTa                                                                                                                <  2:363392/248‑1 (MQ=255)
                                                                                                                               ggATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGGGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTGGCCCTTATCCCGGGGGGAGGTGATGGGATgggg                                                                                                           >  1:114108/1‑248 (MQ=255)
                                                                                                                                      aGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCa                                                                                                                                                                                                             <  2:143058/147‑1 (MQ=255)
                                                                                                                                      aGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCa                                                                                                                                                                                                             >  1:143058/1‑147 (MQ=255)
                                                                                                                                      aGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGc                                                                                                    <  1:263026/251‑1 (MQ=255)
                                                                                                                                           tGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGCTGTCACGCGCCCGTATTTCCGTGGTGCTGATTACCCAATCCTCTTCCGAATACCGTATCAGTTTCTGCGTTCCGCAAAGCGCCTGTGGGCGAGCTGAACGGGCAATGCTGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGAAGGGAAGGCCGGCCATTATCTCGGGGGGAGGGGAGGGGTGGCACACCTTGCGTggg                                                                                               >  2:155773/1‑251 (MQ=255)
                                                                                                                                                ggCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTGCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCAGTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGAtctc                                                                                          <  1:23288/251‑1 (MQ=255)
                                                                                                                                                              ggcttttttaGCGATGTCACGCGTCCGTATTTCCGGGGGGCTGATTACGCAATCCTCGTTCTACTTCAGTATGAGTTTCTGCGTTCCGCAAAGCGACTGGGTGCGAGCTGAACGGGCAATGCGGTAAGAGTGCTACCTGTAACTGAAAGAAGGCTTATTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCCCGGCGAAACTCTTTg                                                                            <  2:107071/242‑1 (MQ=255)
                                                                                                                                                                     gCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACt                                                                     <  1:183381/251‑1 (MQ=255)
                                                                                                                                                                         cGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCCTTATCTCGGGGGTAGGTGATGGGATGCGCACCTTGCGGGGGATCTCGGCGGAATTCTTTTCCGcgcctgcc                                                                 >  2:371728/1‑244 (MQ=255)
                                                                                                                                                                          gATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGAAGGAACGGCTGGCCCTTATCTCGGTGGTAGGTGCTGGTATTCGCACCATGCTTGGGATTTCGGCGGAATTTTTTTCCGcggctgccc                                                                >  2:265921/1‑243 (MQ=255)
                                                                                                                                                                          gATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCGGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGGTCTACCTGGAACTGAAAGAAGGATTGCTGGAGCCGCTGGAGGTGTCGGGACGGCTGGCCATTATATCGGTGGGTGGGGGAGGGGTTGGACACCTGCGTGGGATCTCGGGGAAATTATTTGCCGcggatgccc                                                                >  2:242115/1‑243 (MQ=255)
                                                                                                                                                                               cACGCGCCCGTATTTCCGTGGGGCTGGTTACGCATTCTTCTTCCGATGTCAGGATCAGTTTCTGCGGTCCGCAAATCGACTGTGTGCGAGCTGAACGGGGAATGCAGGAAGAGTTCTACCGGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGGCATTATCGCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCTCGCTGGCCCGCGcc                                                           <  2:319268/251‑1 (MQ=255)
                                                                                                                                                                                    gcCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGGGATGGGATGGGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAatat                                                      >  1:283040/1‑251 (MQ=255)
                                                                                                                                                                                     cccGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGGGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCtt                                                                              >  1:39199/1‑226 (MQ=255)
                                                                                                                                                                                     cccGTATTCCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTCCCTGGAACTGAAAGAAGGCGTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCtt                                                                              <  2:39199/226‑1 (MQ=255)
                                                                                                                                                                                       cGTAGTTCCGAGGTGCTGATTTCTCAATCATCTCCCGAATACAGTAACAGTTTCTGCGGTCCGCAAAGCTACCGTGTGCGAGCTGAGCGGGCAAGGCAGGGAGAGGTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCGCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCaa                                                   <  2:5130/251‑1 (MQ=255)
                                                                                                                                                                                                      gCTGATTACGCATTCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGTGATCTCGGCGAAATTTTTTGCCGCGCTGGGCCGCGTCAATATCAACATTGTCGCCAtt                                      <  1:189230/250‑1 (MQ=255)
                                                                                                                                                                                                      gCTGATTACGCAATCAACTTCCGAAGGCAGTATCAGTTTCTGCGTTCCGCAAGGCGACGTGGTGCGAGCGGAACGGGCAATGCAGGAAGAGTTCTCCCTGGATCTGAAAGAAGGCTCACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCGTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTg                                     <  2:50197/251‑1 (MQ=255)
                                                                                                                                                                                                         gATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAACGCGACTGTGTGCGAGCTGAACGGGCACTGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAAGGAACGGCTGGCCATTATCTCGGGGGGAGGGGATGGTTAACGCACCATGCGGGGGAACACGGCGAAATTCTTTGCCGGGCTGGCCCGCGCCAATATCCACAATGTACCCATTGCTc                                  >  2:224212/1‑251 (MQ=255)
                                                                                                                                                                                                          aTTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCACCCTTGTCGCCATTtctcc                                 >  1:216865/1‑250 (MQ=255)
                                                                                                                                                                                                                     atcTTCCGCATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGGAGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGGATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTTCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGAtcttatg                      >  1:232817/1‑248 (MQ=255)
                                                                                                                                                                                                                                aCAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGTAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTg                                                                                                     <  2:490530/161‑1 (MQ=255)
                                                                                                                                                                                                                                aCAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTg                                                                                                     >  1:490530/1‑161 (MQ=255)
                                                                                                                                                                                                                                       caggttCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCg    <  1:129368/247‑1 (MQ=255)
                                                                                                                                                                                                                                         gTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCCGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCAATATCTCGGTGGTGGGGGGTGGTATGCGAAACTTGCGTGGGTTCTCGGCGAAATTCTTTGCACCGCATGCCCGCGCCACAATCAACATAGTTGCCCTTTCTACGGGGTCTTTTGCACGCCCAATCTCTGTtgtg  >  2:145162/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
ATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTG  >  NC_000913/1180‑1662

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