BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	218,130	T→G	S216R (AGC→CGC)	proS ←	prolyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	218,130	0	T	G	99.4%	558.0 / ‑4.4	158	S216R (AGC→CGC)	proS	prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base G (83/73); total (85/73)
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.71e-01

Only 100 of 161 total aligned reads displayed.

GCAGCAGCGCGCGGTTCTTTCGGCGCGATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCTCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGAT  >  NC_000913/218033‑218226

  agcagcGCGCGGTTCTTTCGGCGCGATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTgc                                                                                              >  2:2940818/1‑100 (MQ=255)

    cagcGCGCGGTTCTTTCGGCGCGATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTgcgc                                                                                            <  2:412623/100‑1 (MQ=255)

    cagcGCGCGGTTCTTTCGGCGCGATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTgcgc                                                                                            >  1:2048891/1‑100 (MQ=255)

       cgcgcgGTTCTTTCGGCGCGATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAg                                                                                         <  1:3355544/100‑1 (MQ=255)

         cgcgGTTCTTTCGGCGCGATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCa                                                                                       <  1:2921297/100‑1 (MQ=255)

              ttCTTTCGGCGCGATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTgg                                                                                  <  1:746303/100‑1 (MQ=255)

                cTTTCGGCGCGATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGaa                                                                                <  1:907817/100‑1 (MQ=255)

                cTTTCGGCGCGATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGaa                                                                                >  1:1965726/1‑100 (MQ=255)

                 nttCGGCGCGATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAAt                                                                               >  1:3244249/2‑100 (MQ=255)

                  ttCGGCGCGATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAAtt                                                                              >  1:1050706/1‑100 (MQ=255)

                    cGGCGCGATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCg                                                                            <  1:165709/100‑1 (MQ=255)

                      gcggGATAGCTTCTGCCGGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTg                                                                          <  1:3306704/100‑1 (MQ=255)

                      gcgcGATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGGCAGCACCTGGATTTCGTg                                                                          <  1:3406787/100‑1 (MQ=255)

                      gcgcGATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTg                                                                          <  1:1989991/100‑1 (MQ=255)

                      gcgcGATAGCTTCTGCCAGCTCAATGTTCGCTGCATAGTCAGAGGGGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCTGCACCTGGAATTCGTg                                                                          <  2:3874762/100‑1 (MQ=255)

                        gcgATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTgag                                                                        <  2:2889140/100‑1 (MQ=255)

                        gcgATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTgag                                                                        <  1:359628/100‑1 (MQ=255)

                        gcgATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTgag                                                                        <  1:1265765/100‑1 (MQ=255)

                        gcgATAGCTGCTGCCATTTCAATCTTCGCTTCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTTCGCCAGCACCTGGAATTCGTgag                                                                        <  2:3352493/100‑1 (MQ=255)

                           aTAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAgg                                                                     <  1:1971519/100‑1 (MQ=255)

                           aTAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAgg                                                                     >  2:2353170/1‑100 (MQ=255)

                           aTAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAgg                                                                     <  1:2779209/100‑1 (MQ=255)

                              gCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGgcgc                                                                  >  1:155291/1‑100 (MQ=255)

                              gCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGgcgc                                                                  >  1:1340368/1‑100 (MQ=255)

                              gCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGCCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGgcgc                                                                  <  1:3469634/100‑1 (MQ=255)

                               cTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCt                                                                 >  1:620811/1‑100 (MQ=255)

                                 tCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTgc                                                               >  2:1241520/1‑100 (MQ=255)

                                  cTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTgcc                                                              >  1:3952655/1‑100 (MQ=255)

                                        gTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGAt                                                        >  1:3975242/1‑100 (MQ=255)

                                        gTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGAt                                                        <  1:220372/100‑1 (MQ=255)

                                        gTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGAACTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGAt                                                        >  2:823947/1‑100 (MQ=255)

                                        gTTCAATGTTCGCTGCATAGGCAGAGGTGTCGGAGAAGAGCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGAt                                                        <  2:3294053/100‑1 (MQ=255)

                                         ttCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATa                                                       <  2:2798319/100‑1 (MQ=255)

                                          tCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAg                                                      <  1:44333/100‑1 (MQ=255)

                                               gTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCg                                                 >  1:3935473/1‑100 (MQ=255)

                                               gTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCg                                                 <  2:2117819/100‑1 (MQ=255)

                                                ttCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCgg                                                <  1:988644/100‑1 (MQ=255)

                                                ttCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCgg                                                <  2:1896189/100‑1 (MQ=255)

                                                  cGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGgtg                                              <  1:2676216/100‑1 (MQ=255)

                                                  cGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGgtg                                              <  1:1168251/100‑1 (MQ=255)

                                                   gCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGgtgt                                             >  1:2321657/1‑100 (MQ=255)

                                                   gCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGgtgt                                             >  2:2992158/1‑100 (MQ=255)

                                                   gCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGgtgt                                             >  1:2740718/1‑100 (MQ=255)

                                                    cTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTc                                            >  1:2015110/1‑100 (MQ=255)

                                                     tGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCg                                           >  1:3098025/1‑100 (MQ=255)

                                                      gCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCgg                                          >  2:214914/1‑100 (MQ=255)

                                                      gCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCgg                                          >  2:524814/1‑100 (MQ=255)

                                                      gCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCgg                                          >  2:889695/1‑100 (MQ=255)

                                                       cATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGc                                         >  1:557493/1‑100 (MQ=255)

                                                        aTAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCt                                        >  2:2128165/1‑100 (MQ=255)

                                                        aTAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCt                                        <  1:533943/100‑1 (MQ=255)

                                                        aTAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCt                                        <  1:2341269/100‑1 (MQ=255)

                                                        aTAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGACGCCGATAGAACCGGTGTCGGCt                                        >  2:1382689/1‑100 (MQ=255)

                                                          aGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTg                                      >  2:1678128/1‑100 (MQ=255)

                                                            tCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTa                                    >  2:3240980/1‑100 (MQ=255)

                                                             cAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTAc                                   <  1:1014791/100‑1 (MQ=255)

                                                             cAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTAc                                   >  2:1304800/1‑100 (MQ=255)

                                                             cAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGGAc                                   >  2:3099437/1‑100 (MQ=255)

                                                                 tgtgtCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGgcg                               <  2:756277/99‑1 (MQ=255)

                                                                 ggTGTCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGgcg                               >  1:1085402/1‑100 (MQ=255)

                                                                  gtgtCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGgcgc                              >  1:2951888/1‑100 (MQ=255)

                                                                  gtgtCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGgcga                              >  1:772458/1‑99 (MQ=255)

                                                                    gtCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCgg                            <  2:3392122/100‑1 (MQ=255)

                                                                     tCGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGa                           <  2:1024363/100‑1 (MQ=255)

                                                                      cGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCTATAGAACCGGTGTCGGCTTGTACGGCGCGGaa                          <  2:2777820/100‑1 (MQ=255)

                                                                      cGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGaa                          >  1:3126478/1‑100 (MQ=255)

                                                                      cGGAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGaa                          >  1:3519740/1‑100 (MQ=255)

                                                                       ggAGAAGACCACGTCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGaaa                         >  1:3300790/1‑100 (MQ=255)

                                                                       ggAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGaaa                         <  1:370308/100‑1 (MQ=255)

                                                                       ggAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGaaa                         >  1:476402/1‑100 (MQ=255)

                                                                       ggAGAAGACCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGaaa                         >  2:377362/1‑100 (MQ=255)

                                                                         agaagaCCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATc                       <  2:1410586/100‑1 (MQ=255)

                                                                          gaagaCCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATcc                      <  1:3537755/100‑1 (MQ=255)

                                                                           aagaCCACATCGTCTTCACCGCGCGGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCa                     <  2:905501/100‑1 (MQ=255)

                                                                            agaCCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAg                    <  2:1760093/100‑1 (MQ=255)

                                                                            agaCCACATCGTCTTCACCGCGCTGCGCCAGCACCGGGAATGCGTGAGAGGCGCTGCCGCCGAGAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAg                    <  2:3357338/100‑1 (MQ=255)

                                                                             gaCCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTATGGCGCGGAAATCCAGc                   >  2:255661/1‑100 (MQ=255)

                                                                             gaCCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGc                   <  1:2669604/100‑1 (MQ=255)

                                                                              aCCACATCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGcc                  >  1:4010328/1‑100 (MQ=255)

                                                                                 acaTCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCa               <  2:1623083/100‑1 (MQ=255)

                                                                                  catcatCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCAt              >  2:3378089/1‑100 (MQ=255)

                                                                                   aTCGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATg             <  2:153969/100‑1 (MQ=255)

                                                                                     tctCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCGCATggg           <  1:272168/98‑3 (MQ=255)

                                                                                     cGTCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCg           <  2:3842909/100‑1 (MQ=255)

                                                                                       tCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGTTTGTACGGCGCGGAAATCCAGCCCCATGCGGc         <  2:1731036/100‑1 (MQ=255)

                                                                                       tCTTCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGc         <  2:1590479/100‑1 (MQ=255)

                                                                                         ttCACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTg       >  1:1361078/1‑100 (MQ=255)

                                                                                           ctccGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGaa     <  2:614280/98‑1 (MQ=255)

                                                                                           cacaGCTCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGaa     >  1:3141297/1‑100 (MQ=255)

                                                                                           cACCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGaa     >  2:444840/1‑100 (MQ=255)

                                                                                            aCCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAg    <  2:3425820/100‑1 (MQ=255)

                                                                                            aCCGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAg    <  1:4068472/100‑1 (MQ=255)

                                                                                             ccGCGCTGCGCCAGCACCTGGAATTCGTGAGCGGCGCTGCCGCCGATAGAACCGGTGGCGGCTTGGACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGa   <  2:1563165/100‑1 (MQ=255)

                                                                                             ccGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGTGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGa   >  1:919734/1‑100 (MQ=255)

                                                                                             ccGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCGCATGGGGCTggatt   <  1:2265928/100‑5 (MQ=255)

                                                                                             ccGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCGCATGGGGCTggatt   <  1:3509928/100‑5 (MQ=255)

                                                                                             ccGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCGCATGGGGCTggatt   <  2:84401/100‑5 (MQ=255)

                                                                                             ccGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGa   >  1:1447396/1‑100 (MQ=255)

                                                                                             ccGCGCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGa   >  1:2484583/1‑100 (MQ=255)

                                                                                              cgcgCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGAGCGGAAATAAAGCCCCATAAGGCTGAAgct  >  2:985140/1‑98 (MQ=255)

GCAGCAGCGCGCGGTTCTTTCGGCGCGATAGCTTCTGCCAGTTCAATGTTCGCTGCATAGTCAGAGGTGTCGGAGAAGACCACATCGTCTTCACCGCTCTGCGCCAGCACCTGGAATTCGTGAGAGGCGCTGCCGCCGATAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGAT  >  NC_000913/218033‑218226

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: `‑`