BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	727,319	A→T	5.2%	W6R (TGG→AGG) ‡	PP_0621 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	727,319	0	A	T	5.2%	1050.1 / 29.9	308	W6R (TGG→AGG) ‡	PP_0621	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (150/142); new base T (8/8); total (158/150)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 313 total aligned reads displayed.

GTGCGACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGGCATC  >  NC_002947/727170‑727468

gTGCGACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCa                                                                                                                                                       >  1:1828921/1‑150 (MQ=255)

 tGCGACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCaa                                                                                                                                                      <  2:5031805/150‑1 (MQ=255)

    gACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGc                                                                                                                                                   >  1:1025931/1‑150 (MQ=255)

    gACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGc                                                                                                                                                   >  1:4190781/1‑150 (MQ=255)

    gACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGc                                                                                                                                                   <  1:5284268/150‑1 (MQ=255)

         ctggctgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTccc                                                                                                                                              >  1:2148784/1‑150 (MQ=255)

          tggctgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCt                                                                                                                                             <  2:3900369/150‑1 (MQ=255)

           ggctgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCTTGGGCCTCCCtt                                                                                                                                            <  1:4694965/150‑1 (MQ=255)

            gctgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCTTGGGCCTCCCttt                                                                                                                                           <  1:4537204/150‑1 (MQ=255)

              tgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCTTGGGCCTCCCttttt                                                                                                                                         <  2:3699306/150‑3 (MQ=255)

              tgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTttca                                                                                                                                         <  2:6778826/150‑1 (MQ=255)

              tgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTttca                                                                                                                                         <  2:726077/150‑1 (MQ=255)

                  tttGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTttcattca                                                                                                                                     >  2:3064090/1‑150 (MQ=255)

                  tttGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTttcattca                                                                                                                                     >  1:4136194/1‑150 (MQ=255)

                     gTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCg                                                                                                                                  >  1:3850233/1‑150 (MQ=255)

                      tATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAATGGAGGCCCAAGGCCTCCCTTTCATTCACCGc                                                                                                                                 >  2:5234338/1‑150 (MQ=255)

                          tGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGaa                                                                                                                             <  1:2264171/150‑1 (MQ=255)

                           gcgcCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAg                                                                                                                            <  2:220289/150‑1 (MQ=255)

                             gcCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGt                                                                                                                          >  2:4533778/1‑150 (MQ=255)

                               cacaGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCa                                                                                                                        <  2:6226637/150‑1 (MQ=255)

                               cacaGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCa                                                                                                                        <  1:444393/150‑1 (MQ=255)

                                  aGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTc                                                                                                                     <  1:5805750/150‑1 (MQ=255)

                                  aGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTc                                                                                                                     <  2:3860409/150‑1 (MQ=255)

                                  aGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTc                                                                                                                     <  2:5697740/150‑1 (MQ=255)

                                   gCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCg                                                                                                                    >  1:3759516/1‑150 (MQ=255)

                                    ccAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCCCTGAAGGTCAGTCGa                                                                                                                   >  2:6508101/1‑150 (MQ=255)

                                     cAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGaa                                                                                                                  <  2:108441/150‑1 (MQ=255)

                                       gAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACa                                                                                                                <  2:4235090/150‑1 (MQ=255)

                                        aaCGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAg                                                                                                               >  2:1513349/1‑150 (MQ=255)

                                         aCGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGa                                                                                                              <  1:338315/150‑1 (MQ=255)

                                          cgcgGCGGCGATCTGGGCGATCAAGATAAGTAGGCGAACGATGGTGTCCATCTTGTCGATTGAAGAAGCCAACAGTGTAGGACGCGTCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGAc                                                                                                             <  2:2346834/150‑1 (MQ=255)

                                            cggcggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCa                                                                                                           <  2:4674361/150‑1 (MQ=255)

                                             ggcggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCaa                                                                                                          >  2:1290785/1‑150 (MQ=255)

                                              gcggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCaaa                                                                                                         <  2:3119907/150‑1 (MQ=255)

                                               cggcgATCAGGGCGATCCAGAAAAGTAGGCTAACCATTTTGTCCATCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCACTTTCATTCACCGCTGAAGTTCAGTCGAACAGACCAAAg                                                                                                        <  2:5308911/150‑1 (MQ=255)

                                               cggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAg                                                                                                        >  1:5095793/1‑150 (MQ=255)

                                                ggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAgg                                                                                                       >  2:1818826/1‑150 (MQ=255)

                                                   gATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCa                                                                                                    >  2:331254/1‑150 (MQ=255)

                                                   gATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCa                                                                                                    >  1:6398746/1‑150 (MQ=255)

                                                   gATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCa                                                                                                    <  2:854935/150‑1 (MQ=255)

                                                    aTCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCAt                                                                                                   <  1:6873206/150‑1 (MQ=255)

                                                      cAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGt                                                                                                 <  1:2857284/150‑1 (MQ=255)

                                                        gggCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAg                                                                                               >  2:4568502/1‑150 (MQ=255)

                                                         ggCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGc                                                                                              >  2:597484/1‑150 (MQ=255)

                                                          gCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCt                                                                                             >  2:2500607/1‑150 (MQ=255)

                                                          gCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCt                                                                                             <  1:1733078/150‑1 (MQ=255)

                                                                cAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCa                                                                                       >  2:1851550/1‑150 (MQ=255)

                                                                cAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCa                                                                                       <  2:2982848/150‑1 (MQ=255)

                                                                cAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCa                                                                                       <  1:1201087/150‑1 (MQ=255)

                                                                cAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCa                                                                                       <  1:4057193/150‑1 (MQ=255)

                                                                      aGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACg                                                                                 <  2:5397194/150‑1 (MQ=255)

                                                                       ggaGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCCATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACgg                                                                                <  2:4081642/148‑1 (MQ=255)

                                                                             gAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTtc                                                                          <  2:4169426/150‑1 (MQ=255)

                                                                               aCCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTtctc                                                                        <  1:6209852/150‑1 (MQ=255)

                                                                                         ccAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTc                                                              <  1:2707562/150‑1 (MQ=255)

                                                                                          cAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCg                                                             >  1:4960374/1‑150 (MQ=255)

                                                                                           aGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCgg                                                            <  1:2215531/150‑1 (MQ=255)

                                                                                           aGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCgg                                                            >  1:2769358/1‑150 (MQ=255)

                                                                                              ttGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTttc                                                                                                                                          <  2:1134305/69‑1 (MQ=255)

                                                                                                gTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCtt                                                       >  1:1593272/1‑150 (MQ=255)

                                                                                                 tCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTc                                                      >  1:5158816/1‑150 (MQ=255)

                                                                                                        aagaagCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACg                                               >  2:6156618/1‑150 (MQ=255)

                                                                                                        aagaagCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACg                                               >  2:5451217/1‑150 (MQ=255)

                                                                                                        aagaagCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCAAGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACg                                               >  2:2220781/1‑150 (MQ=255)

                                                                                                         agaagCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGa                                              >  1:1989051/1‑150 (MQ=255)

                                                                                                          gaagCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGAt                                             >  2:3160433/1‑150 (MQ=255)

                                                                                                           aagCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGAtt                                            >  1:5513454/1‑150 (MQ=255)

                                                                                                               cAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGt                                        <  1:6766999/150‑1 (MQ=255)

                                                                                                               cAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGACTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGt                                        >  2:1069698/1‑150 (MQ=255)

                                                                                                                aaCAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTc                                       >  2:4814048/1‑150 (MQ=255)

                                                                                                                   aGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGcc                                    <  1:908437/150‑1 (MQ=255)

                                                                                                                     tCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGt                                  >  1:4458870/1‑150 (MQ=255)

                                                                                                                     tCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCGTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCGTCTTTCGGGCCTTCCGGCACGAGTGGGTCGCCGt                                  >  2:6409151/1‑150 (MQ=255)

                                                                                                                       tAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTtc                                >  1:1378775/1‑150 (MQ=255)

                                                                                                                       tAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTtc                                >  2:1711739/1‑150 (MQ=255)

                                                                                                                       tAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTtc                                >  1:6511414/1‑150 (MQ=255)

                                                                                                                           aCGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGt                            <  1:1445134/150‑1 (MQ=255)

                                                                                                                             gcgGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTTGGTCGCCGTTCTCGTCt                          >  1:2509995/1‑150 (MQ=255)

                                                                                                                                gCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTa                       <  2:6387195/150‑1 (MQ=255)

                                                                                                                                gCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTa                       <  2:5791285/150‑1 (MQ=255)

                                                                                                                                gCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGCCTTTa                       <  1:3573543/150‑1 (MQ=255)

                                                                                                                                  cAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAgg                     <  2:4351956/150‑1 (MQ=255)

                                                                                                                                  cAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAgg                     >  1:4390155/1‑150 (MQ=255)

                                                                                                                                  cAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAgg                     <  1:5551119/150‑1 (MQ=255)

                                                                                                                                   aaTAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgc                    >  1:4622515/1‑150 (MQ=255)

                                                                                                                                   aaTAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgc                    >  1:6139030/1‑150 (MQ=255)

                                                                                                                                   aaTAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgc                    >  1:4743475/1‑150 (MQ=255)

                                                                                                                                   aaTAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgc                    >  1:5331725/1‑150 (MQ=255)

                                                                                                                                     caaaaaaGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgcag                  <  2:5606830/149‑1 (MQ=255)

                                                                                                                                      atgaaaGGGAGGCCTTGGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgcagc                 >  2:1418778/4‑150 (MQ=255)

                                                                                                                                        aaaaGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGcagcag               >  2:2905998/1‑150 (MQ=255)

                                                                                                                                          aaGGGAGGCCTTGGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGct             >  1:3246586/1‑150 (MQ=255)

                                                                                                                                           aGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGctc            >  1:648898/1‑150 (MQ=255)

                                                                                                                                               aGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTcgc        <  1:6600880/150‑1 (MQ=255)

                                                                                                                                                  cccAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGGc     >  2:6678256/1‑150 (MQ=255)

                                                                                                                                                  cccAAGGCCTCCCTGTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCGTGTAGCTGAACCAGGAACGGTCTTTCTCGGCGTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGGCTTTAGGCAGCAGCTCTCGCGGc     >  2:2674935/1‑150 (MQ=255)

                                                                                                                                                   ccAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGGCa    >  2:3534028/1‑150 (MQ=255)

                                                                                                                                                   ccAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGGCa    >  1:4261665/1‑150 (MQ=255)

                                                                                                                                                   ccAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGGCa    >  2:2278331/1‑150 (MQ=255)

                                                                                                                                                    cAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGACCCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGGCAt   >  2:5171384/1‑150 (MQ=255)

GTGCGACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGGCATC  >  NC_002947/727170‑727468

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