Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,219,752 C→T S386F (TCC→TTC)  aceK → isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,219,7520CT96.9% 199.8 / NA 64S386F (TCC→TTC) aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/32);  minor base G (0/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTC  >  NC_000913/4219506‑4219997
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         
aTCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCGGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTtccc                                                                                                                                                                                                                                                        <  2:487459/250‑1 (MQ=255)
   tACAGGGCTGTAATGAGAAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCCTATTTTCCCGGc                                                                                                                                                                                                                                                     >  2:253283/1‑250 (MQ=255)
    aCAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTAACAACTGGTTAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACGTTGTGCTGGAGAAGCGGCATATTTTCCCGGc                                                                                                                                                                                                                                                     <  2:578470/249‑1 (MQ=255)
     cAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCAtt                                                                                                                                                                                                                                                  >  2:216079/1‑251 (MQ=255)
         gCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATg                                                                                                                                                                                                                                              <  1:546332/251‑1 (MQ=255)
            gTAATGAGCAGTTCATTGAAGCTCCGGGTATTCGTGGAATGTTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCTCGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCATGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGTCATTAATgg                                                                                                                                                                                                                                             <  1:106861/249‑1 (MQ=255)
            gTAATGAGCAGTTCATTGAAGCGCCGGGTCTTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGAATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGa                                                                                                                                                                                                                                            >  1:481174/1‑250 (MQ=255)
            gTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGa                                                                                                                                                                                                                                            <  1:168929/250‑1 (MQ=255)
                 tggCAGTTCATGGAAGCTCCGGGGATTCGTGGAATGGTGATGTTGGTGTGTACGCTGTCTGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCATAAAGATAGGGCTGCCGCTCACGTTCCTGCCTGCTATCAACTGGTGAAAGAGCACGAGCGCGGGGCCCTTATGGCGTACACCCAGGAGTTGGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTAc                                                                                                                                                                                                                                       <  1:563147/248‑1 (MQ=255)
                 gagcctttCATTGAAGCGCCGGGTATTCGGGGAATGGTGATGTTGGTGTATACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCCTCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCGCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCGGGAGAAGCGGCATATTTTCCCGGCATTAATTGAATTACt                                                                                                                                                                                                                                      <  2:157360/245‑1 (MQ=255)
                 gAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACTATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCTGGCATTAACGGAATTAc                                                                                                                                                                                                                                       <  2:301817/250‑1 (MQ=255)
                     cgTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCtt                                                                                                                                                                                                                                  >  2:149697/2‑251 (MQ=255)
                         cATTGAAGCGCAGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGCTCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCCTATTTTCCCGGCATTAATGGAATTACTGCTTCAgg                                                                                                                                                                                                                              >  2:581454/1‑251 (MQ=255)
                          aTTGGAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAgg                                                                                                                                                                                                                              >  2:159554/1‑250 (MQ=255)
                                 cgcCGGGTATTCGGGGCATGGGGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGAGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGTCTGCTATCAACTGGGGAAAGAGCACGATCGCGTGGGCCGCATGGCGGACACCCAGGAGCTTGAAAACTCTGTTCTGGAGCAGCGGCATATTTTCCCGGCATTCATGGAATTTCGGCTTCAGGAAGAAGCg                                                                                                                                                                                                                     >  2:341138/1‑251 (MQ=255)
                                    cGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGaa                                                                                                                                                                                                                  <  1:534913/251‑1 (MQ=255)
                                     gggTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTGCGTGCCTGCTATCAACTGGTTAAAGAGCACGATCGCGTGGGCCGAATGGCGGACATCCATGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGaaa                                                                                                                                                                                                                 <  1:67127/251‑1 (MQ=255)
                                         aTTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAAt                                                                                                                                                                                                             >  1:126684/1‑251 (MQ=255)
                                              tagAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTAAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCCTATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATAACCg                                                                                                                                                                                                        >  1:268162/3‑251 (MQ=255)
                                                     gTGATTTTGGGGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGGCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCTGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCgg                                                                                                                                                                                                 <  2:414641/251‑1 (MQ=255)
                                                     gTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCgg                                                                                                                                                                                                 >  2:534046/1‑251 (MQ=255)
                                                       gATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCg                                                                                                                                                                                               >  1:316981/1‑251 (MQ=255)
                                                       gATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCg                                                                                                                                                                                               >  1:556914/1‑251 (MQ=255)
                                                                gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTGCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAGCAAAtt                                                                                                                                                                                       <  1:113309/250‑1 (MQ=255)
                                                                          ccGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCCGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGAATATTTTCCCGGCATTAATGGAATTACTGCATCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGca                                                                                                                                                                            >  1:94697/1‑249 (MQ=255)
                                                                          ccGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCAGTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCTAACAAATTGTGATTCGCCa                                                                                                                                                                           <  2:360980/251‑1 (MQ=255)
                                                                                 ttGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTta                                                                                                                                                                     <  1:302498/250‑1 (MQ=255)
                                                                                 ttGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTta                                                                                                                                                                     <  1:89443/250‑1 (MQ=255)
                                                                                  tGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTtat                                                                                                                                                                    >  1:168743/1‑250 (MQ=255)
                                                                                              cAAAGTCATCAAAGCCAGGTTCGCGCCGCAGAAAGAGATGACTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAATCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCTGCGGa                                                                                                                                                      <  2:246303/251‑1 (MQ=255)
                                                                                                  gTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGt                                                                                                                                                  <  1:159554/251‑1 (MQ=255)
                                                                                                   tcatcaAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTg                                                                                                                                                 >  2:369720/1‑251 (MQ=255)
                                                                                                    catcaAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTg                                                                                                                                                 <  1:357846/250‑1 (MQ=255)
                                                                                                     atcaAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGcc                                                                                                                                               <  1:276593/251‑1 (MQ=255)
                                                                                                                   cgcgcccaaGAAAGAGATGTCTGCCGCTCACGTTAGTGCCTGCTATCCACTGGAGAAAGCGAACGCTCGCGTGGGCCGACTGGCGGAAACCAAGGAGTTTGAAAACTTTGTGCTGGAGCAGCGGCATATTTTCCCGGCATTCATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGCACAAATTGAGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAatatttgac                                                                                                                                 >  2:532460/9‑246 (MQ=255)
                                                                                                                              aaGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAACGAGCACGCTCGCGTGGGCCGAATGGCGGACACCCAGGAGCTTGAACACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATATTTATATTGAGAGGCGGGTGGTGCCGCTTCATATTTGGCTGGAACAAGTg                                                                                                                      >  2:229442/1‑251 (MQ=255)
                                                                                                                                 agaTGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCTCCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGaa                                                                                                                   <  1:253283/251‑1 (MQ=255)
                                                                                                                                 agaTGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGaa                                                                                                                   <  2:141888/251‑1 (MQ=255)
                                                                                                                                       cggCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTca                                                                                                              <  1:23749/248‑1 (MQ=255)
                                                                                                                                         gCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTcagca                                                                                                          <  1:369720/251‑1 (MQ=255)
                                                                                                                                         gCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTcag                                                                                                             <  1:433175/249‑1 (MQ=255)
                                                                                                                                            gCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATCTTTTCCCGGCATTAATGGAATTCCTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCGTCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGGCAGCAGtt                                                                                                       >  1:76169/1‑251 (MQ=255)
                                                                                                                                               cACGGTCGTGCCCGCTATCAACTGGTGAAAGAGCACGAGCGCGGGGGCCGAATGGCGGACTCCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTGCCCGGCATGAATGGAATTACTGCGTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGAGTCGCCATCTTTATATTGAGCGGCGTATGGTGCCGCTCAATATCGGGCTGGATCATGTTGAAGGTcctcagctgcg                                                                                                     <  2:231443/251‑11 (MQ=255)
                                                                                                                                                 cGTTCGGGCCTGCTTTCAACTGGTGAAAGAGCACGATCGCGTTGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCTGCATATTTTCCCGACATTAATGGAATTACTGCTTCAGGAAGCATCGGAAAAAATCACCGATCGCTGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGAAAGGTCAGCAGTTGcgcg                                                                                                  <  2:266351/251‑1 (MQ=255)
                                                                                                                                                  gTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGcgcg                                                                                                  >  1:360620/1‑250 (MQ=255)
                                                                                                                                                            gCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTg                                                                                         <  2:126684/249‑1 (MQ=255)
                                                                                                                                                               aTCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGGGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTTCGCGACGCCATTgaaga                                                                                     >  1:234655/1‑250 (MQ=255)
                                                                                                                                                                      ggTGAAAGGGCACGATTGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGTCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAAGTGTGATTCGCCAGCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACATGTTGAAGGTCAGCAGTTGCGTGACGCCGTTGAAGAATACGGt                                                                              <  2:248203/250‑1 (MQ=255)
                                                                                                                                                                         gAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACg                                                                          <  2:360620/251‑1 (MQ=255)
                                                                                                                                                                             gagCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGAATTTCTTCCTTAGGAAAGCGCGGGAAAAATTACCGGACTCGGGCAACAAATTTTGGATTGCCATCTTTAAATTTGGCCGCGGGAGGGGGCCGCCCATTATTTGGCTGGAAAAGTGGGAGGGCAGGCAGTTGCCGGACCCCTTGAAGAATAAGGGAAAGCTAtt                                                                     >  2:124156/1‑251 (MQ=255)
                                                                                                                                                                               gCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTc                                                                    >  1:502609/1‑251 (MQ=255)
                                                                                                                                                                               gCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTc                                                                    <  2:5005/251‑1 (MQ=255)
                                                                                                                                                                                         gggggacgacaggCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGGGAAGCGGCATATTTTCCCGGCATTAATGAACTTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGAGCCGCTCAATATCTGGCTGGAACAAGGGGAAGGTCAGCAGTTGCGCGACGCCATTTAAGAATACGGTAACGCCATTCGCCAGCTTGc                                                          >  2:238669/12‑251 (MQ=255)
                                                                                                                                                                                                 aaTGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTa                                                  <  1:547328/251‑1 (MQ=255)
                                                                                                                                                                                                           acCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTccc                                        >  1:219485/1‑251 (MQ=255)
                                                                                                                                                                                                           acCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCACGGCATTAATGGAATTACTGCTTCACGAAGCAGCGGAAACAATCACCGATCTCGGCGAACAAATTGTGATTAGCCATCTTTATCTTGTGCGGCGGATAGTGCCGCCCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGCAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTccc                                        >  1:524711/1‑251 (MQ=255)
                                                                                                                                                                                                                    tttGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACAt                               >  1:16641/1‑251 (MQ=255)
                                                                                                                                                                                                                      tGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATg                              <  1:239768/250‑1 (MQ=255)
                                                                                                                                                                                                                        aaaaCTTTGTGCTGGAGAAGGGGAAGATTTGCCCGGCATAAATGGAAGTACAGCTACAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCGGCTCAATAGCTGGCTGGAACAAGTAGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGACAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTg                           <  2:20381/251‑1 (MQ=255)
                                                                                                                                                                                                                             tttGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATAACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTAAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGGCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTa                       >  1:288469/1‑250 (MQ=255)
                                                                                                                                                                                                                                     ggAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGCCGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAAAAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCg              >  2:369553/1‑251 (MQ=255)
                                                                                                                                                                                                                                       agaAGCGGCATAGTTTCCCGGCATTTATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGAGCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGTAACAAGTGTAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGTCGACATGCTGTTTAAAAACTTCGgt            <  1:233674/251‑1 (MQ=255)
                                                                                                                                                                                                                                               cATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGTAAAAAATCACCGATCTCGGCGAACAACTTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGAGGAAGGTCAGCAGTTGCGCGTCGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACCTTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCg    >  2:102007/1‑251 (MQ=255)
                                                                                                                                                                                                                                                atatTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATTTGGCTGGAACAAGTGGAAGGGCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGACAGCTTGCCGATGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTTGGTGTAAccagt   >  2:514541/1‑248 (MQ=255)
                                                                                                                                                                                                                                                  atTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGGGCCGCGCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTc  <  2:491557/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         
ATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTC  >  NC_000913/4219506‑4219997

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