Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_tpiA 2,795,396 C→A Q382K (CAG→AAG)  gabP → gamma‑aminobutyrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_tpiA2,795,3960CA98.5% 218.7 / ‑5.9 65Q382K (CAG→AAG) gabPgamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (40/24);  total (41/24)
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.88e-01

CCCCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGT  >  NC_000913_3_vch_tpiA/2795250‑2795537
                                                                                                                                                  |                                                                                                                                             
ccccGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAg                                                                                                                                             >  1:182321/1‑149 (MQ=255)
      aCGTGGCGGTGTTACTCTCCACCGGCGCGGCATTATTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCgg                                                                                                                                       >  1:777196/1‑149 (MQ=255)
       cGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGa                                                                                                                                      <  2:289364/149‑1 (MQ=255)
                     tctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTCCGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTg                                                                                                                        >  2:325778/1‑149 (MQ=255)
                     tctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTg                                                                                                                        >  2:407865/1‑149 (MQ=255)
                     tctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTg                                                                                                                        >  1:175791/1‑149 (MQ=255)
                     tctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTg                                                                                                                        <  1:725619/149‑1 (MQ=255)
                       tcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCg                                                                                                                      >  2:713316/1‑149 (MQ=255)
                        ccACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACCCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGa                                                                                                                     <  1:46993/149‑1 (MQ=255)
                                  ggCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaaggca                                                                                                           >  2:24294/1‑147 (MQ=255)
                                   gCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTTCGAACAgaaagaag                                                                                                          >  2:49159/1‑149 (MQ=255)
                                      ttCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGa                                                                                                       >  1:348365/1‑149 (MQ=255)
                                         ttAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAAt                                                                                                    >  2:312958/1‑149 (MQ=255)
                                             ctgttgtggtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTACACCTCGTCGGTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGTAGAAGGAAGCGAAATTCGc                                                                                                <  2:588572/144‑1 (MQ=255)
                                              ggtggtggtGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCt                                                                                               <  2:857529/149‑1 (MQ=255)
                                                      tgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAACGTTGCGGATGCGCAAAATTCAGCGAGCAGAAGGAAGCGAAATTCGCTTGCGaagg                                                                                       >  1:929303/1‑147 (MQ=255)
                                                      tgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtg                                                                                       >  2:496995/1‑149 (MQ=255)
                                                      tgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtg                                                                                       <  1:162939/149‑1 (MQ=255)
                                                       gAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTGGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtgt                                                                                      <  1:675328/149‑1 (MQ=255)
                                                           tATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCGCTTTCAAATTTGCGGATGCGCAAAATTCTGTGAGCAACAGGAAGCGAAATTCGCTTGCTACTGTGGCt                                                                                  >  2:112373/1‑149 (MQ=255)
                                                           tATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCt                                                                                  <  1:975379/149‑1 (MQ=255)
                                                                 gcgcCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATcc                                                                            <  1:231191/149‑1 (MQ=255)
                                                                  cgcCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCg                                                                           <  2:934948/149‑1 (MQ=255)
                                                                   gcCTGCGAAGGTATTTAAATTTCTGATCGACAGATCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTATGCGAGCAGGAGGAAGCCAAATTCGCTTGGGGATGTGGCTTTATCCGt                                                                          >  2:75642/1‑149 (MQ=255)
                                                                      tGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCGGTTTCAAGGTTGCGGATGCGCAAATTTTTGCGGGAAGAAGGAAGCGTATTTCGCTTGCGAATGTGGCTTTATCAGTGGc                                                                       >  1:368337/1‑149 (MQ=255)
                                                                         gAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTAACGAGGCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCa                                                                    >  1:984850/1‑149 (MQ=255)
                                                                             gTATTTAAGTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGACGCGCAAAATTCTGCGAGCAGAAGGAAGCGCAATTCGCTTGCGAATGTGGCTTTCTCCGTGGCTCACCTg                                                                >  1:638869/1‑149 (MQ=255)
                                                                              tATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTTGCTCACCTgg                                                               <  1:1043115/149‑1 (MQ=255)
                                                                                tttAAATTTCTGGTCGACAGGCACGGCGCTATCGCCCTGCTGGTTTATTTAGGCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGCAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCt                                                             <  2:638869/149‑1 (MQ=255)
                                                                                    aaTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTc                                                         <  1:312958/149‑1 (MQ=255)
                                                                                      tttCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTTACCTGGCTAGTCAt                                                       >  1:280378/1‑149 (MQ=255)
                                                                                        tCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAACGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCg                                                     <  1:167107/149‑1 (MQ=255)
                                                                                           gATCGACACCTCCGGCGCTATCGCCCAGCTGGTTTATTTAGTCATCGCCCTTTCACAGTTGCGGAAGCGCAAAATCCTGCGAGCAGAAGGAAGCGCAATTCGCTTGCGCATGTGGCTTTGTCCGTGGATCACCTGGCTAGTCATcccct                                                  >  2:172931/1‑145 (MQ=255)
                                                                                            aTCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTATTCATCGGCtt                                                 >  2:215078/1‑149 (MQ=255)
                                                                                              cGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGGTtta                                               >  1:138872/1‑149 (MQ=255)
                                                                                              cGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTtta                                               <  1:496995/149‑1 (MQ=255)
                                                                                              cGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTtta                                               >  2:627667/1‑149 (MQ=255)
                                                                                               gACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTttat                                              >  1:838635/1‑149 (MQ=255)
                                                                                                aCAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCCCCCGGCTGGTCATCGGCTttatt                                             >  1:4299/1‑149 (MQ=255)
                                                                                                aCAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATACGTGGCTCACCTGGCTAGTCATCGGGTttatt                                             >  2:681676/1‑149 (MQ=255)
                                                                                                aCAGCTCCGGCGCTATCGCCCTGCTAGTTTATTTAGTCATCGCCGTATCAAAGTTGCGGATGCGCACCATTCTGCGAGCAGAAGCAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCACACCTGGCTAGTCATCGGATttatt                                             >  1:944678/1‑149 (MQ=255)
                                                                                                    cTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTtta                                               >  2:59174/1‑143 (MQ=255)
                                                                                                    cTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTtta                                               <  1:59174/143‑1 (MQ=255)
                                                                                                    cTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCt                                         <  2:411576/149‑1 (MQ=255)
                                                                                                             tATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTATTCATCGGCTTTATTACCTTTGTGTTgg                                <  1:215078/149‑1 (MQ=255)
                                                                                                                  cccTGCTGGTTTATTTAGTCATCGCCGTTGCAAAGTTGCGGATGCGCAAAGTTCCGCCAGCGGAAGGAAGAGAAATTCGCTTGCGAGTGTGGCTTTATCCGTGCCTCACCTGGCTAGTCATCGGCCTTATTACCTTTGTGTTggcggtg                           >  1:884166/1‑144 (MQ=255)
                                                                                                                   cccgctgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGa                          <  1:713316/146‑1 (MQ=255)
                                                                                                                   ccTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGa                          >  1:584755/1‑149 (MQ=255)
                                                                                                                     tgctgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTTGTAGTGATg                        >  2:163835/1‑149 (MQ=255)
                                                                                                                      gctgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGc                       >  2:550989/1‑149 (MQ=255)
                                                                                                                               tttaGTCATCGCCGTTTAAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTAGCTTTGTGGTGGTAGTGATGGTATtccgtc              >  1:156127/1‑149 (MQ=255)
                                                                                                                                ttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGCATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTCTtccgtgc             >  1:484722/1‑147 (MQ=255)
                                                                                                                                      aTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCg                                                     >  2:242287/1‑103 (MQ=255)
                                                                                                                                      aTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCg                                                     <  1:242287/103‑1 (MQ=255)
                                                                                                                                      aTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCa       >  2:668897/1‑149 (MQ=255)
                                                                                                                                       tCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTGTTCCGTCCGgcgcca      >  1:878409/1‑147 (MQ=255)
                                                                                                                                       tCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCac      >  2:943589/1‑148 (MQ=255)
                                                                                                                                       tCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGCAGTGATGGTATTCCCTACGGCGGaa      >  1:456452/1‑149 (MQ=255)
                                                                                                                                         gCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGAAGAAGGACGCGAAATTCGCTTGCGAATGTGTCTTTATCCGTGGCTCACCTGGCTAGTACTCGACGTTATTCGCTTTGAGTTGGTCATAATGCTTTTCGATCTCGCGCaatc    >  2:214373/1‑147 (MQ=255)
                                                                                                                                          ccGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTg                           <  1:431814/125‑1 (MQ=255)
                                                                                                                                          ccGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTg                           >  2:431814/1‑125 (MQ=255)
                                                                                                                                           cGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGt  <  1:556492/149‑1 (MQ=255)
                                                                                                                                           cGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGt  <  2:838635/149‑1 (MQ=255)
                                                                                                                                             tttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAg                             >  2:389670/1‑120 (MQ=255)
                                                                                                                                             tttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAg                             <  1:389670/120‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CCCCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGT  >  NC_000913_3_vch_tpiA/2795250‑2795537

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