BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,297,313	G→C	8.2%	Q535H (CAG→CAC)	PP_2893 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,297,313	0	G	C	8.2%	759.4 / 18.1	258	Q535H (CAG→CAC)	PP_2893	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (97/136); new base C (23/2); total (120/138)
Fisher's exact test for biased strand distribution p-value = 7.65e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

ATGCTCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGA  >  NC_002947/3297164‑3297462

    tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCCGCtgg                                                                                                                                                   >  1:5361450/1‑149 (MQ=255)

    tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCtgt                                                                                                                                                   <  2:4402395/150‑1 (MQ=255)

    tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCtgt                                                                                                                                                   >  2:167714/1‑150 (MQ=255)

    tCAAGGCAGAGACCCACGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCAGGGTGACAACCGCATCAACTGCAGTGCCAGCATGCAACACCGCCTGGCTGGGCGAATCGACCTTGACTAGCCGCGCCACGGCCAGcaga                                                                                                                                                   >  1:540655/1‑147 (MQ=255)

      aaGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTgg                                                                                                                                                 <  2:2759324/150‑1 (MQ=255)

       aGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCACGGCCAGCGGTGGc                                                                                                                                                >  1:5015706/1‑150 (MQ=255)

          ccGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGa                                                                                                                                             <  1:3883773/150‑1 (MQ=255)

           cGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAg                                                                                                                                            <  2:3260268/150‑1 (MQ=255)

           cGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCACAGCTGGCCGAg                                                                                                                                            <  1:2014699/150‑1 (MQ=255)

           cGAGACCCAGGGGGCAGGGCAAAACTGGACCCCGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCCGCAGTGGCCGAg                                                                                                                                            >  1:545837/1‑150 (MQ=255)

            gagaCCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCCGCTGTGGCCGAgg                                                                                                                                           >  2:578878/1‑150 (MQ=255)

            gagaCCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAACTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCCGCTGTGGCCGCgg                                                                                                                                           >  1:2237229/1‑150 (MQ=255)

               aCCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTa                                                                                                                                        <  1:857904/150‑1 (MQ=255)

                cccAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTAc                                                                                                                                       <  1:4671576/150‑1 (MQ=255)

                 ccAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACa                                                                                                                                      >  2:5039091/1‑150 (MQ=255)

                  cAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACaa                                                                                                                                     <  2:875580/150‑1 (MQ=255)

                  cAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACaa                                                                                                                                     >  2:4926340/1‑150 (MQ=255)

                   aGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCGGTGGCCGAGGCTACAAg                                                                                                                                    >  2:1121155/1‑150 (MQ=255)

                     ggggCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGc                                                                                                                                  <  1:516796/150‑1 (MQ=255)

                     ggggCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCGGTGGCCGAGGCTACAAGGc                                                                                                                                  >  2:3106496/1‑150 (MQ=255)

                        gcagggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAACGCCgg                                                                                                                               >  1:3562954/1‑150 (MQ=255)

                            ggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCa                                                                                                                           >  1:2740330/1‑150 (MQ=255)

                            ggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCa                                                                                                                           >  2:4097795/1‑150 (MQ=255)

                            ggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCGGTGGCCGAGGCTACAAGGCCGGGTCa                                                                                                                           >  2:2178938/1‑150 (MQ=255)

                             gcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTATAAGGCCGGGTCaa                                                                                                                          <  1:2912647/150‑1 (MQ=255)

                              caAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAg                                                                                                                         <  2:4173617/150‑1 (MQ=255)

                              caAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAg                                                                                                                         >  2:144565/1‑150 (MQ=255)

                               aaaaCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCGGTGGCCGAGGCGACAAGGCCGGGTCCAgg                                                                                                                        >  1:242651/1‑150 (MQ=255)

                               aaaaCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCACGGCCAGCAGTGGCCGAGGCTGCAAGGCCgggggaacg                                                                                                                        >  1:3411408/1‑144 (MQ=255)

                                   cTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCGGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCg                                                                                                                    >  2:2691799/1‑150 (MQ=255)

                                   cTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCCGCTGTTGCCCAGGCGACAAAGCCGGGTCAAAGGCCg                                                                                                                    >  1:2980234/1‑150 (MQ=255)

                                     ggACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGc                                                                                                                  <  1:4506071/150‑1 (MQ=255)

                                      gACCCGGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCt                                                                                                                 <  1:4583178/150‑1 (MQ=255)

                                       aCCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTg                                                                                                                >  1:5550865/1‑150 (MQ=255)

                                         ccTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGa                                                                                                              <  1:4197547/150‑1 (MQ=255)

                                          cTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGAt                                                                                                             <  2:3541021/150‑1 (MQ=255)

                                              gCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAAGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCg                                                                                                         >  1:5574418/1‑150 (MQ=255)

                                                 cacTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCg                                                                                                      <  1:451238/150‑1 (MQ=255)

                                                 cacTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCg                                                                                                      <  1:451373/150‑1 (MQ=255)

                                                  acTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCCGCTGTGGCCGAGGCGACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCgg                                                                                                     >  1:5066508/1‑150 (MQ=255)

                                                  acTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCgg                                                                                                     <  1:4016031/150‑1 (MQ=255)

                                                  acTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCCGCGGGGGGCGCGGCGCCAACGCCGGGTCAAGGGCCGGCGGGGTGTGGCCgg                                                                                                     >  1:1379434/1‑150 (MQ=255)

                                                    tGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGCTGTCGCCGGCa                                                                                                   >  2:777553/1‑150 (MQ=255)

                                                      gCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCCCCGGGGGCCGCGGCGACAAAGCCCGGGCAAGGGCCGGCGGGGTGTCGCCgggccc                                                                                                 >  1:159086/1‑146 (MQ=255)

                                                         aTCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCACGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTAGCCGGCAccaag                                                                                              >  2:1083398/1‑147 (MQ=255)

                                                           ccAGCTGGGTGACACCCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTCGCTGGGCGAAGCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCa                                                                                            <  2:5168325/150‑1 (MQ=255)

                                                           ccAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCa                                                                                            <  2:3858980/150‑1 (MQ=255)

                                                           ccAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCa                                                                                            <  1:4011142/150‑1 (MQ=255)

                                                           ccAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGACACCTTGCa                                                                                            <  2:78520/150‑1 (MQ=255)

                                                            catcTGGGTTAAAACCGCATCAACGGCAGTGGCATCTTTCAACAGCGCCTGGCTGGGCGAATCGACCTTGTCTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAg                                                                                           <  2:4805485/147‑1 (MQ=255)

                                                              gCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCACCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGc                                                                                         >  2:3391661/1‑150 (MQ=255)

                                                                tGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCACCTGTGGCCCAGGCGACAACGCCCGGGCAACGGCCGGCGGGATGTCGCCGGCACCATGCAGGccc                                                                                       >  1:3264980/1‑150 (MQ=255)

                                                                tGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCACGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGccc                                                                                       >  2:2752389/1‑150 (MQ=255)

                                                                 gggTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGcccc                                                                                      <  2:1202946/150‑1 (MQ=255)

                                                                      aCAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTAGGCCAGAGGTGGCCGAGGCTACAAGGCCGGGTCCAGGGCAGGCGGGATGTAGCCGGCACCATGCAGGCCCCCCAAg                                                                                 >  2:2004565/1‑150 (MQ=255)

                                                                        aaCCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAgg                                                                                                                                           >  2:444319/1‑90 (MQ=255)

                                                                        aaCCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAgg                                                                                                                                           <  1:444319/90‑1 (MQ=255)

                                                                        aaCCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGt                                                                               >  2:4554764/1‑150 (MQ=255)

                                                                        aaCCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGt                                                                               <  2:3633539/150‑1 (MQ=255)

                                                                          ccGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCCGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCGGGATGTCGCCGGCACCTTGCAGGCCCCCCCAGGTAc                                                                             >  2:2424277/1‑150 (MQ=255)

                                                                          ccGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCGGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCGGGATGTAGCCGGCACCTTGCAGGCCCCGCCACGTAc                                                                             >  1:4026473/1‑150 (MQ=255)

                                                                          ccGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCGGTGGCCGAGGCGACAAGGCCGGGTCAAGGGCCGGCGGGATGTCGCCGGCCCCGTGCAGGCCCCCCAAGGTAc                                                                             >  1:2710831/1‑150 (MQ=255)

                                                                           cGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACg                                                                            >  2:5053597/1‑150 (MQ=255)

                                                                           cGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCCGCTGGATGTCGCCGGCACCGTGCAGGCCCCGCAAGGTACg                                                                            >  2:349875/1‑150 (MQ=255)

                                                                                     gCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGc                                                                  >  2:5128129/1‑150 (MQ=255)

                                                                                      cAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCa                                                                 <  2:4229016/150‑1 (MQ=255)

                                                                                       aGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAg                                                                <  2:1288073/150‑1 (MQ=255)

                                                                                             aGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCACCGGTGGCCGGGGCGACAAGGCCCGGTCAAGGGCCGGCTGGATGTCGCCGGCACCATTCAGGCCCCGCAAGGTACGCTCACCCGGCAGGGCCac                                                          >  1:358158/1‑149 (MQ=255)

                                                                                                   cAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTg                                                    <  1:3364483/150‑1 (MQ=255)

                                                                                                   cAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTg                                                    <  1:2345906/150‑1 (MQ=255)

                                                                                                    aaCAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTgg                                                   <  1:4356256/150‑1 (MQ=255)

                                                                                                    aaCAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCGGTGGCCCAGGCTACAAGGCCGGGGCAAGGGCCGGCTGGATGTCGCCGGCACCTTTCAGGCCCCGCAAGGTACGCTACCCCAGCAGGGCCCCCGCTTgg                                                   >  2:1916813/1‑150 (MQ=255)

                                                                                                     aCAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGc                                                  >  2:5504391/1‑150 (MQ=255)

                                                                                                       aGCGCCTGTCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCTGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGccc                                                <  1:1857151/150‑1 (MQ=255)

                                                                                                        gcgcCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCa                                               <  2:1270250/150‑1 (MQ=255)

                                                                                                         cgcCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGCTGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAg                                              >  2:3939628/1‑150 (MQ=255)

                                                                                                          gcCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAgg                                             <  1:2938389/150‑1 (MQ=255)

                                                                                                           ccTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGGCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGc                                            >  2:4082781/1‑150 (MQ=255)

                                                                                                           ccTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGc                                            <  1:1297273/150‑1 (MQ=255)

                                                                                                           ccTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCGGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAACGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGc                                            >  2:4535758/1‑150 (MQ=255)

                                                                                                                ctggGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGaaa                                       <  1:2439465/150‑1 (MQ=255)

                                                                                                                   ggCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAAcc                                    <  2:5405995/150‑1 (MQ=255)

                                                                                                                      gAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGcc                                 <  2:1053946/150‑1 (MQ=255)

                                                                                                                      gAATCGACCTTGACTTGCCGCGCCTCGGCCAGCGGTGGGCGAGGCTACAAGGCCGGGTCAATGGCCGGCTGGATGTCGCCGGCACCGTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGCAAACCGcc                                 >  2:1613815/1‑150 (MQ=255)

                                                                                                                        aTCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCT‑ACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTgc                              <  1:1701418/150‑1 (MQ=255)

                                                                                                                             ccTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTgcagca                          <  2:2924625/150‑1 (MQ=255)

                                                                                                                             ccTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTgcagca                          <  2:113391/150‑1 (MQ=255)

                                                                                                                                 gACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTg                      <  2:2440926/150‑1 (MQ=255)

                                                                                                                                        cgcgCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTgg               <  1:4926340/150‑1 (MQ=255)

                                                                                                                                          agcCTCGGCCACAGCTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGAt             >  2:3682205/2‑150 (MQ=255)

                                                                                                                                            ccTCGGCCACAGCTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGc           >  1:1071580/1‑150 (MQ=255)

                                                                                                                                            ccTCGGCCACAGCTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGc           >  1:4104932/1‑150 (MQ=255)

                                                                                                                                              tCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGccc         >  1:698391/1‑150 (MQ=255)

                                                                                                                                              tCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGccc         >  1:2363275/1‑150 (MQ=255)

                                                                                                                                              tCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGccc         >  1:3520864/1‑150 (MQ=255)

                                                                                                                                                ggCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGc       <  2:313813/150‑1 (MQ=255)

                                                                                                                                                 gCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGcc      <  1:5266095/150‑1 (MQ=255)

                                                                                                                                                   cAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTg    <  1:449705/150‑1 (MQ=255)

                                                                                                                                                     gCTGTGGCCGAGGCTACCAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGa  <  2:5302007/150‑1 (MQ=255)

                                                                                                                                                     gCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGa  <  1:1154224/150‑1 (MQ=255)

ATGCTCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGA  >  NC_002947/3297164‑3297462

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