BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,777,050	C→G	5.9%	W166C (TGG→TGC)	PP_2429 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,777,050	0	C	G	5.9%	317.6 / 10.0	104	W166C (TGG→TGC)	PP_2429	membrane protein
Reads supporting (aligned to +/- strand): ref base C (55/40); new base G (2/4); total (59/45)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 105 total aligned reads displayed.

GCGAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTT  >  NC_002947/2776904‑2777199

gcgAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCAc                                                                                                                                                      >  2:149094/1‑150 (MQ=255)

gcgAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGGTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCAc                                                                                                                                                      >  1:1563009/1‑150 (MQ=255)

   aGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTTCCTGCAGCCACCgg                                                                                                                                                   <  1:859634/150‑1 (MQ=255)

      cGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAACCCCCCGACAACACACGGTTGAGCACCCGCAACCGCCCCGGTGGCTGCAGGCACCggggc                                                                                                                                                <  1:1307072/150‑3 (MQ=255)

      cGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGGCTGCAGGCACCggggc                                                                                                                                                <  1:1312018/150‑3 (MQ=255)

      cGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGgcg                                                                                                                                                <  1:282743/150‑1 (MQ=255)

         tGACACGTTGAGCCGAGGCCGCTCAGGACAATTGTCGGTCGATGACAGGCTCAGGCCAATACTGCAGCCCATGGCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTc                                                                                                                                             <  1:917160/150‑1 (MQ=255)

           acacGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTccc                                                                                                                                           <  1:1471734/150‑1 (MQ=255)

            cacGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTcccc                                                                                                                                          >  1:1880010/1‑150 (MQ=255)

                    gCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCaa                                                                                                                                  >  2:917772/1‑150 (MQ=255)

                     ccGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAac                                                                                                                                 >  2:388111/1‑150 (MQ=255)

                         ggCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCg                                                                                                                             <  1:1686837/150‑1 (MQ=255)

                           ccGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGcc                                                                                                                           <  1:388111/150‑1 (MQ=255)

                              cTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAg                                                                                                                        >  2:1832795/1‑150 (MQ=255)

                              cTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAg                                                                                                                        <  1:176050/150‑1 (MQ=255)

                               tCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAgg                                                                                                                       >  1:1482892/1‑150 (MQ=255)

                                   gCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCAccc                                                                                                                                      >  1:64344/1‑131 (MQ=255)

                                   gCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCAccc                                                                                                                                      <  2:64344/131‑1 (MQ=255)

                                   gCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAggcggc                                                                                                                   <  2:1880010/150‑1 (MQ=255)

                                     cccTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCa                                                                                                                 >  2:1111586/1‑150 (MQ=255)

                                      ccTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCggg                                                                                                                                                  >  2:1560948/1‑116 (MQ=255)

                                      ccTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCggg                                                                                                                                                  <  1:1560948/116‑1 (MQ=255)

                                      ccTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCAt                                                                                                                >  1:1085772/1‑150 (MQ=255)

                                          gTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑c                                                                                                          >  1:1300524/1‑150 (MQ=255)

                                            cGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑Cgg                                                                                                        <  1:1776735/150‑1 (MQ=255)

                                              gTCAAGTAGAGGCTCAGGCCAATACTGCNGCCCATGCNGATATCAACNAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCa                                                                                                      >  1:562057/1‑150 (MQ=255)

                                              gTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCa                                                                                                      <  1:1832795/150‑1 (MQ=255)

                                              gTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCa                                                                                                      <  2:1635323/150‑1 (MQ=255)

                                              gTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGCCGACTTCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCa                                                                                                      >  1:1070493/1‑150 (MQ=255)

                                               tCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCCACACCGCCCAGGCGGCCATGCT‑‑CGGCag                                                                                                     >  2:694658/1‑150 (MQ=255)

                                                 aaGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCagag                                                                                                   <  2:1770846/150‑1 (MQ=255)

                                                 aaGTAGAGGCTAAGGACAATAATGAAGCCTATGACGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCagag                                                                                                   <  2:1070493/150‑1 (MQ=255)

                                                  aGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATCTCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGc                                                                                                  <  1:17184/150‑1 (MQ=255)

                                                   gTAGAGGCTGAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCg                                                                                                 >  2:1348704/1‑150 (MQ=255)

                                                      gagGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGcc                                                                                                                           >  2:578632/1‑123 (MQ=255)

                                                      gagGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGcc                                                                                                                           <  1:578632/123‑1 (MQ=255)

                                                         gCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGc                                                                                           >  1:836842/1‑150 (MQ=255)

                                                         gCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGc                                                                                           >  2:1007790/1‑150 (MQ=255)

                                                         gCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGc                                                                                           >  2:1875219/1‑150 (MQ=255)

                                                         gCTCACGCCAATACTGTAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTCAACACCCGCAACCGCCCCGGTGCCGGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGagccgc                                                                                           >  1:1899665/1‑145 (MQ=255)

                                                           tCAGGCCAATACTGCAGCGCATGCCGATATCAACAAAGCCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGAtaaacga                                                                                         <  1:1036770/150‑7 (MQ=255)

                                                              ggCCAATACTGCCGCCCATGCCGATATCAACACAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGaa                                                                                      <  2:562057/150‑1 (MQ=255)

                                                                 cAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAacac                                                                                   >  2:1386433/1‑150 (MQ=255)

                                                                 cAATAATGCAGCACATCCGGAGATCGAAGAAACGCGCGCCAACATTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACCCCGCCCAGGCCGCCCAGCT‑‑CGGCAGACCGCTCAGCAAGAacac                                                                                   >  2:610879/1‑150 (MQ=255)

                                                                  aaTACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACAcc                                                                                  <  2:393659/150‑1 (MQ=255)

                                                                     aCTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGCACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAgc                                                                               <  2:1300524/150‑1 (MQ=255)

                                                                       tGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAgcgc                                                                             >  2:1186022/1‑150 (MQ=255)

                                                                         cAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCa                                                                           >  2:388611/1‑150 (MQ=255)

                                                                          aGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCa                                                                                             <  2:1845304/131‑1 (MQ=255)

                                                                          aGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCa                                                                                             >  1:1845304/1‑131 (MQ=255)

                                                                           gCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAgc                                                                               >  1:1792314/1‑144 (MQ=255)

                                                                           gCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAgc                                                                               <  2:1792314/144‑1 (MQ=255)

                                                                            cccATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCt                                                                        >  1:948616/1‑150 (MQ=255)

                                                                            cccATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCG‑CCCGGTGGCTGCAGGCAc                                                                                                                                                      <  2:1576027/73‑1 (MQ=38)

                                                                                   cGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCacac                                                                 >  1:277346/1‑150 (MQ=255)

                                                                                     atatCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGc                                                                                           >  2:557249/1‑122 (MQ=255)

                                                                                     atatCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGc                                                                                           <  1:557249/122‑1 (MQ=255)

                                                                                      tatCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTgg                                                              >  2:1599127/1‑150 (MQ=255)

                                                                                       atCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGc                                                             >  1:1902709/1‑150 (MQ=255)

                                                                                          aacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAac                                                                                     >  1:719763/1‑123 (MQ=255)

                                                                                          aacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAac                                                                                     <  2:719763/123‑1 (MQ=255)

                                                                                          aacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAac                                                          >  1:1105725/1‑150 (MQ=255)

                                                                                           acaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAacg                                                         >  1:1337908/1‑150 (MQ=255)

                                                                                            caaAACCCCCGCCCACACTTGGTTGAACACCAGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCGCACTGGCAacga                                                        <  2:1070041/150‑1 (MQ=255)

                                                                                            caaAACCCCCGCCAACACTTGGTTTAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGCGCGTCGAGAACACCAGCGCCAGCTGCCGCACTGGCAacga                                                        >  1:1070041/1‑150 (MQ=255)

                                                                                             aaaaCCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAacgac                                                       >  1:914840/1‑150 (MQ=255)

                                                                                                aCCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑AGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGgc                                                    >  1:25200/1‑150 (MQ=255)

                                                                                                 cccccGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGgcg                                                   <  1:149094/150‑1 (MQ=255)

                                                                                                  ccccGCCAACACTTGGTTGAACACCCGCACCCGACCCCCCTCTTTTTCTGGCCGGGTGCTCCCCACCCCCAACAACTCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGgcgc                                                  <  1:917772/150‑1 (MQ=255)

                                                                                                  ccccGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCCGAGAGAGGGAACGCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCCGCACGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGgcgc                                                  >  1:1505836/1‑150 (MQ=255)

                                                                                                  ccccGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCTTTGCTTACTGC‑CAGC‑ATGATGATTAATGCGATCGCCAGCTGCCACACTGGCAACGACGgcgc                                                  <  1:529728/150‑1 (MQ=255)

                                                                                                  ccccGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGgcgc                                                  >  2:584884/1‑150 (MQ=255)

                                                                                                  ccccGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGgcgc                                                  >  2:692536/1‑150 (MQ=255)

                                                                                                      gCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCa                                              >  2:736522/1‑150 (MQ=255)

                                                                                                              ttGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCagag                                                                                                   <  1:1155341/89‑1 (MQ=255)

                                                                                                              ttGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAacgacg                                                      >  1:1834084/1‑134 (MQ=255)

                                                                                                              ttGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAacgacg                                                      <  2:1834084/134‑1 (MQ=255)

                                                                                                              ttGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATAcccc                                      >  2:1264636/1‑150 (MQ=255)

                                                                                                              ttGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATAcccc                                      <  1:1186022/150‑1 (MQ=255)

                                                                                                                 gTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGac                                   <  2:1085772/150‑1 (MQ=255)

                                                                                                                       cacCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGc                             >  2:1438673/1‑150 (MQ=255)

                                                                                                                        acCCGCAACCGCCCCGGTGCCTGCAGTCACCGGGCGCTCACCACGCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCCATCAGCAAGAACACCCGCCCCAGCTGCCACACTGGCAACGACGGCGCGGCACATACCCCGACCgccgcc                            >  2:276436/1‑150 (MQ=255)

                                                                                                                          ccGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCAc                          >  1:1175080/1‑150 (MQ=255)

                                                                                                                          ccGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCAc                          >  2:1584904/1‑150 (MQ=255)

                                                                                                                             cAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGc                       <  2:1375453/150‑1 (MQ=255)

                                                                                                                               aCCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGCAGACACGGGCCAACTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCa                     >  1:717917/1‑150 (MQ=255)

                                                                                                                                 cGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCACGCCACCGCCCAGGCGGCCATGCT‑‑CTGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATc                   >  2:1036842/1‑150 (MQ=255)

                                                                                                                                    cccGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAc                <  2:596913/150‑1 (MQ=255)

                                                                                                                                    cccGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAc                <  1:1584904/150‑1 (MQ=255)

                                                                                                                                    cccGGGGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAc                <  2:1314350/150‑1 (MQ=255)

                                                                                                                                         tGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGAc           <  1:388611/150‑1 (MQ=255)

                                                                                                                                           ccTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCCGCGCATGCTGTAAGTCGTGCAACAACGACTCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCt         >  1:717861/1‑150 (MQ=255)

                                                                                                                                             tGCAGGCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCa                     <  1:1862383/136‑1 (MQ=255)

                                                                                                                                             tGCAGCCACCGGGCGCTCCCCACCCCGAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAGATACCCCGACCACCGCCACGGCCATCAACCAGACCTTg       >  1:1558257/1‑150 (MQ=255)

                                                                                                                                             tGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGGGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTg       >  1:1080594/1‑150 (MQ=255)

                                                                                                                                             tGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCtag       >  2:582069/1‑148 (MQ=255)

                                                                                                                                              gCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTgg      <  2:1482892/150‑1 (MQ=255)

                                                                                                                                               cAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTggg     <  1:692536/150‑1 (MQ=255)

                                                                                                                                                aGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTgggg    <  2:277346/150‑1 (MQ=255)

                                                                                                                                                  ccACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGtt  <  1:1599127/150‑1 (MQ=255)

GCGAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCT‑‑CGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTT  >  NC_002947/2776904‑2777199

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