Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,501,720 T→A 8.4% Y53* (TAT→TAA)  tehB → tellurite, selenium methyltransferase, SAM‑dependent; tellurite, selenium resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,501,7200TA 8.4% 149.4 / 2.7 45Y53* (TAT→TAA) tehBtellurite, selenium methyltransferase, SAM‑dependent; tellurite, selenium resistance protein
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base A (2/2);  total (24/21)
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.96e-01

TCTTTACCAATTTTATTATTGCAATACTGCTCATCCGTACTTTTGCGCTTCTGATGCAGGGAAAATTGTTAGTCAGAACCGAGCGCGCCGTTTTAATGAAAGCAGAGGACAAAGAATGATCATTCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCG‑CTGATTATCCATGTAC  >  NC_000913/1501447‑1501992
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                 
tCTTTACCAATTTTATTATTGCAATACTGCTCATCCGTACTTTTGCGCTTCTGATGCAGGGAAAATTGTTAGTCAGAACCGAGCGCGCCGTTTTAATGAAAGCAGAGGACAAAGAATGATCATTCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACCCTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGGATgg                                                                                                                                                                                                                                                                    >  2:282172/1‑289 (MQ=255)
           tttATTATTGCAATACTGCTCATCCGTACTTTTGCGCTTCTGATGCAGGGAAAATTGTTAGTCAGAACCGAGCGCGCCGTTTTAATGAAAGCAGAGGACAAAGAATGATCATTCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATg                                                                                                                                                                                                                                                          >  2:405087/1‑288 (MQ=255)
                             cTCATCCGTACTTTTGCGCTTCTGATGCAGGGAAAATTGTTAGTCAGAACCGAGCGCGCCGTTTTAATGAAAGCAGAGGACAAAGAATGATCATTCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACg                                                                                                                                                                                                                                        >  1:180605/1‑288 (MQ=255)
                             cTCATCCGTACTTTTGCGCTTCTGATGCAGGGAAAATTGTTAGTCAGAACCGAGCGCGCCGTTTTAATGAAAGCAGAGGACAAAGAATGATCATTCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTAAGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACg                                                                                                                                                                                                                                        >  2:459233/1‑288 (MQ=255)
                                                            gAAAATTGTTAGGCAGAACCGAGCGCGCCGTTTTAATGAAAGCAGAGGACAAAGAATGATCATTCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGa                                                                                                                                                                                                         <  1:261886/288‑1 (MQ=255)
                                                                            aaCCGAGCGCGCCGTTTTAATGAAAGCAGAGGACAAAGAATGATCATTCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAg                                                                                                                                                                                        >  2:464490/1‑289 (MQ=255)
                                                                              ccGAGCGCGCCGTTTTAATGAAAGCAGAGGACAAAGAATGATCATTCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTc                                                                                                                                                                                      >  2:332888/1‑289 (MQ=255)
                                                                                 aGCGCGCCGTTTTAATGAAAGCAGAGGACAAAGAATGATCATTCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAg                                                                                                                                                                                                                                   >  1:147371/1‑241 (MQ=255)
                                                                                 aGCGCGCCGTTTTAATGAAAGCAGAGGACAAAGAATGATCATTCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAg                                                                                                                                                                                                                                   <  2:147371/241‑1 (MQ=255)
                                                                                                   aaGCAGAGGACAAAGAATGATCATTCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTcac                                                                                                                                                                  <  1:464490/288‑1 (MQ=255)
                                                                                                    aGCAGAGGACAAAGAATGATCATTCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACAt                                                                                                                                                                >  1:93078/1‑289 (MQ=255)
                                                                                                    aGCAGAGGACAAAGAATGATCATTCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTAACCCCCCGAGTCGTTGATCTGAATAACCTCACAt                                                                                                                                                                >  2:272365/1‑289 (MQ=255)
                                                                                                         agGACAAAGAATGATCATTCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGAt                                                                                                                                                           >  2:241354/1‑289 (MQ=255)
                                                                                                            aCAAAGAATGATCATTCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATa                                                                                                                                                          <  1:285109/287‑1 (MQ=255)
                                                                                                                          ttCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTAt                                                                                                                                          <  1:24518/289‑1 (MQ=255)
                                                                                                                                    aaaCTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACt                                                                                                                                <  2:93078/289‑1 (MQ=255)
                                                                                                                                     aaCTAATTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACtg                                                                                                                               <  2:190772/289‑1 (MQ=255)
                                                                                                                                                  aTAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTcac                                                                                                                                                                  <  2:40825/241‑1 (MQ=255)
                                                                                                                                                  aTAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTcac                                                                                                                                                                  >  1:40825/1‑241 (MQ=255)
                                                                                                                                                  aTAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGAGTGTGGTGCTGATGtt                                                                                                                  >  1:76413/1‑289 (MQ=255)
                                                                                                                                                   tAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGtt                                                                                                                  >  1:388838/1‑288 (MQ=255)
                                                                                                                                                      atatGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGGGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCtt                                                                                                              <  2:277869/289‑1 (MQ=255)
                                                                                                                                                            aTTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCt                                                                                                        >  1:461902/1‑289 (MQ=255)
                                                                                                                                                                 cccGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTaaaa                                                                                                    >  1:159147/1‑288 (MQ=255)
                                                                                                                                                                 cccGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTaaaa                                                                                                    >  2:28472/1‑288 (MQ=255)
                                                                                                                                                                      acacacTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTAAGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATc                                                                                               <  1:459233/288‑1 (MQ=255)
                                                                                                                                                                               gAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGa                                                                                     <  2:461902/289‑1 (MQ=255)
                                                                                                                                                                                aaGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGAt                                                                                    >  1:167680/1‑289 (MQ=255)
                                                                                                                                                                                      tGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTgg                                                                                                                                                                                                          <  2:449659/165‑1 (MQ=255)
                                                                                                                                                                                      tGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTgg                                                                                                                                                                                                          >  1:449659/1‑165 (MQ=255)
                                                                                                                                                                                        gAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAata                                                                            >  1:12623/1‑289 (MQ=255)
                                                                                                                                                                                          aGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATg                                                                          <  1:332888/289‑1 (MQ=255)
                                                                                                                                                                                                gAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACg                                                                     <  1:241354/288‑1 (MQ=255)
                                                                                                                                                                                                       gTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACt                                                              <  2:388838/288‑1 (MQ=255)
                                                                                                                                                                                                       gTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACt                                                              <  2:389239/288‑1 (MQ=255)
                                                                                                                                                                                                       gTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTACCAATATGCAACGTTGCACt                                                              <  2:167680/288‑1 (MQ=255)
                                                                                                                                                                                                                 gTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTTAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCtggtgg                                                   >  1:295180/1‑289 (MQ=255)
                                                                                                                                                                                                                       cGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTAAGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACa                                              >  2:190449/1‑288 (MQ=255)
                                                                                                                                                                                                                          tGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACg                                                                     <  1:96436/262‑1 (MQ=255)
                                                                                                                                                                                                                          tGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACg                                                                     >  2:96436/1‑262 (MQ=255)
                                                                                                                                                                                                                          tGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACaa                                             >  1:298262/1‑286 (MQ=255)
                                                                                                                                                                                                                                                 cGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCg                   <  1:28472/289‑1 (MQ=255)
                                                                                                                                                                                                                                                 cGTAACAGTCTTTACCTGGCAGCCAATGGTTAAGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCg                   <  1:190449/289‑1 (MQ=255)
                                                                                                                                                                                                                                                             tACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCG‑CTGATTATCCAt      <  1:282172/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gCAGCCAATGGCTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATACCCTCACATTTTATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATTGTCTTTGAGGCTATAAACATCCCCGGGTTGATTGCCATTATGAAACGTTGCACTAACCCTGGTGGTTCCAACCTGATTGTGGCGGCAATGGA‑ACCGTCTGCTTATCCATgatc  >  2:69084/1‑284 (MQ=255)
                                                                                                                                                                                                                                                                                 tgatgTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGAt                                                                                    >  2:427627/1‑192 (MQ=255)
                                                                                                                                                                                                                                                                                 tgatgTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGAt                                                                                    <  1:427627/192‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                 
TCTTTACCAATTTTATTATTGCAATACTGCTCATCCGTACTTTTGCGCTTCTGATGCAGGGAAAATTGTTAGTCAGAACCGAGCGCGCCGTTTTAATGAAAGCAGAGGACAAAGAATGATCATTCGTGACGAAAACTATTTTACTGATAAATATGAATTAACCCGCACACACTCTGAAGTACTGGAAGCGGTGAAAGTGGTTAAACCGGGTAAAACGCTGGATCTGGGCTGTGGCAATGGTCGTAACAGTCTTTACCTGGCAGCCAATGGTTATGATGTTGACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCG‑CTGATTATCCATGTAC  >  NC_000913/1501447‑1501992

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑