BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,389,952	C→G	8.9%	A23G (GCC→GGC)	PP_1212 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,389,952	0	C	G	8.9%	306.7 / 18.0	103	A23G (GCC→GGC)	PP_1212	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (44/48); new base G (4/5); total (49/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

Only 100 of 107 total aligned reads displayed.

GCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGC  >  NC_002947/1389816‑1390099

gctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATc                                                                                                                                        <  2:2102286/150‑1 (MQ=255)

gctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATc                                                                                                                                        <  2:1965586/150‑1 (MQ=255)

gctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATc                                                                                                                                        <  2:1697409/150‑1 (MQ=255)

gctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATc                                                                                                                                        >  2:1471588/1‑150 (MQ=255)

  tAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATctt                                                                                                                                      <  2:1792970/150‑3 (MQ=255)

  tAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGa                                                                                                                                      >  1:1431904/1‑150 (MQ=255)

  tAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGa                                                                                                                                      >  1:1776203/1‑150 (MQ=255)

   aaCCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGAt                                                                                                                                     >  1:87119/1‑150 (MQ=255)

     ccATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCTGCGCCGCTGATCGATTg                                                                                                                                   >  1:946285/1‑150 (MQ=255)

       atatatAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGcc                                                                                                                                 >  1:721507/1‑150 (MQ=255)

         atataATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGTAAGTGCTGCAGAAGATCAGCGTCGCGCCGCTGATCGATTGCCCg                                                                                                                               <  2:788606/150‑1 (MQ=255)

         atataATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCg                                                                                                                               <  1:1895082/150‑1 (MQ=255)

           ataGTCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGc                                                                                                                             <  2:1027783/150‑1 (MQ=255)

            taATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGCCACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGcc                                                                                                                            >  2:2141383/1‑150 (MQ=255)

             aaTCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCt                                                                                                                           <  2:946285/150‑1 (MQ=255)

               tCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGAtt                                                                                                                                    <  2:386884/139‑1 (MQ=255)

               tCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGAtt                                                                                                                                    >  1:386884/1‑139 (MQ=255)

                 ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATc                                                                                                                                        >  2:142330/1‑133 (MQ=255)

                 ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATc                                                                                                                                        <  1:142330/133‑1 (MQ=255)

                 ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGAGCTGCAGAAGATCAGCGCCGCGCCGGTGATCGATTGCACGAccgtgcc                                                                                                                       >  2:1116951/1‑145 (MQ=255)

                 ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGCCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCgccgcgccgc                                                                                                                                             <  2:2298498/128‑1 (MQ=255)

                 ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGCCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCgccgcgccgc                                                                                                                                             >  1:2298498/1‑128 (MQ=255)

                  cGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGAt                                                                                                                                     >  2:1072451/1‑135 (MQ=255)

                  cGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGAt                                                                                                                                     <  1:1072451/135‑1 (MQ=255)

                  cGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAg                                                                                                                      >  2:875178/1‑150 (MQ=255)

                   gCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGACCGATTGCCCGGCCTGTCAgg                                                                                                                     >  2:1252656/1‑150 (MQ=255)

                          gctgcCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGc                                                                                                              >  1:1773051/1‑150 (MQ=255)

                            tgcCCACCGCCAGCAAGCGTTATCGTTGATGAAAGACATCCATGCCCCTTTATGACTAGCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCAc                                                                                                            >  2:1234054/1‑150 (MQ=255)

                                ccccGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGaa                                                                                                        >  1:1015213/1‑150 (MQ=255)

                                ccccGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGaa                                                                                                        >  2:1988513/1‑150 (MQ=255)

                                   cGGCAGCCAGCGTTATGGTTGACGAGAGGCACCCATGCCCCTTTATGACTATCAATGTGCGTCCGGCCAGCAGCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGaaaaa                                                                                                     <  1:1852657/150‑1 (MQ=255)

                                   cGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGaaaaa                                                                                                     >  1:1241410/1‑150 (MQ=255)

                                     gCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAgc                                                                                                   <  1:2251957/150‑1 (MQ=255)

                                        gccagcATTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAgctgc                                                                                                <  2:375050/150‑1 (MQ=255)

                                         ccagcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGcc                                                                                                                            <  1:1671000/121‑1 (MQ=255)

                                         ccagcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGcc                                                                                                                            >  2:1671000/1‑121 (MQ=255)

                                           agcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGt                                                                                             >  1:2325361/1‑150 (MQ=255)

                                           agcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGt                                                                                             >  1:886928/1‑150 (MQ=255)

                                               ttATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGt                                                                                         <  1:1116951/150‑1 (MQ=255)

                                                    gTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAg                                                                                    <  2:2325361/150‑1 (MQ=255)

                                                     ttGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAgg                                                                                   >  1:708983/1‑150 (MQ=255)

                                                      tGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGc                                                                                  <  1:1344293/150‑1 (MQ=255)

                                                          gagagaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTc                                                                              <  1:1471588/150‑1 (MQ=255)

                                                               acacCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTg                                                                         >  2:888327/1‑150 (MQ=255)

                                                                  cccATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGc                                                                      <  2:87119/150‑1 (MQ=255)

                                                                   ccATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCg                                                                     >  2:852412/1‑150 (MQ=255)

                                                                   ccATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCg                                                                     >  2:2163335/1‑150 (MQ=255)

                                                                    cATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCTGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCgg                                                                    <  1:1234054/150‑1 (MQ=255)

                                                                         cccTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGa                                                                                                                                          <  1:2199807/75‑1 (MQ=255)

                                                                         cccTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGa                                                                                                                                          >  2:2199807/1‑75 (MQ=255)

                                                                         cccTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGGGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACg                                                               >  1:1630339/1‑150 (MQ=255)

                                                                           cTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGt                                                             <  1:1988513/150‑1 (MQ=255)

                                                                           cTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGt                                                             <  2:1015213/150‑1 (MQ=255)

                                                                           cTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGt                                                             <  2:708983/150‑1 (MQ=255)

                                                                            tttATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGtt                                                            <  1:852412/150‑1 (MQ=255)

                                                                             ttATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTg                                                           <  2:36669/150‑1 (MQ=255)

                                                                                tGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTa                                                        <  2:1431904/150‑1 (MQ=255)

                                                                                 gACTATCAATGTGCGTCCTGCGAGTACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCAGGCCTGTCAGGCGCCGGCAATGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTAc                                                       <  2:351317/150‑1 (MQ=255)

                                                                                    tATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGaa                                                    >  2:2257911/1‑150 (MQ=255)

                                                                                     aTCAATGTGCGTCCTGCGAGCACCGCATGGCAGAGCTGCAGAACAACGGCACAGCGGCGCTTGTCGATTGTGTGGTCTGCTAGGAGGCGGCAGTGAAATCACCGCTTTCTTGTCCGGTCTTTCGCCTGAGCGGTAACGGTTGGTACGaac                                                   <  1:1252656/150‑2 (MQ=255)

                                                                                             gCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTc                                           <  2:1773051/150‑1 (MQ=255)

                                                                                             gCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTc                                           <  1:888327/150‑1 (MQ=255)

                                                                                               gTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCaa                                         >  1:801291/1‑150 (MQ=255)

                                                                                               gTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCCACTTCac                                         <  2:1241410/150‑2 (MQ=255)

                                                                                                  cTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGTCAACCAACGGTATGTGTTCGCCCGGCCTGTCAGGCTCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGAc                                      >  1:420750/1‑150 (MQ=255)

                                                                                                  cTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCa                                                                                     <  1:582538/103‑1 (MQ=255)

                                                                                                  cTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCa                                                                                     >  2:582538/1‑103 (MQ=255)

                                                                                                   tGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAgg                                                                                   >  2:1258925/1‑104 (MQ=255)

                                                                                                   tGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAgg                                                                                   <  1:1258925/104‑1 (MQ=255)

                                                                                                   tGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGAcc                                     >  2:1910877/1‑150 (MQ=255)

                                                                                                    gCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTg                                                           >  2:2159861/1‑127 (MQ=255)

                                                                                                    gCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTg                                                           <  1:2159861/127‑1 (MQ=255)

                                                                                                      gAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGGAGGCAACCAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCggg                                  <  2:1891109/150‑1 (MQ=255)

                                                                                                       atcACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCgggg                                 <  1:478401/148‑1 (MQ=255)

                                                                                                       aGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGAACTGAAAAAGCTGCTGTCGGTGCCAGGCCTTCGCGTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCgggg                                 <  2:721507/150‑1 (MQ=255)

                                                                                                           ccGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTAc                                                       >  1:2392126/1‑124 (MQ=255)

                                                                                                           ccGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTAc                                                       <  2:2392126/124‑1 (MQ=255)

                                                                                                            cGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaaa                            >  1:503246/1‑150 (MQ=255)

                                                                                                              cATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaaaaa                          <  2:886928/150‑1 (MQ=255)

                                                                                                                  gAAGTGCTGCATAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAAt                      <  1:1910877/150‑1 (MQ=255)

                                                                                                                    aGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCt                    <  2:801291/150‑1 (MQ=255)

                                                                                                                        ctgcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCTGCACTGAAAAAGCTGCTGTAGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTTCGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCa                <  1:2297526/150‑1 (MQ=255)

                                                                                                                          gcAGAAGATCAGCGCCGCGCCGCTGATGGATTGCCCGGCCTGTCAGGCGCCCGCAATGTACACCCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTACCGGTTGGTACGAAAACGGCTTTAAGACCGGGGCGAAAAAGAATATGGCAgg              >  2:19473/1‑150 (MQ=255)

                                                                                                                            agaagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggcg            >  1:248136/1‑150 (MQ=255)

                                                                                                                             gaagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGt                                                             <  2:2016076/100‑1 (MQ=255)

                                                                                                                             gaagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGt                                                             >  1:2016076/1‑100 (MQ=255)

                                                                                                                              tcgaTCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggcggc          >  1:536406/3‑150 (MQ=255)

                                                                                                                              aagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGc                                                                            <  1:529197/84‑1 (MQ=255)

                                                                                                                              aagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGc                                                                            >  2:529197/1‑84 (MQ=255)

                                                                                                                              aagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggcggc          >  2:372332/1‑150 (MQ=255)

                                                                                                                              aagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggcggc          <  1:2141383/150‑1 (MQ=255)

                                                                                                                                gaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGaaa                                                   >  1:838803/1‑107 (MQ=255)

                                                                                                                                gaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGaaa                                                   <  2:838803/107‑1 (MQ=255)

                                                                                                                                 aTCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAAcc                                                 <  1:519961/108‑1 (MQ=255)

                                                                                                                                 aTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGAc       <  2:248136/150‑1 (MQ=255)

                                                                                                                                  tCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAAt                      >  2:1845698/1‑134 (MQ=255)

                                                                                                                                  tCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACa      >  2:2119344/1‑150 (MQ=255)

                                                                                                                                    aGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAg    >  2:1817672/1‑150 (MQ=255)

                                                                                                                                     gcgccgcgccGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAgg   >  2:1378199/1‑150 (MQ=255)

                                                                                                                                      cgccgcgccgCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGc  >  2:927815/1‑150 (MQ=255)

GCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGC  >  NC_002947/1389816‑1390099

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