Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 388,155 T→G S121A (TCG→GCG)  tauD → taurine dioxygenase, 2‑oxoglutarate‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913388,1550TG93.3% 510.7 / 18.1 165S121A (TCG→GCG) tauDtaurine dioxygenase, 2‑oxoglutarate‑dependent
Reads supporting (aligned to +/- strand):  ref base T (3/8);  new base G (87/67);  total (90/75)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 165 total aligned reads displayed.

TAACGATAATCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGC  >  NC_000913/388058‑388251
                                                                                                 |                                                                                                
tAACGATAATCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGGTTCTGGCAGCTAAAGAGTTACCTGCg                                                                                                >  2:2407198/1‑100 (MQ=255)
tAACGATAATCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACACCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCg                                                                                                <  2:2109983/100‑1 (MQ=255)
  aCGATAATCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGAc                                                                                              <  1:1636952/100‑1 (MQ=255)
  aCGATAATCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGAc                                                                                              <  2:3112527/100‑1 (MQ=255)
   cGATAATCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGAcc                                                                                             <  1:3389975/100‑1 (MQ=255)
   cGATAATCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGAcc                                                                                             >  1:4075790/1‑100 (MQ=255)
   cGATAATCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGAcc                                                                                             >  1:651035/1‑100 (MQ=255)
   cGATAATCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGAcc                                                                                             >  1:1959249/1‑100 (MQ=255)
   cGATAATCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGAcc                                                                                             >  2:3178554/1‑100 (MQ=255)
   cGATAATCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGAcc                                                                                             >  2:621749/1‑100 (MQ=255)
    taTAATCCGCCAGATAACGACAACCGGCATACCGATGTGACATTTATTGAAACGCCGCCCTCAGGGGCGATTCTTGTAGCTAAAGAGTTATCTTCGACcg                                                                                            <  2:666559/99‑1 (MQ=255)
    gATAATCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACcg                                                                                            >  2:3512728/1‑100 (MQ=255)
     aTAATCCGTCAGATAACGAGAAATGGCATACCGATGTGACATCTAGTGAAACGCCTACCGCCGGGGCGATTCTGGCAGCTAAAGAGTGACCTGCGACcgg                                                                                           <  2:1540621/100‑1 (MQ=255)
       tctCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACcggcg                                                                                         <  1:1762113/98‑1 (MQ=255)
        aTCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACcggcgg                                                                                        <  2:1926845/100‑1 (MQ=255)
         tCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGt                                                                                       <  2:1346046/100‑1 (MQ=255)
         tCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGt                                                                                       <  1:3837983/100‑1 (MQ=255)
          ccGGCAGATAACGACAACTGGCATACCGATGTCACATTTATTGAAACGCCACACGCCCCGCCCCTCCTGGCATCTAAAGAGTTACCTGCCACAGGCGGTg                                                                                      >  2:64128/1‑100 (MQ=255)
          ccGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACcggcgaga                                                                                      >  2:1733419/1‑97 (MQ=255)
          ccGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTg                                                                                      <  2:486280/100‑1 (MQ=255)
          ccGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTg                                                                                      <  2:1106123/100‑1 (MQ=255)
          ccGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTg                                                                                      <  1:1128989/100‑1 (MQ=255)
             ccAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATa                                                                                   <  2:2428100/100‑1 (MQ=255)
             ccAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATa                                                                                   <  1:3983851/100‑1 (MQ=255)
                  tAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGctc                                                                              >  2:3520589/1‑100 (MQ=255)
                    aCGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTg                                                                            >  1:4001362/1‑100 (MQ=255)
                     cGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTgg                                                                           >  1:3579301/1‑100 (MQ=255)
                     cGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTgg                                                                           >  2:922251/1‑100 (MQ=255)
                     cGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTgg                                                                           >  1:2518654/1‑100 (MQ=255)
                       aCAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGAc                                                                         >  1:851655/1‑100 (MQ=255)
                           cTGGCATACCGATGTGACATTTATTGAAACGCTACCCGCAGGGGGGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGc                                                                     <  1:3742233/100‑1 (MQ=255)
                           cTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGc                                                                     >  2:3482980/1‑100 (MQ=255)
                           cTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGc                                                                     <  2:2179718/100‑1 (MQ=255)
                           cTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGc                                                                     >  1:2789027/1‑100 (MQ=255)
                             ggCATACCGATGTGACATTTATTGAAACGCCACCCGCGGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCgg                                                                   >  2:2377626/1‑100 (MQ=255)
                             ggCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCgg                                                                   <  1:3892626/100‑1 (MQ=255)
                             ggCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCgg                                                                   >  1:1360036/1‑100 (MQ=255)
                             ggCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCgg                                                                   >  2:1646726/1‑100 (MQ=255)
                             ggCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCgg                                                                   <  1:2422324/100‑1 (MQ=255)
                               cATACCGATGTGACATTTATTGAAACGCCACCCGCTGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTa                                                                 >  2:3897908/1‑100 (MQ=255)
                                aTACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTAt                                                                <  1:573726/100‑1 (MQ=255)
                                aTACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTAt                                                                >  2:1666453/1‑100 (MQ=255)
                                 tACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTAtt                                                               >  2:3070468/1‑100 (MQ=255)
                                   ccGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGc                                                             >  2:401625/1‑100 (MQ=255)
                                     gATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCgg                                                           >  1:1798060/1‑100 (MQ=255)
                                       tgtgACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGcc                                                         <  1:2876516/100‑1 (MQ=255)
                                        gtgACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCt                                                        >  2:1536642/1‑100 (MQ=255)
                                        gtgACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCt                                                        <  2:2903118/100‑1 (MQ=255)
                                         tgACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTa                                                       <  2:1941037/100‑1 (MQ=255)
                                           aCATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATg                                                     <  1:2018783/100‑1 (MQ=255)
                                           aCATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATg                                                     <  1:2225443/100‑1 (MQ=255)
                                            cATTTATTGAAACGCCACCCGCAGGGGCGAGTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGa                                                    <  1:251526/100‑1 (MQ=255)
                                             aTTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAg                                                   >  1:2886704/1‑100 (MQ=255)
                                              tttATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAgg                                                  <  1:2536495/100‑1 (MQ=255)
                                              tttATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAgg                                                  <  1:2353541/100‑1 (MQ=255)
                                              tttATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAgg                                                  >  1:273415/1‑100 (MQ=255)
                                               ttATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGgc                                                 <  2:577038/100‑1 (MQ=255)
                                                tATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGgcg                                                >  1:2756633/1‑100 (MQ=255)
                                                  ttGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGct                                              >  1:1849123/1‑100 (MQ=255)
                                                   tGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGctc                                             >  1:679449/1‑100 (MQ=255)
                                                   tGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGctc                                             >  1:3393201/1‑100 (MQ=255)
                                                    gAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGctct                                            >  2:886757/1‑100 (MQ=255)
                                                    gAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGctct                                            >  1:2522337/1‑100 (MQ=255)
                                                     aaaCGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCCATGAGGCGctctc                                           >  2:3263060/1‑100 (MQ=255)
                                                      aaCGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGctctct                                          >  1:387073/1‑100 (MQ=255)
                                                        cGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGt                                        <  2:1986730/100‑1 (MQ=255)
                                                         gccgccCGCAGGGGCGATTCTGGCAGCTAGAGGGTTACCTGCGACCGGCGGTGAGACGCTCTAGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGtt                                       <  2:3962192/100‑1 (MQ=255)
                                                          ccACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTc                                      <  1:861492/100‑1 (MQ=255)
                                                          ccACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTc                                      <  1:511830/100‑1 (MQ=255)
                                                          ccACCCGCAGGGGCGATTCTGGCAGCGAAAGAGGTACCTGCGACCGGCGGGGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCCCGCTGTTc                                      >  2:2002598/1‑100 (MQ=255)
                                                            aCCCGCAGGGGCGATTATGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTccc                                    >  2:3877884/1‑100 (MQ=255)
                                                             cccGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCt                                   >  1:375464/1‑100 (MQ=255)
                                                             cccGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCt                                   >  2:2900378/1‑100 (MQ=255)
                                                              ccGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCtt                                  <  2:883723/100‑1 (MQ=255)
                                                               cGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTc                                 <  1:2057686/100‑1 (MQ=255)
                                                                   ggggCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGcc                             <  2:2870403/100‑1 (MQ=255)
                                                                      gCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAgc                          <  1:609677/100‑1 (MQ=255)
                                                                      gCGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAgc                          <  2:1335619/100‑1 (MQ=255)
                                                                       cGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAgct                         >  1:2093493/1‑100 (MQ=255)
                                                                       cGATTCTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAgct                         >  1:1892741/1‑100 (MQ=255)
                                                                            cTGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGa                    >  2:2940682/1‑100 (MQ=255)
                                                                             tGGCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAg                   <  1:2633124/100‑1 (MQ=255)
                                                                              ggCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGt                  >  1:674036/1‑100 (MQ=255)
                                                                              ggCAGCTAAAGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGt                  <  1:1838319/100‑1 (MQ=255)
                                                                                     aaaGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCtgc           >  1:1019806/1‑100 (MQ=255)
                                                                                     aaaGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCtgc           >  1:3527713/1‑100 (MQ=255)
                                                                                      aaGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCtgcg          >  1:1973358/1‑100 (MQ=255)
                                                                                      aaGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCtgcg          >  1:3777100/1‑100 (MQ=255)
                                                                                      aaGAGTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCtgcg          >  1:3097394/1‑100 (MQ=255)
                                                                                        gagTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCtgcgtg        >  2:3609441/1‑100 (MQ=255)
                                                                                         agTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCtgcgtgc       <  1:1390058/100‑1 (MQ=255)
                                                                                          gTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCtgcgtgcg      >  2:2826465/1‑100 (MQ=255)
                                                                                          gTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCtgcgtgcg      >  1:2319446/1‑100 (MQ=255)
                                                                                          gTTACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCtgcgtgcg      >  1:1668873/1‑100 (MQ=255)
                                                                                           ttACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCgg     <  1:2789723/100‑1 (MQ=255)
                                                                                            tACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGa    >  1:3557629/1‑100 (MQ=255)
                                                                                            tACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGa    >  1:3180199/1‑100 (MQ=255)
                                                                                            tACCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGa    >  1:303086/1‑100 (MQ=255)
                                                                                             aCCTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAg   <  2:1486484/100‑1 (MQ=255)
                                                                                              ccTGCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCCGCTGCTGAGTGGGCGGCGTGCGGAGc  >  2:2701336/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TAACGATAATCCGCCAGATAACGACAACTGGCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGC  >  NC_000913/388058‑388251

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