Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,218,853 C→T Y86Y (TAC→TAT aceK → isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,218,8530CT98.4% 207.2 / NA 65Y86Y (TAC→TATaceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (32/31);  minor base G (0/1);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAAT  >  NC_000913/4218608‑4219099
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
cTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCg                                                                                                                                                                                                                                                   >  2:512737/1‑251 (MQ=255)
cTGGAATTATTGATTGCTCAAAACATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAATGAGCATTATAcc                                                                                                                                                                                                                                                     >  1:83523/1‑249 (MQ=255)
          tgattgCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTACCGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCg                                                                                                                                                                                                                                         >  2:59536/1‑251 (MQ=255)
          tgattgCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGCGCATTATACCCGGCTGTTGCCg                                                                                                                                                                                                                                         >  1:324459/1‑251 (MQ=255)
            attgCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCgg                                                                                                                                                                                                                                        >  1:492872/1‑250 (MQ=255)
              tgCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCGCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGAt                                                                                                                                                                                                                                      <  2:144205/250‑1 (MQ=255)
               gCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGAtt                                                                                                                                                                                                                                     >  1:103699/1‑250 (MQ=255)
                  cAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTAc                                                                                                                                                                                                                                   <  1:239883/249‑1 (MQ=255)
                  cAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTGGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGTCAAAGTACGGACGCGGCATTTTTACTACGTGGTAAATAGCATTACGccctgtttttcgctgattacc                                                                                                                                                                                                                                  <  2:219533/250‑19 (MQ=255)
                        aTTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCtt                                                                                                                                                                                                                           >  1:413687/1‑251 (MQ=255)
                          tttGCAAGGCTTCGATGCTCAGTATGGTCGCTTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCg                                                                                                                                                                                                                         >  1:260994/1‑251 (MQ=255)
                          tttGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTc                                                                                                                                                                                                                          <  2:412162/250‑1 (MQ=255)
                                aGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCGTCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTTCCCGCGCTTCGAGATTg                                                                                                                                                                                                                   <  2:493847/251‑1 (MQ=255)
                                      cGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGgag                                                                                                                                                                                                              <  2:277949/250‑1 (MQ=255)
                                          gCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCtt                                                                                                                                                                                                         <  1:75822/251‑1 (MQ=255)
                                              aGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTGCGAGATTGCGGAGAGCtttttt                                                                                                                                                                                                     <  1:165761/251‑1 (MQ=255)
                                                         ttCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGACAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCgtgt                                                                                                                                                                                           >  1:437952/1‑250 (MQ=255)
                                                           ccTCGAAGTGACCTCCGGTTCGCAGCAGCGTGTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTTGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTAc                                                                                                                                                                                         <  1:475837/250‑1 (MQ=255)
                                                            cTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGACGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTg                                                                                                                                                                                       >  2:390528/1‑251 (MQ=255)
                                                                aaGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCCGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGCGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGCGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCgg                                                                                                                                                                                   >  2:530561/1‑251 (MQ=255)
                                                                aaGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCgg                                                                                                                                                                                   >  2:446132/1‑251 (MQ=255)
                                                                    gACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTa                                                                                                                                                                                <  2:111835/250‑1 (MQ=255)
                                                                    gACCTCCGGTGCGCAGCAGCGGTTCGAACAGGCCGAGTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTGTACGATCAACAAGGTGGTCTGGTCGTGGAGCAAATGCGCTGCATTACTAACGGCCAAAGGACGTACGCGGCATTTTTACTACGTGTTAAAGAGCATTAGACCCGGCTGTTGCCGGATTACCCGCGCGTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTAt                                                                                                                                                                               <  1:19324/251‑1 (MQ=255)
                                                                                gcagcaGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCgg                                                                                                                                                                                                                >  2:65538/1‑206 (MQ=255)
                                                                                gcagcaGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCgg                                                                                                                                                                                                                <  1:65538/206‑1 (MQ=255)
                                                                                            cGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTccc                                                                                                                                                        <  2:413687/250‑1 (MQ=255)
                                                                                            cGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTccc                                                                                                                                                        <  1:309048/250‑1 (MQ=255)
                                                                                              aaCAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCg                                                                                                                                                       >  1:421377/1‑249 (MQ=255)
                                                                                                      gACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCtttt                                                                                                                                             >  1:315098/1‑251 (MQ=255)
                                                                                                       aCTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTTTTTGACCACCGCTCGCTTACTCCCGAGCGGCttttt                                                                                                                                            >  1:293204/1‑251 (MQ=255)
                                                                                                            cATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCtt                                                                                                                                       <  2:103699/251‑1 (MQ=255)
                                                                                                               gCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAg                                                                                                                                    <  1:512737/251‑1 (MQ=255)
                                                                                                                cTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGc                                                                                                                                   >  1:263175/1‑251 (MQ=255)
                                                                                                                    ccAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGctct                                                                                                                                <  2:263175/250‑1 (MQ=255)
                                                                                                                        cagGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGATGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGCTTGCGGAGAGCTTTTTTAACTCCGGGTACTGTCGGTTATTTGACCACCGCTCGCTTATTCCCGAGCGGTTTTTTTTCTTTAGCTCTCaccc                                                                                                                           >  2:487240/1‑248 (MQ=255)
                                                                                                                            cGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTGCCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGCGCGGCTTTTTATCTTTAGCTCTCAGCCag                                                                                                                         >  2:344989/1‑249 (MQ=255)
                                                                                                                                gAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAgcgc                                                                                                                    <  1:515965/250‑1 (MQ=255)
                                                                                                                                 aaaaaCCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGCGCGGCTTTTTATCTTTATCTCTCAGCCAGAGcgccg                                                                                                                  >  2:344131/1‑251 (MQ=255)
                                                                                                                                 aaaaaCCGTATCCATCTTTACGATCATCACGTTGGACTGGTAGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGAGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTACCTCCGTGTACTGTCGGTTATTTGACCAACGCTCCCCTACTCCCGGGCGGCTTTTTATCTTTaaattccagccacaacgccg                                                                                                                  >  2:547346/1‑232 (MQ=255)
                                                                                                                                       agtatcctacTTTACGATCAATCCGTTGGGCTGGTCGTGGAGCATCTGCGCTGCCTTACTAACGGCCAAAGTACGGACGTGGCATGTTTTCCACGTGGTAAAGAGCATTATCCCCGGCTGTTGCCGGATTACCCGCGCTTCGATATTGCGGAGCGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCGCTCTGACTCCCTAGCGGCTTTTTACCTTTATCTCGCAGCCATAGCGCCGCtctc                                                                                                             <  2:373680/241‑4 (MQ=255)
                                                                                                                                        gTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGATAGCTTTTTTaa                                                                                                                                                                                                   <  2:316614/163‑1 (MQ=255)
                                                                                                                                        gTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTaa                                                                                                                                                                                                   >  1:316614/1‑163 (MQ=255)
                                                                                                                                         tATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTCTCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGt                                                                                                           >  1:291578/1‑250 (MQ=255)
                                                                                                                                            ccATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTAccc                                                                                                       >  2:133218/1‑250 (MQ=255)
                                                                                                                                            ccATCTTTACGATCATCACGTTGGTCTGGGCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATGTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCGTCGAGATTGCGGAGAGCTTATTTATCTCCGTGTACTTTCTGTTAGTGGGCCACCGCGCTCCTTCACCCCAGCGGCTTTTTATCCTTATCTCTCAGTCACATcgccgccttgctacca                                                                                                       <  2:529247/251‑11 (MQ=255)
                                                                                                                                                      gATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCcccgccc                                                                                             <  1:294998/251‑1 (MQ=255)
                                                                                                                                                       atcatcACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTCCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCcccgcccg                                                                                            >  2:338154/1‑251 (MQ=255)
                                                                                                                                                                   ggtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGaa                                                                                    >  1:210386/1‑248 (MQ=255)
                                                                                                                                                                        ggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTc                                                                            <  1:31999/251‑1 (MQ=255)
                                                                                                                                                                           cGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGGCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTAGCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCAcc                                                                         <  2:27973/251‑1 (MQ=255)
                                                                                                                                                                                  cAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCa                                                                  <  1:272657/251‑1 (MQ=255)
                                                                                                                                                                                    aCTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAACGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCCTTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGCTCACg                                                                >  2:332853/1‑251 (MQ=255)
                                                                                                                                                                                                            aaGTACGGACGCGGCATTTTTACTACGTGTTCAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAACCAGAGCGCCGCTTTCGTACCATTCCCCCCCCGCTGGCGAAAGACTTTCACCCCCATCACGCCTGGGCATCTCTACTGATGggcg                                        >  1:187952/1‑251 (MQ=255)
                                                                                                                                                                                                              gTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGtt                                      >  1:231587/1‑251 (MQ=255)
                                                                                                                                                                                                                 cGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTATACCCGGCTGGTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTCTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATc                                   <  2:128020/251‑1 (MQ=255)
                                                                                                                                                                                                                            tttttACTACGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTAcc                        >  1:116562/1‑251 (MQ=255)
                                                                                                                                                                                                                                 actacGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCCCGGCTGGGAATCTCTACTAATGCGGGTTATCAGCGACCTACCGctgc                   >  2:290316/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cGTGTTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACGCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcct               <  2:293204/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gtTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcctgc             <  2:262325/250‑1 (MQ=255)
                                                                                                                                                                                                                                        gtTAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCAGTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcctgc             <  2:125232/250‑1 (MQ=255)
                                                                                                                                                                                                                                          tAAAGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGctgcgcctgcgc           <  1:133218/250‑1 (MQ=255)
                                                                                                                                                                                                                                           aaaGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTg         <  2:315098/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aaGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACgg                                                               <  2:439637/196‑1 (MQ=255)
                                                                                                                                                                                                                                            aaGAGCATTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACgg                                                               >  1:439637/1‑196 (MQ=255)
                                                                                                                                                                                                                                                  aTTATACCCGGCTGTTGCCGGATTACCCGCGCTTCGATATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTCCCGCTGCGCCTGCGCTGGCAGAAt  <  2:73391/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
CTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAAT  >  NC_000913/4218608‑4219099

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: