Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 125,083 C→T 37.6% Y689Y (TAC→TAT aceE → pyruvate dehydrogenase E1 component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913125,0830CT37.6% 39.0 / 43.6 58Y689Y (TAC→TATaceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (14/22);  new base T (14/8);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

TCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGC  >  NC_000913/124948‑125225
                                                                                                                                       |                                                                                                                                              
tCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTa                                                                                                                                         >  1:706355/1‑143 (MQ=255)
      cgacaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACACTATCACTACGCTGa                                                                                                                                   >  2:628937/3‑143 (MQ=255)
       cacaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGTGCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAACAACAAGAGAACGTTTACTACTGTATCACTACGCTCaa                                                                                                                                  >  2:515162/1‑143 (MQ=255)
       cacaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTATATCACTACGCTGaa                                                                                                                                  <  1:115295/143‑1 (MQ=255)
       cacaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTATATCACTACGCTGaa                                                                                                                                  >  1:94220/1‑143 (MQ=255)
       cacaTTCAGTCGCTGACTATCCCGAACTGTATCTATTAAGATCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTATATCACTACGGTgga                                                                                                                                  >  1:484821/1‑141 (MQ=255)
          aTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGa                                                                                                                               >  2:567208/1‑143 (MQ=255)
            tcgttCGCTCACGCTCCCGAACTGTAGCCCTTACGGCCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCCGGAGCGTATGTACGGTGAAAAACACGAGAACGTTTACTACTACATCACTACGCTGAACGaaa                                                                                                                             <  2:114780/139‑1 (MQ=255)
                    tGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCCCTACGCTGAACGAAAACTAGcac                                                                                                                     >  2:932797/1‑143 (MQ=255)
                    tGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACcac                                                                                                                     >  2:87553/1‑143 (MQ=255)
                       cTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTAATACATCACTACGCTGAACGAAAACTACCACATg                                                                                                                  >  1:588427/1‑143 (MQ=255)
                          tCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCg                                                                                                               >  1:610165/1‑143 (MQ=255)
                          tCCCGAACTGTATCTCTTACGACCCGGCTTACGCTGACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCg                                                                                                               <  2:438565/143‑1 (MQ=255)
                            ccGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGCCGGTCTGGGGCGTCTGTACGGTCAAAAACACGAGAACGTTTAGTACTATATCACTGCGCTGAAGAAAAGCAACCTCTTGCCgac                                                                                                             >  2:189195/1‑141 (MQ=255)
                            ccGAACTGTATCTCTTACGACCCGGCTTACGCCACAGGTCGTGCAGTCATCATGCATGACGGTCTGGAGCGTATAGACGGTGAAATACAAGAGACAGTTTACTACTATATCACTACGCTGAACTAAAACTACCACATGCCGGc                                                                                                             >  2:332625/1‑143 (MQ=255)
                               aaCTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACCTCACTACGCTGAACGAAGACTCCCACATGCCGGCAAt                                                                                                          >  1:505877/1‑143 (MQ=255)
                               aaCTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCGGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAAt                                                                                                          <  1:638838/143‑1 (MQ=255)
                               aaCTGTAGCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAAt                                                                                                          <  1:443753/143‑1 (MQ=255)
                                          ttACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTg                                                                                               <  1:191415/143‑1 (MQ=255)
                                          ttACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAACGCCGGAAgggg                                                                                               >  2:172213/1‑141 (MQ=255)
                                           tACGACCCGGCTTACGCTTACGAAGTTGCTGGCATCTTACATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGc                                                                                              <  1:889526/143‑1 (MQ=255)
                                              gatccGGCTTACGCGTACGAAGTTGCTGCCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGACCGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGa                                                                                           <  1:395891/140‑1 (MQ=255)
                                              gACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGCCTGGAGCGTATGTCCGGTGAAAAACAAGAGAACGGTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGa                                                                                           <  1:52766/143‑1 (MQ=255)
                                                cccGTCTTATGCTCCCGAGGTTCCTGTGATCAGGCATGACGGTCGGGAGCGTATGTACGGGGAAAAACCAGAGAACGTTTACTACTATAGCGCTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAgg                                                                                         <  2:1991/143‑1 (MQ=255)
                                                cccGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTATATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAgg                                                                                         <  1:382074/143‑1 (MQ=255)
                                                cccGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCAAATGCCGGCCATGCCGGCAGGTGCTGAgg                                                                                         >  2:584821/1‑143 (MQ=255)
                                                   ggCTTACGCTTACGACGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACGACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAg                                                                                      <  1:546421/143‑1 (MQ=255)
                                                   ggCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTATATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAg                                                                                      <  1:603248/143‑1 (MQ=255)
                                                      ttacgcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTATATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTa                                                                                   >  1:265449/1‑143 (MQ=255)
                                                      ttacgcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTa                                                                                   <  1:550264/143‑1 (MQ=255)
                                                         cgcttacgAAGTTGCTGTCAGCGCGCTTGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATcc                                                                                <  2:610165/143‑1 (MQ=255)
                                                           cttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGt                                                                              >  1:399344/1‑143 (MQ=255)
                                                               cgAAGTCGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAg                                                                          <  1:549860/143‑1 (MQ=255)
                                                                        tGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTATATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTGAAGGTGTCTACa                                                                 >  1:76809/1‑143 (MQ=255)
                                                                        tGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTATATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCGGAGGAAGGTATCCGTAAAGGTATCTACa                                                                 >  1:350851/1‑143 (MQ=255)
                                                                               atgcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGGACGAAAACTACCCCTTGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGa                                                          <  2:173384/143‑1 (MQ=255)
                                                                                    tgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTATATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTa                                                     >  2:178393/1‑143 (MQ=255)
                                                                                     gACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTATATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTAt                                                    >  2:601500/1‑143 (MQ=255)
                                                                                     gACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGGTTACTACATTATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTAt                                                    <  1:213185/143‑1 (MQ=255)
                                                                                      ttggTCTGGAGCGTATGTACTGCGAAAAACAAGAGAACGTTTACTACTATATCACTACGCTGTACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTAtt                                                   <  1:471069/141‑1 (MQ=255)
                                                                                         gTCTGGAGCGTATGTGCGGTGAAAAACAAGAGAACGTTTACTACTACCTCGCTACGCTGAACGAAAACTACCACACGCCGGCAATGCCGGCAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGaa                                                <  2:279698/143‑1 (MQ=255)
                                                                                             ggAGCGTATGTACGGAGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGCAAGGTGCTGACGAAGGTACCCGTAAAGGAATCTACAAACACGAAACAAATGAAgcca                                            >  2:478122/1‑140 (MQ=255)
                                                                                                  gtttGTACGGTGAAAAACAAGAGAACGTTTAGTACTACATCCCTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCaaa                                       <  2:676583/140‑1 (MQ=255)
                                                                                                  ggaTGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCaaa                                       <  1:444141/141‑1 (MQ=255)
                                                                                                     tGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGt                                    <  2:30227/143‑1 (MQ=255)
                                                                                                        aCGGTGAAAACAAAGAGAACGTTTACTACTATATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTaaa                                 <  2:403973/143‑1 (MQ=255)
                                                                                                        aCGGTGAAAAACAAGAGAACGTTTACTACTATATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTaaa                                 <  1:153224/143‑1 (MQ=255)
                                                                                                            tGAAACACACGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATCCCCGCAATGCCGGTAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTc                             <  1:34206/143‑1 (MQ=255)
                                                                                                             gAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTCAAGGTATATCTAATCTCGAGACTATTGAAGGTAGCCGAGGTAAAGTTCa                            >  2:412828/1‑143 (MQ=255)
                                                                                                              aaaaaCAAGAGAACGTTTACTACTATATCACTACGCTGAACGAAAACTACCACATGCCGGCAATACCGGAAGAGGCTGAGGAAGGTATCCGTAAAGGTATCTACAAAACAGAAAAGATTGGAGGTAGCAGAGGTACACTTCAg                           >  1:207613/1‑143 (MQ=255)
                                                                                                              aaaaaCAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGACGGTAGCAAAGGTAAAGTTCAg                           >  2:495700/1‑143 (MQ=255)
                                                                                                                    aaGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTc                     >  1:609381/1‑143 (MQ=255)
                                                                                                                      gagaACGTTTACTACTATATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCgg                   >  1:573580/1‑143 (MQ=255)
                                                                                                                      gagaACGTTTACTACTACATCACTACGCAGCTTCCAGACTGCTACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACGCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCgg                   <  2:192643/143‑1 (MQ=255)
                                                                                                                       agaACGTTTACTCCGACATCACGACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGc                  <  2:741896/143‑1 (MQ=255)
                                                                                                                         aaCGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTc                <  1:38250/143‑1 (MQ=255)
                                                                                                                            ggttACTACTACAGCCCTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCgg             <  2:468236/141‑1 (MQ=255)
                                                                                                                             tttACTACTATATCACTACGCTGAACGAAAACGACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGCCTCCggg            >  2:494223/1‑142 (MQ=255)
                                                                                                                                 ctactatATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTa        <  1:515162/143‑1 (MQ=255)
                                                                                                                                 ctactacATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCACAGGTAAAGTTAAGCTACCCGTCTCCGGTTCTa        >  1:904773/1‑143 (MQ=255)
                                                                                                                                   accacATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATc      <  1:567208/140‑1 (MQ=255)
                                                                                                                                    ctacATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATcc     <  1:565883/143‑1 (MQ=255)
                                                                                                                                     tacATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCt    >  2:162104/1‑143 (MQ=255)
                                                                                                                                       cATCACTACGCTGAACGAAAACTTTGGCGCGCCAACAAAGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGc  <  1:389498/143‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                              
TCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGC  >  NC_000913/124948‑125225

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