BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,918,011	T→C	10.1%	V122A (GTC→GCC) ‡	PP_4328 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,918,011	0	T	C	10.1%	318.1 / 18.7	118	V122A (GTC→GCC) ‡	PP_4328	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (50/55); new base C (4/8); total (55/63)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

Only 100 of 118 total aligned reads displayed.

CTTCAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTC  >  NC_002947/4917863‑4918159

cTTCAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCTCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTTGACACCAGCAAGGTGCCGGTc                                                                                                                                                      <  2:424227/150‑1 (MQ=255)

cTTCAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTc                                                                                                                                                      <  2:294916/150‑1 (MQ=255)

 ttCAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACCACAGCAAGGTGCCGGTCg                                                                                                                                                     <  1:695289/150‑1 (MQ=255)

 ttCAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCg                                                                                                                                                     >  1:259765/1‑150 (MQ=255)

 ttCAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAACAAGGTGCCGGTcc                                                                                                                                                     >  2:2295298/1‑149 (MQ=255)

    aGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCa                                                                                                                                                  <  2:135166/150‑1 (MQ=255)

       ccAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCt                                                                                                                                               <  2:1562498/150‑1 (MQ=255)

       ccAACGCCAGCCAGCCGCTCCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCt                                                                                                                                               <  1:1006888/150‑1 (MQ=255)

        cAACGCCAGCCAGCCGCTGCCCCAGGGCAGCTAGATCACCGTCACCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGGCAGGTCGCCTCCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCtt                                                                                                                                              <  1:6845/150‑1 (MQ=255)

        cAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCtt                                                                                                                                              <  2:927176/150‑1 (MQ=255)

         aaCGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtg                                                                                                                                             <  1:132972/150‑1 (MQ=255)

           cGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTtgct                                                                                                                                           <  2:446735/150‑1 (MQ=255)

           cGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtgct                                                                                                                                           <  1:2268490/150‑1 (MQ=255)

            gccagccaGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCAGGACCTtgcct                                                                                                                                          >  1:1066942/1‑148 (MQ=255)

            gccagccaGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtgctg                                                                                                                                          <  1:118274/150‑1 (MQ=255)

                gccagccGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGACCCCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTgcaag                                                                                                                                      <  1:1979715/150‑1 (MQ=255)

                  cagccGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTgcaaggc                                                                                                                                    >  2:2018130/1‑150 (MQ=255)

                   agccGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTgcaaggca                                                                                                                                   <  1:2312625/150‑1 (MQ=255)

                       gCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGt                                                                                                                               >  2:880169/1‑150 (MQ=255)

                           ccccAGGGCAGCCAGATCACCGTCAGCCAGCCCGAAAGCAATCGCCTGGTCCTCGTGCTTGAATATTGCTGTGCCAGCCAGGTCGCCACCCTCACCCTGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTg                                                                                                                           <  1:892446/150‑1 (MQ=255)

                            cccAGTGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGa                                                                                                                          <  2:269537/150‑1 (MQ=255)

                            cccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGa                                                                                                                          <  1:1861651/150‑1 (MQ=255)

                            cccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGa                                                                                                                          <  1:2285080/150‑1 (MQ=255)

                            cccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCACCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtgc                                                                                                                                            >  1:2149098/1‑132 (MQ=255)

                                gggCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGaccaa                                                                                                                      <  2:1951194/150‑1 (MQ=255)

                                 ggCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCTCCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGaccaac                                                                                                                     <  2:914265/150‑1 (MQ=255)

                                 ggCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGaccaac                                                                                                                     <  2:692850/150‑1 (MQ=255)

                                   cAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGaccaacca                                                                                                                   >  2:2188846/1‑150 (MQ=255)

                                   cAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGaccaacca                                                                                                                   <  2:642238/150‑1 (MQ=255)

                                    aGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑g                                                                                                                 >  2:2035395/1‑150 (MQ=255)

                                     gCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑gc                                                                                                                <  1:880169/150‑1 (MQ=255)

                                      ccAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑gct                                                                                                               >  1:2072359/1‑150 (MQ=255)

                                        aGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑gctgc                                                                                                             <  1:2018130/150‑1 (MQ=255)

                                        aGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑gctgc                                                                                                             >  1:151764/1‑150 (MQ=255)

                                             aCCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGa                                                                                                                          <  2:1440914/133‑1 (MQ=255)

                                               cGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCa                                                                                                      >  2:2116205/1‑150 (MQ=255)

                                                gTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAg                                                                                                     <  1:382123/150‑1 (MQ=255)

                                                 tCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAgg                                                                                                    >  1:1884084/1‑150 (MQ=255)

                                                     ccAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAg                                                                                                     <  2:1351475/145‑1 (MQ=255)

                                                     ccAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAg                                                                                                     >  1:1351475/1‑145 (MQ=255)

                                                      cAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCAccc                                                                                               >  2:252317/1‑150 (MQ=255)

                                                        gACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGaccaac                                                                                                                     <  1:380744/127‑1 (MQ=255)

                                                        gACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGaccaac                                                                                                                     >  2:380744/1‑127 (MQ=255)

                                                          ccGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCAGGCTCACCCAGCTAACGGCGATCAAGGTGCCGGTAGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAAGCACCCGGcc                                                                                           >  1:1741072/1‑150 (MQ=255)

                                                            gAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCaa                                                                                         >  1:1660328/1‑150 (MQ=255)

                                                            gAAAGCAATCGCCTGCCGACCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCaa                                                                                         >  2:598447/1‑150 (MQ=255)

                                                                  aaTCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGcc                                                                                   <  2:2025305/150‑1 (MQ=255)

                                                                   aTCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGa                                                                                                                          >  2:238464/1‑111 (MQ=255)

                                                                   aTCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCa                                                                                                      <  2:2165045/130‑1 (MQ=255)

                                                                   aTCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCa                                                                                                      >  1:2165045/1‑130 (MQ=255)

                                                                         tGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTg                                                                                                                           <  1:1338148/104‑1 (MQ=255)

                                                                         tGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTAcc                                                                            >  2:1224913/1‑150 (MQ=255)

                                                                          ggCGATCATCCTGCAACCGGCCAGCTCCAGCCCGGTCGCCACTCTCACCCAACTGGACACCAGCAAGGTAGCGGGCGACACCTTCCTGTAAGGAAAGATCCTGACCAACCA‑GCTGCTGCCCCAGGCACCCGGCCAACCGGCCGAGTAcca                                                                           >  2:786276/1‑149 (MQ=255)

                                                                           gCGATCATGCGGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCgg                                                                          >  1:2154011/1‑150 (MQ=255)

                                                                                 aTGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTg                                                                    <  2:2072359/150‑1 (MQ=255)

                                                                                 aTGCTGCAACAGGCCAGCGCCAGCCACGTCCACACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCGTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTg                                                                    <  1:1843715/150‑1 (MQ=255)

                                                                                     tgcAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCTAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCa                                                                <  1:2295298/150‑1 (MQ=255)

                                                                                     tgcAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGaccaa                                                                                                                      >  1:708864/1‑97 (MQ=255)

                                                                                     tgcAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGaccaa                                                                                                                      <  2:708864/97‑1 (MQ=255)

                                                                                     tgcAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCa                                                                <  1:1224913/150‑1 (MQ=255)

                                                                                      gcAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAg                                                               <  1:1450138/150‑1 (MQ=255)

                                                                                        aaCAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCC‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGcc                                                             >  1:630541/1‑150 (MQ=255)

                                                                                        aaCAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCg                                                                                              <  2:1020177/117‑1 (MQ=255)

                                                                                        aaCAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCg                                                                                              >  1:1020177/1‑117 (MQ=255)

                                                                                         aCAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGcaaggcaagg                                                                                                                                <  1:1417849/83‑1 (MQ=255)

                                                                                         aCAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGcaaggcaagg                                                                                                                                >  2:1417849/1‑83 (MQ=255)

                                                                                             gccagcgccagcCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTc                                                        >  1:547429/1‑150 (MQ=255)

                                                                                                 gcgccagcCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACa                                                    <  2:547429/150‑1 (MQ=255)

                                                                                                  cgccagcCAGGTCGCCACCCTCCCCCAGCTGGACAACAGCAAGGCGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAg                                                   <  2:1741072/150‑1 (MQ=255)

                                                                                                  cgccagcCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAg                                                   <  2:1797724/150‑1 (MQ=255)

                                                                                                     cagccagGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCaa                                                >  2:1858973/1‑150 (MQ=255)

                                                                                                      agccagGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAAc                                               <  2:1884084/150‑1 (MQ=255)

                                                                                                       gccagGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAAcc                                              >  1:2356457/1‑150 (MQ=255)

                                                                                                         cagGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCggg                                                                         <  2:1254622/121‑1 (MQ=255)

                                                                                                         cagGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCggg                                                                         >  1:1254622/1‑121 (MQ=255)

                                                                                                           ggTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCGTGCCCAGGCACCCGTGCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGc                                          >  2:1173213/1‑150 (MQ=255)

                                                                                                            gTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtgctg                                                                                                                                          >  1:628456/1‑54 (MQ=21)

                                                                                                                ccACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTg                                     >  2:1117900/1‑150 (MQ=255)

                                                                                                                 cACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACATCAACCAGGCCGGTGc                                    <  2:2154011/150‑1 (MQ=255)

                                                                                                                  aCCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGc                                                                                    <  2:747967/101‑1 (MQ=255)

                                                                                                                  aCCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGcc                                   >  1:262241/1‑150 (MQ=255)

                                                                                                                  aCCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGcc                                   >  1:2236569/1‑150 (MQ=255)

                                                                                                                   cccTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCca                                  <  1:1117900/150‑1 (MQ=255)

                                                                                                                   cccTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCca                                  <  2:810803/150‑1 (MQ=255)

                                                                                                                    ccTCACCCAGCTGGACACCAGCAAGGTGCCGGTGGGCACCTTGCTGCAAAGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCcac                                 >  1:1961361/1‑150 (MQ=255)

                                                                                                                    ccTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACTATGGTATAGCAATTGCCAACCTG‑CAATTTTCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGGGCcac                                 >  1:782611/1‑150 (MQ=255)

                                                                                                                       cACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACt                              <  2:151764/150‑1 (MQ=255)

                                                                                                                        aCCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTg                             >  1:1434017/1‑150 (MQ=255)

                                                                                                                         cccAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGa                            <  1:292850/150‑1 (MQ=255)

                                                                                                                          ccAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAg                                                               >  2:1228483/1‑114 (MQ=255)

                                                                                                                             gCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCAt                        <  1:2188846/150‑1 (MQ=255)

                                                                                                                                   ttgcAGCAAGGTGCCGACCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGc                                                                                            >  2:1701397/4‑76 (MQ=37)

                                                                                                                                   ttgcAGCAAGGTGCCGACCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGc                                                                                            <  1:1701397/73‑1 (MQ=37)

                                                                                                                                    acCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgc                 <  2:2236569/150‑1 (MQ=255)

                                                                                                                                      cAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgccc               >  1:824603/1‑150 (MQ=255)

                                                                                                                                       aGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgcccg              >  1:911495/1‑150 (MQ=255)

                                                                                                                                         cAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgcccgcg            >  2:618180/1‑150 (MQ=255)

                                                                                                                                                ccGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTg     <  2:630541/150‑1 (MQ=255)

                                                                                                                                                   gTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGCGCATCGACAGCCCGCGCCCGCTGGTc  >  2:692139/1‑150 (MQ=255)

                                                                                                                                                   gTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAAACAGGCCGGTGCCACATTGAGCGCTGACAGCCCGCCCCCGCTGGTc  >  2:21647/1‑150 (MQ=255)

CTTCAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCA‑GCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTC  >  NC_002947/4917863‑4918159

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: `‑`