BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,962,966	T→A	5.7%	Q64L (CAG→CTG)	ygdB ←	DUF2509 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,962,966	0	T	A	5.7%	290.1 / 2.8	107	Q64L (CAG→CTG)	ygdB	DUF2509 family protein
Reads supporting (aligned to +/- strand): ref base T (56/43); new base A (4/3); total (60/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

TCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGC  >  NC_000913/2962824‑2963112

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)

      aaTGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCa                                                                                                                                        >  1:374254/1‑149 (MQ=255)

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

       aTGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTACACGCAAACACACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGACCTGCTGGCTGCGTCTGCCAGCAGAGCAt                                                                                                                                       >  1:799360/1‑149 (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)

                       gCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTACCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAgc                                                                                                                       >  1:603207/1‑149 (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)

                                         cATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGcc                                                                                                     >  1:1552798/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)

                                                 gCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATGGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGGCTGCCCGCAGTGCATTTTTCCCCACGCCACCGCCGCCTGAACGCTGGCCAggcggcc                                                                                             >  2:388522/1‑148 (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:1214606/1‑149 (MQ=255)

                                                         tAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAtgact                                                                                     >  1:916180/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)

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

                                                                                     aaacaaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCaaa                                                         >  2:449598/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)

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

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

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

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

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

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

                                                                                                         ttCAGTGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCCTGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCa                                     >  2:1258160/1‑149 (MQ=255)

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

                                                                                                            aGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGGTCATTc                                  >  1:66474/1‑149 (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)

                                                                                                                tACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTg                              >  1:1334/1‑149 (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)

                                                                                                                                 tgctgGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCa             >  2:613734/1‑149 (MQ=255)

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

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

                                                                                                                                   ctggctgCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCaaa           <  2:1466071/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)

TCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGC  >  NC_000913/2962824‑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: `‑`