BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,031,923	T→A	6.0%	L5* (TTG→TAG) ‡	PP_0889 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,031,923	0	T	A	6.0%	349.7 / 9.7	114	L5* (TTG→TAG) ‡	PP_0889	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (48/56); new base A (4/3); total (55/59)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

ATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTC  >  NC_002947/1031776‑1032062

atatCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGc                                                                                                                                           >  1:596230/1‑150 (MQ=255)

acatCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGc                                                                                                                                           >  1:1510687/3‑150 (MQ=255)

   tctcGGCCACCGCACGGGCGCTGGGCATGCACAGACGGCCGTGCCTTCGCAGTGTGTAGAAGCGGCCGGTTCGCTGTTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCcg                                                                                                                                        <  1:435322/150‑1 (MQ=255)

     tcGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAACCCGGGCT‑GGCCcgcg                                                                                                                                      >  2:354618/1‑150 (MQ=255)

       ggCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGaa                                                                                                                                    <  1:898345/150‑1 (MQ=255)

        gCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAg                                                                                                                                   <  1:75812/150‑1 (MQ=255)

            ccGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCt                                                                                                                               >  1:1706995/1‑150 (MQ=255)

                 cGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCTCGAAGGGCTGCaaa                                                                                                                          <  2:1734882/150‑1 (MQ=255)

                    gcgcTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCa                                                                                                                       <  2:1571211/150‑1 (MQ=255)

                          ggCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGCTGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCa                                                                                                                 >  2:92782/1‑150 (MQ=255)

                           gCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCaa                                                                                                                <  2:1706995/150‑1 (MQ=255)

                           gCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCaa                                                                                                                <  2:596230/150‑1 (MQ=255)

                                cACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAg                                                                                                           >  2:1410177/1‑150 (MQ=255)

                                cACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAg                                                                                                           <  2:326107/150‑1 (MQ=255)

                                  ccGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGc                                                                                                         >  1:402756/1‑150 (MQ=255)

                                  ccGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCA‑AGCCcgcgctcaca                                                                                                                                >  2:1104407/1‑121 (MQ=255)

                                  ccGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCA‑AGCCcgcgctcaca                                                                                                                                <  1:1104407/127‑7 (MQ=255)

                                   cGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCg                                                                                                        >  1:719012/1‑150 (MQ=255)

                                    gacggacgTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGa                                                                                                       >  1:532860/1‑150 (MQ=255)

                                       ggacgTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCt                                                                                                    <  2:781361/150‑1 (MQ=255)

                                           gTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGAt                                                                                                <  2:92068/150‑1 (MQ=255)

                                           gTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGAt                                                                                                <  1:1715323/150‑1 (MQ=255)

                                              gCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTTTTGCAACTCTGTG‑GCCGCCTGGGTTGGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGaa                                                                                             >  1:717500/1‑150 (MQ=255)

                                                 gcgcAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGa                                                                                                 >  2:1700786/1‑143 (MQ=255)

                                                 gcgcAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGa                                                                                                 <  1:1700786/143‑1 (MQ=255)

                                                       gCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAg                                                                                            >  1:192751/1‑142 (MQ=255)

                                                       gCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAg                                                                                            <  2:192751/142‑1 (MQ=255)

                                                       gCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGtt                                                                                    >  1:1246256/1‑150 (MQ=255)

                                                       gCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCA‑AGCCcgcgc                                                                                                                                     >  1:1734028/1‑100 (MQ=255)

                                                       gCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCA‑AGCCcgcgc                                                                                                                                     <  2:1734028/101‑2 (MQ=255)

                                                         tGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTgg                                                                                  <  2:1510687/150‑1 (MQ=255)

                                                          gCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCA‑AGCCcgcgct                                                                                                                                    >  2:943396/1‑97 (MQ=255)

                                                          gCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCA‑AGCCcgcgct                                                                                                                                    <  1:943396/99‑3 (MQ=255)

                                                                  ggCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTg                                                                         <  2:1246256/150‑1 (MQ=255)

                                                                   gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgc                                                                        <  2:35027/150‑1 (MQ=255)

                                                                   gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgc                                                                        >  1:1860713/1‑150 (MQ=255)

                                                                     cGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgcg                                                                       >  2:1080548/1‑149 (MQ=255)

                                                                     cGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgcg                                                                       <  1:1080548/149‑1 (MQ=255)

                                                                       gTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTg                                                                    >  1:684075/1‑150 (MQ=255)

                                                                          cGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGtt                                                                 <  2:402756/150‑1 (MQ=255)

                                                                           ggCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTAATCCTTGCGCTAGtgt                                                                >  2:1379197/1‑148 (MQ=255)

                                                                             cgcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCa                                                              <  2:1753769/150‑1 (MQ=255)

                                                                              gcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                             <  1:102393/150‑1 (MQ=255)

                                                                              gcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                             <  2:684075/150‑1 (MQ=255)

                                                                               cTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCACTCTCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAgg                                                            <  1:726871/150‑1 (MQ=255)

                                                                                 gATTCCTTATCGATGAGGTTGGCCTGTTCGCGGGGTAATCCGCGCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTa                                                          >  1:1341050/1‑150 (MQ=255)

                                                                                  aTTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGGAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAg                                                         >  1:1794969/1‑150 (MQ=255)

                                                                                  aTTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATcc                                                                            >  2:1225022/1‑131 (MQ=255)

                                                                                  aTTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATcc                                                                            <  1:1225022/131‑1 (MQ=255)

                                                                                     ccTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAgcgc                                                      <  1:1870304/150‑1 (MQ=255)

                                                                                      cTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGtt                                                                                    >  2:650033/1‑119 (MQ=255)

                                                                                      cTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTc                                                                                <  1:342906/123‑1 (MQ=255)

                                                                                      cTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTc                                                                                >  2:342906/1‑123 (MQ=255)

                                                                                           cGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGACT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGc                                                <  2:38955/150‑1 (MQ=255)

                                                                                            gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCt                                                                                                                               <  1:484093/70‑1 (MQ=255)

                                                                                            gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCt                                                                                                                               >  2:484093/1‑70 (MQ=255)

                                                                                            gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCtt                                                                          >  2:1043673/1‑123 (MQ=255)

                                                                                            gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCtt                                                                          <  1:1043673/123‑1 (MQ=255)

                                                                                            gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTACTCCAACTATGACCCTTACCC‑GGTCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCt                                               >  1:530133/1‑150 (MQ=255)

                                                                                             atgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTg                                              >  1:1667052/1‑150 (MQ=255)

                                                                                              tgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGtt                                                                 <  2:1561680/130‑1 (MQ=255)

                                                                                              tgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGtt                                                                 >  1:1561680/1‑130 (MQ=255)

                                                                                               gatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGCTGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGaa                                            <  1:354618/150‑1 (MQ=255)

                                                                                                  gTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATc                                                                                                              >  1:1898661/1‑81 (MQ=255)

                                                                                                  gTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATc                                                                                                              <  2:1898661/81‑1 (MQ=255)

                                                                                                   ttGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACt                                        <  2:530133/150‑1 (MQ=255)

                                                                                                      gCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                     <  2:719012/150‑1 (MQ=255)

                                                                                                      gCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                     <  2:1860713/150‑1 (MQ=255)

                                                                                                       ccTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTgg                                                                   <  1:110154/119‑1 (MQ=255)

                                                                                                       ccTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTgg                                                                   >  2:110154/1‑119 (MQ=255)

                                                                                                        cTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCt                                   >  2:805281/1‑150 (MQ=255)

                                                                                                          gTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGa                                 <  1:1410177/150‑1 (MQ=255)

                                                                                                                gggggAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                             >  2:870398/1‑116 (MQ=255)

                                                                                                                gggggAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                             <  1:870398/116‑1 (MQ=255)

                                                                                                                    gAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCCGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTgg                       >  1:242575/1‑150 (MQ=255)

                                                                                                                    gAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgcgc                                                                      <  2:1909921/103‑1 (MQ=255)

                                                                                                                    gAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgcgc                                                                      >  1:1909921/1‑103 (MQ=255)

                                                                                                                     agtCCGCTCCTGCAACTGTGAGGGCGGGCT‑TGCCCGCGAAGGACTTCAAAGCAGTGCCAATCAAGGCGAGCTCGATGAAGGTACTGTTTGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCGGACCCGCAGTGGc                      <  2:242575/148‑1 (MQ=255)

                                                                                                                      aTCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCa                     <  1:909896/150‑1 (MQ=255)

                                                                                                                       tCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAg                    >  2:421488/1‑150 (MQ=255)

                                                                                                                          gCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgc                                                                        <  2:1123783/95‑1 (MQ=255)

                                                                                                                          gCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgc                                                                        >  1:1123783/1‑95 (MQ=255)

                                                                                                                          gCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGa                                          >  2:929776/1‑125 (MQ=255)

                                                                                                                          gCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGa                                          <  1:929776/125‑1 (MQ=255)

                                                                                                                          gCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAAc                 >  2:173256/1‑150 (MQ=255)

                                                                                                                            tCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGa               <  1:92782/150‑1 (MQ=255)

                                                                                                                             ccTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGAc              >  1:883188/1‑150 (MQ=255)

                                                                                                                              cTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGAccc                              >  1:1707715/1‑133 (MQ=255)

                                                                                                                              cTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGAccc                              <  2:1707715/133‑1 (MQ=255)

                                                                                                                              cTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGAcc                               <  2:116585/132‑1 (MQ=255)

                                                                                                                              cTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGAcc                               >  1:116585/1‑132 (MQ=255)

                                                                                                                                gCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGa                                          <  1:1118120/119‑1 (MQ=255)

                                                                                                                                gCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGa                                          >  2:1118120/1‑119 (MQ=255)

                                                                                                                                  aaCTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGc         >  2:115645/1‑150 (MQ=255)

                                                                                                                                  aaCTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGc         <  2:717500/150‑1 (MQ=255)

                                                                                                                                  aaCTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGc         >  1:1353607/1‑150 (MQ=255)

                                                                                                                                      cccttcgcgGGCA‑AGCCCGCGCAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCAcc     >  2:594022/6‑150 (MQ=255)

                                                                                                                                        gAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACt                                                                                       <  2:714297/66‑1 (MQ=255)

                                                                                                                                        gAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACt                                                                                       >  1:714297/1‑66 (MQ=255)

                                                                                                                                         aGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTc  <  1:805281/150‑1 (MQ=255)

ATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCT‑TGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTC  >  NC_002947/1031776‑1032062

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