BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,361,724	C→G	6.5%	R718P (CGG→CCG)	PP_2076 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,361,724	0	C	G	6.5%	337.9 / 11.0	108	R718P (CGG→CCG)	PP_2076	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (52/49); new base G (4/3); total (56/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

Only 100 of 119 total aligned reads displayed.

TGGCTGCTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTCCAGCAG  >  NC_002947/2361577‑2361872

tGGCTGCTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcg                                                                                                                                                      >  3:124205/1‑149 (MQ=255)

   ctgctTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcgcgc                                                                                                                                                   <  4:124205/149‑1 (MQ=255)

   ctgctTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcgcgc                                                                                                                                                   <  5:32976/149‑1 (MQ=255)

   ctgctTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcgcgc                                                                                                                                                   <  7:239033/149‑1 (MQ=255)

     gctTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGa                                                                                                                                                 >  6:15506/1‑149 (MQ=255)

           gTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGacca‑cc                                                                                                                                           <  1:317250/148‑1 (MQ=255)

           gTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCt                                                                                                                                          >  2:111135/1‑149 (MQ=255)

                 cTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCca                                                                                                                                    <  6:172320/149‑1 (MQ=255)

                    aaTCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCcaca                                                                                                                                  >  7:106276/1‑148 (MQ=255)

                    aaTCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCcaca                                                                                                                                  >  2:122711/1‑148 (MQ=255)

                      tCGAGTCAGTGAAATCCCACCACGCTAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCTGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATaa                                                                                                                               <  6:273002/149‑1 (MQ=255)

                             aGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTc                                                                                                                        >  6:35831/1‑149 (MQ=255)

                             aGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTc                                                                                                                        >  6:74925/1‑149 (MQ=255)

                              gTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcgcg                                                                                                                                                    >  7:194164/1‑121 (MQ=255)

                              gTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcgcg                                                                                                                                                    <  8:194164/121‑1 (MQ=255)

                                 aaaTCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCaa                                                                                                                     <  5:15506/148‑1 (MQ=255)

                                    tCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTc                                                                                                                 <  1:122711/149‑1 (MQ=255)

                                      ccaccaCGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCg                                                                                                                >  4:322004/1‑148 (MQ=255)

                                       caccacGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGaa                                                                                                              <  3:219338/149‑1 (MQ=255)

                                       aaccacGCGAGTACCAACCCCCCGCTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTAGTCCGGTCACACTGTTCGCGCCGAACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACATCg                                                                                                                >  2:77616/2‑147 (MQ=255)

                                        accacGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCg                                                                                                                >  7:69686/1‑146 (MQ=255)

                                        accacGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGaa                                                                                                              <  8:69686/148‑1 (MQ=255)

                                          cacGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGTGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAAATGGCTCCAACTTCGAAGcc                                                                                                           >  2:242505/1‑149 (MQ=255)

                                           acGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGccc                                                                                                          >  5:23043/1‑149 (MQ=255)

                                             gcgAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcgcgc                                                                                                                                                   >  6:39254/1‑107 (MQ=255)

                                             gcgAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACcgcgc                                                                                                                                                   <  5:39254/107‑1 (MQ=255)

                                             gcgAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGacca‑cc                                                                                                                                           <  1:122494/114‑1 (MQ=255)

                                             gcgAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGacca‑cc                                                                                                                                           >  2:122494/1‑114 (MQ=255)

                                                  tACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAg                                                                                                   >  2:38394/1‑149 (MQ=255)

                                                  tACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAg                                                                                                   >  1:231455/1‑149 (MQ=255)

                                                    ccAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGcc                                                                                                 <  7:5529/149‑1 (MQ=255)

                                                     cAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGTGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCg                                                                                                >  6:319483/1‑149 (MQ=255)

                                                     cAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGccc                                                                                                          <  4:67245/139‑1 (MQ=255)

                                                     cAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGccc                                                                                                          >  3:67245/1‑139 (MQ=255)

                                                      aaCCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGaa                                                                                                              >  1:232665/1‑134 (MQ=255)

                                                      aaCCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGaa                                                                                                              <  2:232665/134‑1 (MQ=255)

                                                       aCCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTa                                                                                              <  2:231455/149‑1 (MQ=255)

                                                        ccccccGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGCAGCCCGCCGCAGCCGTAg                                                                                             >  5:128939/1‑149 (MQ=255)

                                                          ccccGGTGCCATCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGt                                                                                            <  1:77616/149‑1 (MQ=255)

                                                          ccccGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGt                                                                                            >  8:257553/1‑148 (MQ=255)

                                                          ccccGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTACAACTTCGAAGCCCGCCGCAGCCGTAGTc                                                                                           <  1:193193/149‑1 (MQ=255)

                                                           cccGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCTACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCa                                                                                          <  8:183134/149‑1 (MQ=255)

                                                           cccGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTc                                                                                           >  2:299064/1‑148 (MQ=255)

                                                                   cAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGgcc                                                                                  <  7:75257/149‑1 (MQ=37)

                                                                    aGCTAAACTGGCGCCCCATAGTTTCAAACGACCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTAACGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGgccg                                                                                 >  5:115211/1‑149 (MQ=255)

                                                                      cTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGaa                                                                                                              <  8:262456/118‑1 (MQ=39)

                                                                      cTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCccgc                                                                                                        <  6:270852/124‑1 (MQ=38)

                                                                         aaCTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAg                                                                                             <  1:325402/132‑1 (MQ=37)

                                                                         aaCTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAg                                                                                             >  2:325402/1‑132 (MQ=37)

                                                                         aaCTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAg                                                                             <  1:111135/148‑1 (MQ=37)

                                                                          aCTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAg                                                                             >  4:313058/1‑147 (MQ=37)

                                                                               cgcCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCg                                                                      <  7:257553/149‑1 (MQ=39)

                                                                                gcCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCgg                                                                     <  7:87361/149‑1 (MQ=39)

                                                                                 ccccATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGt                                                                    <  3:313058/149‑1 (MQ=39)

                                                                                    cATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACtt                                                                                                                  <  3:282946/100‑1 (MQ=40)

                                                                                     aTAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGTGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTc                                                                <  5:319483/149‑1 (MQ=255)

                                                                                     aTAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGgccgcc                                                                               >  7:189065/1‑134 (MQ=37)

                                                                                     aTAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGgccgcc                                                                               <  8:189065/134‑1 (MQ=37)

                                                                                     aTAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTc                                                                <  5:35831/149‑1 (MQ=39)

                                                                                        gTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCAc                                                              <  5:74925/148‑1 (MQ=39)

                                                                                         tttCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCAccc                                                            >  8:50159/1‑149 (MQ=39)

                                                                                           tCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCg                                                                      <  4:223229/137‑1 (MQ=40)

                                                                                           tCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCcgc                                                          >  1:1589/1‑149 (MQ=39)

                                                                                               aCGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCAt                                                      >  7:208819/1‑149 (MQ=39)

                                                                                                 gCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCa                                                    >  4:76560/1‑149 (MQ=39)

                                                                                                   ccTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACt                                                  >  7:159175/1‑149 (MQ=39)

                                                                                                     tAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGt                                                >  6:268185/1‑149 (MQ=39)

                                                                                                      agagCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGt                                                >  4:177564/1‑148 (MQ=38)

                                                                                                       gagCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTc                                               <  2:1589/148‑1 (MQ=38)

                                                                                                       gagCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCa                                              <  1:299064/149‑1 (MQ=38)

                                                                                                       gagCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCa                                              <  8:159175/149‑1 (MQ=38)

                                                                                                       gagCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCa                                              <  2:5282/149‑1 (MQ=38)

                                                                                                        agCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAg                                             <  3:177564/149‑1 (MQ=37)

                                                                                                         gCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCAcc                                                             >  4:134443/1‑132 (MQ=37)

                                                                                                         gCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCAcc                                                             <  3:139302/132‑1 (MQ=37)

                                                                                                           ggCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCa                                          <  3:76560/149‑1 (MQ=37)

                                                                                                                      ctcGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACt                                                  >  3:144017/1‑130 (MQ=38)

                                                                                                                      ctcGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACt                                                  <  4:144017/130‑1 (MQ=38)

                                                                                                                       tcGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTcgcg                                   >  6:17238/1‑144 (MQ=39)

                                                                                                                       tcGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTcgcg                                   <  5:17238/144‑1 (MQ=39)

                                                                                                                            cGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGCCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGaa                          <  8:208819/148‑1 (MQ=255)

                                                                                                                               tCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGc                       >  3:95831/1‑148 (MQ=39)

                                                                                                                                 acacTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCAt                    >  8:4331/1‑149 (MQ=39)

                                                                                                                                  cacTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATc                   <  3:322004/149‑1 (MQ=39)

                                                                                                                                   acTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCa                  <  5:321839/149‑1 (MQ=39)

                                                                                                                                   acTGTTAGCTCCGGAGCGCGCGACCA‑CCTGGTCCACATAACTGTCTCCAACTTCGAAGCACGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCTCGCATCACTGTCAGTCAGTTCGCTCGCTTTGAATGCCATCa                  <  4:50340/149‑1 (MQ=14)

                                                                                                                                    gtggtCGCGCGGTACGGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAg                 >  4:215629/5‑149 (MQ=255)

                                                                                                                                    gtggtCGCGCGGTACGGCGCGAACA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAg                 >  3:281666/5‑149 (MQ=255)

                                                                                                                                     tggtCGCGCGGTACGGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGaa                          <  3:191448/136‑1 (MQ=255)

                                                                                                                                     tggtCGCGCGGTACGGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCAGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGaa                          >  4:191448/4‑139 (MQ=255)

                                                                                                                                       ttCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTcgc                                    <  3:91362/127‑1 (MQ=40)

                                                                                                                                         cgcgCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATcc            >  6:319563/1‑149 (MQ=39)

                                                                                                                                          gcgcCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCGGCATccc           >  7:71660/1‑149 (MQ=255)

                                                                                                                                            gcCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCt         >  1:201709/1‑149 (MQ=39)

                                                                                                                                               gTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTCca      <  5:319563/149‑1 (MQ=39)

                                                                                                                                               gTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTCca      <  6:128939/149‑1 (MQ=39)

                                                                                                                                                tACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTCcag     >  2:126517/1‑149 (MQ=39)

                                                                                                                                                  ccGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTCcagca   <  2:201709/149‑1 (MQ=39)

                                                                                                                                                   cgcgcgACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTCcagcag  >  7:250653/1‑149 (MQ=39)

                                                                                                                                                   cgcgcgACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTCcagcag  >  3:294075/1‑149 (MQ=39)

TGGCTGCTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCA‑CCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCCCCTCCAGCAG  >  NC_002947/2361577‑2361872

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