BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,962,964	T→G	7.2%	T65P (ACG→CCG) ‡	ygdB ←	DUF2509 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,962,964	0	T	G	7.2%	255.2 / 4.4	107	T65P (ACG→CCG) ‡	ygdB	DUF2509 family protein
Reads supporting (aligned to +/- strand): ref base T (56/42); new base G (4/5); total (60/47)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

CAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGC  >  NC_000913/2962817‑2963112

cAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTc                                                                                                                                                     >  2:719837/1‑149 (MQ=255)

       tCCATCAATGACTTCGCCTGTTCGCCACAACGCAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCcagc                                                                                                                                              >  2:1496001/1‑149 (MQ=255)

       tCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGGTTCAGCGTAATGCGAGCACTGAACTGCTGGCTGCGTCTGCcagc                                                                                                                                              >  1:1447342/1‑149 (MQ=255)

        ccATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAATCGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCcagca                                                                                                                                             >  1:673366/1‑149 (MQ=255)

        ccATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCcagca                                                                                                                                             <  1:147565/149‑1 (MQ=255)

         cATCAATGACTTCGCCTGTTCGCCACAACGCAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGAAAACAAACCTGGGCATAGGTTTCAGCGTACTGCGAGCAATGAACTGCTGGCTGCGTCTGAcaggat                                                                                                                                            >  2:1483906/1‑146 (MQ=255)

           tCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCAGACGCAGCCAGCAGtt                                                                                                                                          <  1:440494/149‑2 (MQ=255)

             aaTGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCAGACGCAGCCAGCAGTTCa                                                                                                                                        <  1:1472895/149‑1 (MQ=255)

                gACTTCGCCTGTTCGCCACAACGAAACGCCGTCATAACCGGCACTCGAGAAGGCGTCATTCTCTGCCAGTAAACGCAAACAAACCGGGGCATCGGGTTCAGCGTACTGCGAGCACTGCACTGCTGGCTGCGTCTGCCAGCAGTGCAttt                                                                                                                                     <  1:1125534/149‑1 (MQ=255)

                  cTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCAttttt                                                                                                                                   >  1:63027/1‑149 (MQ=255)

                  cTTCGCCTGTTAGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCTGCAGTGCAatttt                                                                                                                                   >  2:514357/1‑149 (MQ=255)

                   ttCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCGTCATTATCTGCCAGTAAACGCCAACACCCCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTc                                                                                                                                  >  1:1550671/1‑149 (MQ=255)

                   ttCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGCCAATAAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTc                                                                                                                                  <  2:565753/149‑1 (MQ=255)

                     cGCCTGTGCGCCACTACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTccc                                                                                                                                <  1:685700/149‑1 (MQ=255)

                      gCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTc                                                                                                                                  >  2:1552832/1‑146 (MQ=255)

                      gCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTc                                                                                                                                  <  1:1552832/146‑1 (MQ=255)

                          gTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTTTCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGc                                                                                                                           >  2:1438030/1‑149 (MQ=255)

                          gTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAAAAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGc                                                                                                                           <  2:1393898/149‑1 (MQ=255)

                                 acaACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGcc                                                                                                                    <  2:481496/149‑1 (MQ=255)

                                   aaCTAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACACACCTGGGCATGGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGa                                                                                                                  <  1:277219/149‑1 (MQ=255)

                                   aaCGAAACGCCTTCATACCCGGCACTCAATAACGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGTTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGa                                                                                                                  <  2:884051/149‑1 (MQ=255)

                                     cGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACt                                                                                                                <  2:913830/149‑1 (MQ=255)

                                      gAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTg                                                                                                               >  1:1466071/1‑149 (MQ=255)

                                            ccTTCATAACCGGCAATAAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGAt                                                                                                         <  2:433773/149‑1 (MQ=255)

                                             cTTCATCACCGGCACTCCATAAGGCTTCATTATCTGCCAGTCACCGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGCCCGGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATg                                                                                                        <  2:625243/149‑1 (MQ=255)

                                                 aTAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCCGCGCCGACTGAACGATGGCCt                                                                                                    >  2:1034969/1‑149 (MQ=255)

                                                    aCCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCCGCGCCGACTGAACGATGGCCTggc                                                                                                 >  1:831528/1‑149 (MQ=255)

                                                       ggcatacaatAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCTAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTggcggc                                                                                              <  1:1163990/143‑1 (MQ=255)

                                                       ggCAATCAATAAGGCTTCATTATCTGCCTGTAAACGCAAACAACCCTGGGCTTCGCTTTCAGCGTCCTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCCGTGCAGTTTTCCCCACGCCAGCGCGGACTGACCGATGGCCTggcggc                                                                                              >  1:128765/1‑149 (MQ=255)

                                                       ggCAATCAATAAGGCTTCATTATCTGCCACTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGAGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGACTcgcggc                                                                                              >  2:473715/1‑149 (MQ=255)

                                                        gCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGTGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGgcggcg                                                                                             >  1:566826/1‑149 (MQ=255)

                                                         caatcaatAAGGCTTCATTATCTGCCAGTAAACGCAACCAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGcc                                                                                                                    <  1:785809/125‑1 (MQ=255)

                                                         caatcaatAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGcc                                                                                                                    >  2:785809/1‑125 (MQ=255)

                                                         caatcaatAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGgcgc                                                                                            <  1:627543/149‑1 (MQ=255)

                                                          aatcaatATGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACAGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGTCTGAACGATGGCCTGGCGGCGCa                                                                                           >  2:411583/1‑149 (MQ=255)

                                                              aatAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCTTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAtga                                                                                       >  1:1263271/1‑149 (MQ=255)

                                                               atAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAtgac                                                                                      >  1:427029/1‑149 (MQ=255)

                                                                tAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAtgact                                                                                     >  1:916180/1‑149 (MQ=255)

                                                                tAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAtgact                                                                                     >  1:1214606/1‑149 (MQ=255)

                                                                    gCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATgactgact                                                                                 >  1:285354/1‑149 (MQ=255)

                                                                        cATTATCTGCAAGTAAACGCAAACAACCCTGGGCAGCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTcc                                                                             <  2:1359771/149‑1 (MQ=255)

                                                                          ttATCTGCCAGTAAACGCAAACAACCCTGGGCCTCGGTCTCCGCGTACTGCGAGCACTGAACTGCTGGCTGCGGCGGCCAGCAGTGCATTGTTCCCCACGCCAGCGCCGACTGAACGATGGCCGGGCGGCGCAATGACTGACTTTCCAt                                                                           <  1:1393846/149‑1 (MQ=255)

                                                                             tcagCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCt                                                                        <  1:195209/146‑1 (MQ=255)

                                                                                gCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGAACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTACCCCACGCCACCGCCGACTGCCCGAAGGCCTGGCGGCCCAATGACAGACTTTACATGCacac                                                                     >  1:1025227/1‑149 (MQ=255)

                                                                                   aGTAAACGCAAACAAACCTGGGCATCGGTTTCCGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGCCTGACTTTCCATGCTCAcgcg                                                                  >  1:1459854/1‑149 (MQ=255)

                                                                                    gTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCga                                                                 <  1:522292/149‑1 (MQ=255)

                                                                                    gTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTACCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAAGACTGACTTTCCATGCTCAcgcca                                                                 >  1:1307762/1‑147 (MQ=255)

                                                                                     tAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCgag                                                                <  1:901573/149‑1 (MQ=255)

                                                                                       aaCGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGagag                                                              <  1:1319240/149‑1 (MQ=255)

                                                                                        aCGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCTATGACTGACTTTCCATGCTCACGCGAGAgg                                                             >  1:1339934/1‑149 (MQ=255)

                                                                                          gCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCCCGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCGCACGCGAGAGGCa                                                           >  1:1568771/1‑149 (MQ=255)

                                                                                            aaacaagcCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCaaa                                                         >  2:700031/1‑149 (MQ=255)

                                                                                            aaacaaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCaaa                                                         >  2:449598/1‑149 (MQ=255)

                                                                                            aaacaaacCTGGGCATCGGATTCAGCGTACTGCGAGCACTGAACAGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACAGAACGATGGCCTGGCGGCGCAATGCCTGCATTTCCATGCTCACGCGAGAGGCaaa                                                         >  2:1534416/1‑149 (MQ=255)

                                                                                             aacaaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCaaaa                                                        >  1:1063076/1‑149 (MQ=255)

                                                                                             aacaaaaGTGGGCATCGGTTTCAGCGTACTGCGAGCTCTGAACTGCTGGCTGCGTCTTCCAGCAGTGAATTTTTCACCACGCCAGCGCCTTCTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCaaaa                                                        <  1:569825/149‑1 (MQ=255)

                                                                                              acaaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAAc                                                       >  1:41133/1‑149 (MQ=255)

                                                                                              acaaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAAc                                                       >  1:610881/1‑149 (MQ=255)

                                                                                               caaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAACCt                                                      >  2:477042/1‑149 (MQ=255)

                                                                                                    cTGGGCATCGGTTTCAGTGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGAt                                                 <  2:566826/149‑1 (MQ=255)

                                                                                                    cTGGGCATCGGTTGCAGCGTACTGCGAGCACTGAACTGCTGGATGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGAt                                                 <  2:890806/149‑1 (MQ=255)

                                                                                                         cATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCAGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCATTGACTGACTTTCCATGCTCAAGCGAGAGGCAAAACTGCGCTCCTGt                                            >  2:1241295/1‑149 (MQ=255)

                                                                                                          aTCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGtt                                           <  2:603207/149‑1 (MQ=255)

                                                                                                          aTCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGtt                                           <  2:982930/149‑1 (MQ=255)

                                                                                                          aTCGGTTGCCGCGTACTGCGAGCACTGACCTGCTGGCTGCGTCAGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGtt                                           <  2:935712/149‑1 (MQ=255)

                                                                                                            cGGGTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTgg                                         <  2:10996/149‑1 (MQ=255)

                                                                                                                  cAGCGTACTGCGAGCACTGAACTGCTGTCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCACTGACTGCCTTTCCATGCACACGCGAGAGGCAAAACGTCGATACTGTTGTCTCAtt                                   >  1:1064395/1‑149 (MQ=255)

                                                                                                                   aGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGGTCATTc                                  >  1:66474/1‑149 (MQ=255)

                                                                                                                   aGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTc                                  <  1:638617/149‑1 (MQ=255)

                                                                                                                   aGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCGCATTc                                  >  1:624145/1‑149 (MQ=255)

                                                                                                                    gCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATTGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTcc                                 <  1:582545/149‑1 (MQ=255)

                                                                                                                    gCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTcc                                 <  1:1438030/149‑1 (MQ=255)

                                                                                                                     cGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCCCGCCAGCGCCGACTGACCGATGGCGTGGCGGCGCAATGAGTGACTTACCTTGCGAACGCGAGAGGCAAAACTGCGATCCTGTTGGCTAATTCCt                                >  1:16155/1‑149 (MQ=255)

                                                                                                                       tACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCCGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTg                              <  2:1144604/149‑1 (MQ=255)

                                                                                                                       tACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTg                              <  1:1488480/149‑1 (MQ=255)

                                                                                                                        aCTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGt                             >  2:786210/1‑149 (MQ=255)

                                                                                                                        aCTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGt                             <  1:668604/149‑1 (MQ=255)

                                                                                                                         cTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTa                            <  1:547620/149‑1 (MQ=255)

                                                                                                                          tGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTATCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCAAGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTACTTGTaa                           <  2:1447342/149‑1 (MQ=255)

                                                                                                                           gCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCGTGTTGGCTCATTCCTTGTAAt                          >  1:505824/1‑149 (MQ=255)

                                                                                                                            cGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCATTGACTGACTTTCCATGCTCACGCGAGAGGCAAACCTGCGATCCTGTTGGCTCATTCCTTGTAATa                         <  1:899990/149‑1 (MQ=255)

                                                                                                                                 aCTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCaaa                    <  1:268388/149‑1 (MQ=255)

                                                                                                                                 aCTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGAATTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCaaa                    >  1:362147/1‑149 (MQ=255)

                                                                                                                                 aCTGAACTGCTGGCTGCCTCAGCCAGCAGTGCATTTTTCCCCAAGCCAGAGCCGACTGAACGATTGCCTGCCGGCGCAATGAATGCCTATCCATGATCACGCGATAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGaaaa                    >  2:494796/1‑149 (MQ=255)

                                                                                                                                 aCTGAACTGCTGGCGGCGGCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGGTGGCTCATTCCTTGTAATAGCaaa                    <  1:486195/149‑1 (MQ=255)

                                                                                                                                   tGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGACGACTGAACGATGGCCTGGCGGCGCAATTACTGACTATCCATGATCAAGCTAGAGGCAAAACTGCGCTCGTGTTGGCTCATGCCTTGAAATAGCaaagt                  >  2:869131/1‑147 (MQ=255)

                                                                                                                                    gcaCTGCTGGCAGACGCAGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAATGATGGCCTGGCGCCGCTATGACTGACTTTCCATGCTCACGCGAGCGGCATCACTGCGAGCCTGTTGGCTCTTTCCTTGTAATAGGAAACTa                 >  1:178714/3‑149 (MQ=255)

                                                                                                                                    gcaCTGCTGGCAGACGCAGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAAGGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATGCCGTGTAATTGCCAACTa                 >  1:184229/3‑149 (MQ=255)

                                                                                                                                    gcaCTGCTGGCAGACGCAGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTGGATATAGCACACTa                 >  2:340150/3‑149 (MQ=255)

                                                                                                                                     aaCTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTAc                <  1:1190034/149‑1 (MQ=255)

                                                                                                                                     aaCTGCTGGCTGCGTCTGACAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCATTGCCTGACGTTCCATGCTGCCGCGAGAGGCAAAACTGCCATCCTCTAGGCTCATTCCTTGTAATAGCAAACTAc                >  2:303701/1‑149 (MQ=255)

                                                                                                                                      aCTGCTGGCTGCGGCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTAcc               <  1:1483906/149‑1 (MQ=255)

                                                                                                                                        tgctgGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTAcacc             >  1:915984/1‑146 (MQ=255)

                                                                                                                                        tgctgGCTGCGTCTGCCAGCAGTGCATATATCCCCACGCCCGCGCCGACTGAGCGATGGCCTGGCGGCGCAATGACTGACTTTCCATCCACACGCGAGAGGCAAAACTGCGATCCTGGTGGCTCAATCCTTGTAAAAGCAAACTACCCa             >  1:36597/1‑149 (MQ=255)

                                                                                                                                          gtggctgCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCaaa           <  2:1214606/148‑1 (MQ=255)

                                                                                                                                          ctggctgCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCaaa           <  2:1466071/149‑1 (MQ=255)

                                                                                                                                            ggctgcggctgcCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAAc         <  1:440567/149‑1 (MQ=255)

                                                                                                                                             gctgCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATCGCCTGGCGGCGCAATGCGTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCACAAcc        >  1:964687/1‑149 (MQ=255)

                                                                                                                                               tgCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACcag      <  1:719837/149‑1 (MQ=255)

                                                                                                                                                   tCTGCCAGCCGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCTATGACTGACTTTCCCTGCTCACGCGAGAGGCAAAACTCCGATCCTGTAGGCTCAGTCCTTGTAATAGCACACTACCCAAACCCagcagc  >  1:1027695/1‑149 (MQ=255)

CAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGC  >  NC_000913/2962817‑2963112

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`