Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,219,424 G→A 6.4% G277S (GGT→AGT)  aceK → isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,219,4240GA 6.4% 140.2 / 2.4 65G277S (GGT→AGT) aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base G (24/34);  new base A (2/2);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

AGACGGACGACGGCGAGTTATTTATTGATACCTGCCTGACGACGACCGCCGAAGCGAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGT  >  NC_000913/4219278‑4219569
                                                                                                                                                  |                                                                                                                                                 
aGACGGACGACGGCGAGTTAATTATAGCTCCCAGCCAGACGCCGCCCTCCGAAGCGAGAAGTGATGTTGGCATTGCGCGTTCTGATATTATGGTTTATGCGCCGCTGCCCCCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGt                                                                                                                                                 <  2:92378/149‑1 (MQ=255)
    ggacgACGGCGAGTTATTTATTGATACCTGCCTGACGACGACCGCCGAAGCGAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCTCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTaaaa                                                                                                                                             >  2:839048/1‑149 (MQ=255)
    ggacgACGGCGAGTTATTTATTGATACCTGCCTGACGACGACCGCCGAAGCGAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTaaaa                                                                                                                                             >  1:171392/1‑149 (MQ=255)
      acgacgGCGAGTTATTTATTGATACCTGCCTGACGACGACCGCCGAAGCGAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGATGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAacc                                                                                                                                           >  2:1020369/1‑149 (MQ=255)
      acgacgGCGAGTTATATATTGATCCCTGCCTGACGACGACCGCCGAAGCGAGCATTGTCTATGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAacc                                                                                                                                           <  1:165248/149‑1 (MQ=255)
           ggCGAGTTATTTATTGATACCTGCCTGACGACGACCGCCGAAGCGAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGc                                                                                                                                      >  2:510499/1‑149 (MQ=255)
            gCGAGTTATTTATTGATACCTGCCTGACGACGACCGCCGAAGCGAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCt                                                                                                                                     <  1:482686/149‑1 (MQ=255)
            gCGAGTTATTTATTGATACCTGCCTGACGACGACCGCCGAAGCGAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAAGTAAAACCACCGCt                                                                                                                                     >  2:655075/1‑149 (MQ=255)
            gCGAGTTATTTATTGATACCTGCCTGACGACGACCGCCGAAGCGAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAAGTAAAACCACCGCt                                                                                                                                     >  2:187810/1‑149 (MQ=255)
                 ttatttatTGATACCTGCCTGACGACGACCGCCGAAGCGAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAAtt                                                                                                                                <  1:368378/149‑1 (MQ=255)
                   atttattGATACCTGCCTGACGACGACCGCCGAAGCGAGCATTGTTTTTGGCTTTGCGCGTTGTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTCCCAGGTAAAACCACCGCGGAATTGt                                                                                                                              >  2:114890/1‑149 (MQ=255)
                    tttattGATACCTGCCTGACGACGACCGCCGAAGCGAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAaccac                                                                                                                                         <  1:482384/137‑1 (MQ=255)
                    tttattGATACCTGCCTGACGACGACCGCCGAAGCGAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAaccac                                                                                                                                         >  2:482384/1‑137 (MQ=255)
                     ttattGATACCTGCCTGACGACGACCGCCGAAGCGAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGtat                                                                                                                            >  1:398685/1‑149 (MQ=255)
                       attGATACCTGCCTGACGACGACCGCCGAAGCGAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTatat                                                                                                                          >  2:1005278/1‑149 (MQ=255)
                                 gcctgACGACGACCGCCGAACCGAGCATTGTGTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCAAGTGGATACGGGAAATTCTCCCAGGTAAAACCCCCGCTGCATTGTATAAGGCTATCGGc                                                                                                                >  2:585719/1‑149 (MQ=255)
                                     tacgacgacCGCAGAAGCGAGCATTGTTTTTTGATTTGCGCGTTCTTATTTTATGGTTTATGCTCCGCTGCCCGCAGAGCTGGTCGAGTTGCTACGGTAAATTCTGCCATGTAAAACCACCGCTGAATTGTATATGGATATCggattcc                                                                                                            >  2:883726/2‑144 (MQ=255)
                                       cgccgacCGCCGAAGCGAGCTTTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAg                                                                                                          <  1:423366/146‑1 (MQ=255)
                                             ccgccgATGAGAGTATTGTTTTTGGCTTTCCGCGTTCTTATTTAATGGTTTATGCGCCGCTGCCCGCATCGCTGGTCGAGTGGCTACGGGAAATTCTGGCAGGTAAATCCACCGCTGAATTGTATATGGCTATCGGCTGCCTGAAGCAc                                                                                                    <  2:338843/149‑1 (MQ=255)
                                               gccgAAGCGAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGc                                                                                                  >  2:110209/1‑149 (MQ=255)
                                                ccgAAGCGTGCAGTGTTTTTGGCGTTGCGCGTTCTTATATTATGGTTTATGCGCCGCTGCCCGCAGCGCGTGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGcc                                                                                                 <  2:476559/149‑1 (MQ=255)
                                                ccgAAGCGAGCAGTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGcc                                                                                                 <  2:593349/149‑1 (MQ=255)
                                                   aaGCGAGCAGTG‑TTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGGGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCaaaa                                                                                             <  2:130520/149‑1 (MQ=255)
                                                    aGCGCGCATGGTGTTTGGCTGTGCGCGTGCTTATTTTCTGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGGGGCTCCGGGACATTCTGCCAGGTAAAACCACCGCTGACTTGTATCTGGCTATCCGCTGCCAGAAGCACGCCaaaa                                                                                             <  2:853260/149‑1 (MQ=255)
                                                     gCGAGCAGTGTTTTTGGCTTTTCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAAc                                                                                            <  1:829725/149‑1 (MQ=255)
                                                      cGAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGTTCTAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCAGCCAGAAGCACGCCaactcc                                                                                           >  1:901613/1‑145 (MQ=255)
                                                       gAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCg                                                                                          >  1:148922/1‑149 (MQ=255)
                                                         gCATTGTTTGTGGCTATGCGCGTTCATATTTTATGGTTTATGAGGCGCTGGACGCAGCGCT‑GTAGAGTGGCTACGTGAAATTCTGCCAGGTAAAACCACCTCTGAATTTTATATGGCTATAGGCTGCCAGAAGCACTCCAAAACCGaaa                                                                                       <  2:539280/149‑1 (MQ=255)
                                                          cATTGTTTTTTGCTTTGCGCGCTCTTATTATATGGTTAATGCGCCGCTGCCCTCAGCGCAGGTCGAGTCGCTACTGGAAAGTCTGCCATGTAACACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGaaa                                                                                       <  1:488383/149‑1 (MQ=255)
                                                           aTTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAg                                                                                      >  1:425406/1‑149 (MQ=255)
                                                            ttGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGATTTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGc                                                                                     >  1:92246/1‑149 (MQ=255)
                                                             tGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCt                                                                                    >  2:981558/1‑149 (MQ=255)
                                                              gTTTTTGTCTTTTCTCTTTCTTATTTTATGGTTTATGCGCCGATGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACGGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTa                                                                                   <  2:484279/149‑1 (MQ=255)
                                                               tttttGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAACCCACCGCTGAATTGTATATGGCTATCGGCTGCCATAAGCACGCCAAAACCGAAAGCTAc                                                                                  >  1:900218/1‑149 (MQ=255)
                                                               attttGGCTTTGCGCGATCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTCCGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGTCTGCCAGCAGCACGCCATCACCGCAAGCTAc                                                                                  >  2:676979/2‑149 (MQ=255)
                                                                     gCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGATTTGTATATGGCTATCGTCTGCCAGAAGCACGCCAACACCGAAAGCTACCGCGaa                                                                            >  1:219395/1‑149 (MQ=255)
                                                                            gcgTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAATTTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCg                                                                     >  2:159586/1‑149 (MQ=255)
                                                                              gTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGtt                                                                   <  1:359436/149‑1 (MQ=255)
                                                                                   tATTTTTTGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCCCGAATATCTCGTATATCt                                                              >  2:44364/1‑149 (MQ=255)
                                                                                    tttttATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTa                                                             <  2:1021119/148‑1 (MQ=255)
                                                                                      tttATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACa                                                           >  2:855634/1‑149 (MQ=255)
                                                                                         aTGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAggg                                                        <  2:485215/149‑1 (MQ=255)
                                                                                         aTGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAggg                                                        <  1:759493/149‑1 (MQ=255)
                                                                                           ggTTTATGCGACGCTGCCCGCAGCGCTGGTCGAGTCGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCAACCGCGAATATCTCGTTTATCTACAGGGCt                                                      >  2:253512/1‑149 (MQ=255)
                                                                                             tttATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGt                                                    <  1:981558/149‑1 (MQ=255)
                                                                                                aTGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAAt                                                 <  2:405276/149‑1 (MQ=255)
                                                                                                 ttcggcgcTGGCCGGATCGGTTGTCGAGTGGCTTCGGGAAATTCTTGCAGGTTAAACCACCGGTGAATTTTATATGGCTTTCGGCTTCCAGAAGCACGGCAAAACCGAAAGCTACCGCGTATATGTCGTTTATCTACAGGGCTGTAATg                                                <  2:545785/144‑1 (MQ=255)
                                                                                                   cgccgcTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAATACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGTGCTGTAATGAg                                              <  1:772160/149‑1 (MQ=255)
                                                                                                    gccgcTGCCCGAAGTGCTGGTCGAGTGGCTACGGGTAATTCTGCCTGGGAAAACCACCGCTGAATTGTTTATGTCTAGCGGCTGCCAGAAGCTCGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGc                                             <  2:67056/149‑1 (MQ=255)
                                                                                                     ccgcTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCa                                            <  2:171392/149‑1 (MQ=255)
                                                                                                              tcAGAGCTGGTCGATTGGCTAAGGGAAATTCTTCCATGTAAAACCACCGCTTAATTGTATATGGATATAGGATGACAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACATGGCTGTAATGAGCAGTTCATTGa                                   <  1:883726/148‑1 (MQ=255)
                                                                                                              gCAGCGCTGGTCGAGGGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCGGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGa                                   <  2:365925/149‑1 (MQ=255)
                                                                                                               cAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGaa                                  <  1:153324/149‑1 (MQ=255)
                                                                                                                      ggTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCgg                           <  1:1020369/149‑1 (MQ=255)
                                                                                                                      ggTCGAGTGGCTAAGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCgg                           <  1:839048/149‑1 (MQ=255)
                                                                                                                             gggCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCg                    <  2:607372/148‑1 (MQ=255)
                                                                                                                              ggCTACGGGAAATTCTGCCAAGTAAAACCACCGCTGAATTGTATATGGCTATCTGCTGCCAGAATCACGCCAAAACCGAAAGCTTCCGAGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGt                   <  1:187810/149‑1 (MQ=255)
                                                                                                                              ggCTACGGGAAATTCTGCCAAGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGt                   <  1:655075/149‑1 (MQ=255)
                                                                                                                                   cGGGACATTCTGCCAGGTAAAACCACCGCGGAATTGTATATGGCTATCGGCTGGCAGAAGCACGCCTAAACCGAAAGCGAGCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAAt              <  1:676979/149‑1 (MQ=255)
                                                                                                                                      gAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGt           <  1:510499/149‑1 (MQ=255)
                                                                                                                                       aaaTTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGATTATCTACAGGGCCGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTg          <  2:425406/149‑1 (MQ=255)
                                                                                                                                              gCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTgg   >  2:493566/1‑149 (MQ=255)
                                                                                                                                               ccAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTc                     >  2:724614/1‑130 (MQ=255)
                                                                                                                                               ccAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTc                     <  1:724614/130‑1 (MQ=255)
                                                                                                                                               ccAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGgt  <  1:1005278/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGACGGACGACGGCGAGTTATTTATTGATACCTGCCTGACGACGACCGCCGAAGCGAGCATTGTTTTTGGCTTTGCGCGTTCTTATTTTATGGTTTATGCGCCGCTGCCCGCAGCGCTGGTCGAGTGGCTACGGGAAATTCTGCCAGGTAAAACCACCGCTGAATTGTATATGGCTATCGGCTGCCAGAAGCACGCCAAAACCGAAAGCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGT  >  NC_000913/4219278‑4219569

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