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,7520CT98.4% 209.9 / NA 63S386F (TCC→TTC) aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (29/33);  minor base G (0/1);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCA  >  NC_000913/4219520‑4219998
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
aattaGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAAtt                                                                                                                                                                                                                                      <  1:408847/247‑1 (MQ=255)
       agtttcatGAAGCGCCGGGTATTCTTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATATGGTATTCAAAGTCATCAAAGACAGGTGCGCGACGCAGAAGGAGATGTCTGCCGATCACGTTCGTGCCTGCTATCAACTGTTGAAAGAGCACGAGCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCtt                                                                                                                                                                                                                               <  2:576604/244‑1 (MQ=255)
              tGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGaa                                                                                                                                                                                                                         <  2:228840/250‑1 (MQ=255)
                aaGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGAGCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAACCTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAagca                                                                                                                                                                                                                      <  1:687622/251‑1 (MQ=255)
                   cgcCGGGTATTCGTGGATTGGTGATGGTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCGCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGTTCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGTAATTACTGCTGCAGGAAGCAGCg                                                                                                                                                                                                                   <  2:437965/251‑1 (MQ=255)
                          tATTCGTGGAATGGTGATGTTGGTGTTTACGCCGCCGGGCTTTGATCGGGTATTCAAAGTCATAAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGGGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTAATGCTTCAGGAAGCAGCGGaaaaa                                                                                                                                                                                                             >  2:182973/1‑250 (MQ=255)
                                    aTGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGAAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGAtct                                                                                                                                                                                                  >  1:273351/1‑251 (MQ=255)
                                      ggTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCg                                                                                                                                                                                                >  1:131307/1‑251 (MQ=255)
                                      ggTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCg                                                                                                                                                                                                >  2:498988/1‑251 (MQ=255)
                                        tGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGc                                                                                                                                                                                              >  1:552217/1‑251 (MQ=255)
                                        ggATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCATCGGAAAAAATCACCGATCTCGGc                                                                                                                                                                                              <  2:493223/250‑1 (MQ=255)
                                         gATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGc                                                                                                                                                                                              <  1:498988/250‑1 (MQ=255)
                                               ggTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACaa                                                                                                                                                                                        >  2:639572/1‑250 (MQ=255)
                                                        gCTGCCGGGCTTTGCTCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTc                                                                                                                                                                              >  2:679971/1‑251 (MQ=255)
                                                             cGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGCTCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCAt                                                                                                                                                                         >  2:723565/1‑251 (MQ=255)
                                                                  tttGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCttt                                                                                                                                                                     >  1:204230/1‑250 (MQ=255)
                                                                   ttGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCttt                                                                                                                                                                     <  2:177530/249‑1 (MQ=255)
                                                                    tGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTtat                                                                                                                                                                   >  1:574256/1‑250 (MQ=255)
                                                                    tGATCGGGTATTCAAAGTCATAAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTtat                                                                                                                                                                   >  2:329650/1‑250 (MQ=255)
                                                                      aTCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATAtt                                                                                                                                                                >  1:390750/1‑251 (MQ=255)
                                                                                           aaGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCCGCTATCAACTGGTGAAAGCGCACGATCGCGTGGGCCGAATGGCAGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACAGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGCTGGAGACGCTc                                                                                                                                           >  1:318728/1‑251 (MQ=255)
                                                                                                       cgccgcAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATctggctg                                                                                                                               >  1:585680/1‑251 (MQ=255)
                                                                                                         ccgcAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGa                                                                                                                             <  1:558117/251‑1 (MQ=255)
                                                                                                            cAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGCAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAATTCAACGATCTCGGCGAACAAATTGTGATTCGACATCTTTATATTGAGCGGCGGATGGTACCGCTCACTATctggctggcac                                                                                                                           >  2:287252/1‑247 (MQ=255)
                                                                                                            cAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACa                                                                                                                          <  1:201137/251‑1 (MQ=255)
                                                                                                                  gagaTGTCTGCCGATCCCGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGAAACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGCGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTgg                                                                                                                     >  2:700226/1‑250 (MQ=255)
                                                                                                                     aTGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAgg                                                                                                                 >  2:582067/1‑251 (MQ=255)
                                                                                                                       ggcTGCCGCCCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTc                                                                                                               <  2:194875/249‑1 (MQ=255)
                                                                                                                                 cACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTgc                                                                                                      <  1:8115/250‑1 (MQ=255)
                                                                                                                                        gTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGcc                                                                                              <  1:23899/251‑1 (MQ=255)
                                                                                                                                           ccTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCAtt                                                                                           <  1:679971/251‑1 (MQ=255)
                                                                                                                                            cTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAACTCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTg                                                                                          >  2:566339/1‑251 (MQ=255)
                                                                                                                                                 agcAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTGCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGTCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATTGTGCCGCTCAATATCAGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTgaagaa                                                                                     <  2:318250/249‑1 (MQ=255)
                                                                                                                                                    aaCTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATAc                                                                                  >  2:244777/1‑251 (MQ=255)
                                                                                                                                                        ggTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTc                                                                                                                                                                              >  2:579723/1‑155 (MQ=255)
                                                                                                                                                        ggTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTc                                                                                                                                                                              <  1:579723/155‑1 (MQ=255)
                                                                                                                                                          tGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAAc                                                                            >  2:69320/1‑251 (MQ=255)
                                                                                                                                                               gagcacgctcgagtgaggagAATGGCCGAAACCCCGGAGTTTGAAAACTTTGTGCAGGCGAAGCGGCATCTTTTCCCGGCATTAATGGACTTACTGCTTAAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACACATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCAGCTCAATATCTGGCTGGAACACGTGGACGGTAAGCAGTTGCGCGGCGCCATTGAAGAATACGGTAACGCTAt                                                                       >  2:115722/20‑251 (MQ=255)
                                                                                                                                                                agCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTAtt                                                                      <  1:723565/251‑1 (MQ=255)
                                                                                                                                                                 gCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTc                                                                     <  1:69320/251‑1 (MQ=255)
                                                                                                                                                                  cACGCTCGCGTGGGCCGAATGGCAGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGCAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCCGACAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGTTCAGAAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCg                                                                    <  2:168431/251‑1 (MQ=255)
                                                                                                                                                                  cACGATCGCGTGGGCCGAATGGCGGACACCCATGATTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATCGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGTTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCg                                                                    <  1:624974/251‑1 (MQ=255)
                                                                                                                                                                        cgcgTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAAGCACCGATCGCGGCGAACAAATTGTTATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGGGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCt                                                              <  2:732952/251‑1 (MQ=255)
                                                                                                                                                                           gTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGc                                                           <  2:363993/251‑1 (MQ=255)
                                                                                                                                                                             gggCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTAGTCGCCAGCTTGcc                                                          <  2:131307/250‑1 (MQ=255)
                                                                                                                                                                                 cGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCgctgc                                                     >  2:567526/1‑251 (MQ=255)
                                                                                                                                                                                  gAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCAGTAATTGAATTACTGCTTCAGGATGCAGCGGAAAAAGTCACCGATCTCGGCGAACAATTTGTGATTCTCGATCTTTGTATTGAGCGGCGGATGTTTCCGCTCAATATCTGGCTGGAACAGGGGGAAGGTCAGCAGTTGCGCAACACCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCgctgct                                                    <  1:430282/251‑1 (MQ=255)
                                                                                                                                                                                        cGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACACGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACAtt                                              >  1:665514/1‑251 (MQ=255)
                                                                                                                                                                                        cGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACAtt                                              <  1:329650/251‑1 (MQ=255)
                                                                                                                                                                                                    aGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGAc                                  >  2:315177/1‑251 (MQ=255)
                                                                                                                                                                                                        tGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATg                               <  1:536098/250‑1 (MQ=255)
                                                                                                                                                                                                             aCTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGttt                         <  1:182973/251‑1 (MQ=255)
                                                                                                                                                                                                             aCTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGttt                         >  2:271526/1‑251 (MQ=255)
                                                                                                                                                                                                                gtgtgCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTa                        <  1:642574/248‑1 (MQ=255)
                                                                                                                                                                                                                 tgtgCTGGAGAAGCGGCATACTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATGCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGTTCAGCAGTTGCGCGATGCCATTGAAGAATACGGTAACGCGATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCGGTTTaaa                      <  2:286521/250‑1 (MQ=255)
                                                                                                                                                                                                                    gCTGGAGAAGCGGCATATTTTCCCGGCATTAATGTAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGCCGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAAc                   <  1:274325/250‑1 (MQ=255)
                                                                                                                                                                                                                        gagaAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCgg              <  2:60180/251‑1 (MQ=255)
                                                                                                                                                                                                                          gaAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGAGAACAACTTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGTCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGACTACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGCCATGCTGTTTAAAAACTTCGgtg            >  2:38780/1‑251 (MQ=255)
                                                                                                                                                                                                                            aGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTc          <  2:585680/251‑1 (MQ=255)
                                                                                                                                                                                                                              cGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCAACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCAc        <  2:215898/251‑1 (MQ=255)
                                                                                                                                                                                                                              cGGCAGATTTTCCCGGCACTAATGGAATTACTGCTTCATGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGGTGGTGGCGCGCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCGTTGAAGAATACTGTAACGGTATTCGCCAGCTTTCCTCTGCTAACATTTTCCCTGTCGACATGCTGTTTAAAAACTTCGGTGTCAc        <  1:698942/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCg     >  2:39213/1‑251 (MQ=255)
                                                                                                                                                                                                                                    atTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCCGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTCAAAACTTCGGTGTCACCCGTcc  >  1:92263/1‑250 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
AATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCA  >  NC_000913/4219520‑4219998

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