BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,223,232	T→A	5.1%	S95T (TCG→ACG)	PP_4598 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,223,232	0	T	A	5.1%	875.7 / 22.2	255	S95T (TCG→ACG)	PP_4598	membrane protein
Reads supporting (aligned to +/- strand): ref base T (134/108); new base A (8/5); total (142/113)
Fisher's exact test for biased strand distribution p-value = 7.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 273 total aligned reads displayed.

ACCCGGCGTATTTCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCGGGCCGAT  >  NC_002947/5223083‑5223381

aCCCGGCGTATTTCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCt                                                                                                                                                       >  1:3224466/1‑150 (MQ=255)

     gCGTATTTCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCtcgttc                                                                                                                                                  <  2:1842227/150‑1 (MQ=255)

      cGTATTTCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCtcgttcg                                                                                                                                                 >  1:2774099/1‑150 (MQ=255)

         aTTTCTACCCGGTCAGCGCCACGCCGTTCTCCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGc                                                                                                                                              <  1:3382543/150‑1 (MQ=255)

         aTTTCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGc                                                                                                                                              >  2:4648926/1‑150 (MQ=255)

          tttCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCg                                                                                                                                             >  2:5320741/1‑150 (MQ=255)

           ttCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTGGGTGTCGATCTTCATGAGCTTTCGCAACATCGCGTGCTATGACCGTTGGTGGGTGGGGCGCAAGGCCTTGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCgg                                                                                                                                            <  1:1168072/150‑1 (MQ=255)

           ttCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCgg                                                                                                                                            <  2:5005720/150‑1 (MQ=255)

            tCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCggg                                                                                                                                           >  2:3446723/1‑150 (MQ=255)

               aCCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGaaa                                                                                                                                        <  2:808948/150‑1 (MQ=255)

               aCCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGaaa                                                                                                                                        >  1:3531978/1‑150 (MQ=255)

               aCCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGaaa                                                                                                                                        >  2:515371/1‑150 (MQ=255)

               aCCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCACGAACGAGCgggttt                                                                                                                                        <  2:1489674/150‑4 (MQ=255)

                    gTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCgtg                                                                                                                                   <  2:4565897/150‑1 (MQ=255)

                     tCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCgtgg                                                                                                                                  <  2:3198906/150‑1 (MQ=255)

                      cAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCgtggt                                                                                                                                 <  1:2951367/150‑1 (MQ=255)

                      cAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCgtggt                                                                                                                                 <  2:1025032/150‑1 (MQ=255)

                            cACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTAc                                                                                                                           <  2:3998504/150‑1 (MQ=255)

                               gCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACcga                                                                                                                        <  1:244117/150‑1 (MQ=255)

                                    tCACCTTGCTGGGCCTGTCGCTGTCAATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGaaa                                                                                                                   <  1:1847184/150‑1 (MQ=255)

                                      aCCTTGCTGGGCCTGTCTCTGTCGATCTTCATGAGCTTTCTCAACAACGCTTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAaca                                                                                                                 <  1:4361492/150‑1 (MQ=255)

                                      aCCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAaca                                                                                                                 <  1:1317803/150‑1 (MQ=255)

                                      aCCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAaca                                                                                                                 <  1:4120140/150‑1 (MQ=255)

                                        cTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACt                                                                                                               >  1:4282366/1‑150 (MQ=255)

                                              gggCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCAc                                                                                                         >  2:3895158/1‑150 (MQ=255)

                                                   tgtcgctgtcgATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAAc                                                                                                    >  1:4445996/1‑150 (MQ=255)

                                                   tgtcgctgtcgATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAAc                                                                                                    >  2:2655970/1‑150 (MQ=255)

                                                       gctgtcgATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAActgct                                                                                                >  1:3834260/1‑150 (MQ=255)

                                                       gctgtcgATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAActgct                                                                                                >  1:1390112/1‑150 (MQ=255)

                                                       gctgtcgATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAActgct                                                                                                >  2:4833863/1‑150 (MQ=255)

                                                       gctgtcgATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAActgct                                                                                                >  1:3423360/1‑150 (MQ=255)

                                                         tgtcgATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACtgctgc                                                                                              >  1:1786056/1‑150 (MQ=255)

                                                         tgtcgATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACtgctgc                                                                                              <  2:1153838/150‑1 (MQ=255)

                                                          gtcgATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTgcg                                                                                             >  1:3965020/1‑150 (MQ=255)

                                                          gtcgATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTgcg                                                                                             >  1:4687259/1‑150 (MQ=255)

                                                           tcgATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTgcgc                                                                                            <  1:2292166/150‑1 (MQ=255)

                                                           tcgATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTgcgc                                                                                            >  2:3861114/1‑150 (MQ=255)

                                                              aTCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAAc                                                                                         <  2:2250452/150‑1 (MQ=255)

                                                                  tCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCtgt                                                                                     >  2:4490046/1‑150 (MQ=255)

                                                                  tCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCtgt                                                                                     <  1:1079898/150‑1 (MQ=255)

                                                                    aTGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCtgtgt                                                                                   <  2:4926572/150‑1 (MQ=255)

                                                                    aTGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCtgtgt                                                                                   <  2:1253928/150‑1 (MQ=255)

                                                                    aTGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCtgtgt                                                                                   <  2:106166/150‑1 (MQ=255)

                                                                     tGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCtgtgtg                                                                                  <  2:796147/150‑1 (MQ=255)

                                                                     tGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCtgtgtg                                                                                  >  1:1985335/1‑150 (MQ=255)

                                                                      gAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTgg                                                                                 >  2:3790461/1‑150 (MQ=255)

                                                                        gCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGGCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGTGAAAGCGTGGTGATTACCGACTAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGtt                                                                               >  2:1766679/1‑150 (MQ=255)

                                                                          tttCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTc                                                                             <  1:309382/150‑1 (MQ=255)

                                                                           ttCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCACGAACGAGCgg                                                                                                                                            >  1:477916/1‑86 (MQ=255)

                                                                            tCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGc                                                                           <  1:460583/150‑1 (MQ=255)

                                                                             cGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGcc                                                                          <  1:4044622/150‑1 (MQ=255)

                                                                              gCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGccc                                                                         >  2:4237966/1‑150 (MQ=255)

                                                                              gCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGccc                                                                         >  2:434310/1‑150 (MQ=255)

                                                                               caacaaCGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCa                                                                        <  1:5597458/150‑1 (MQ=255)

                                                                                    acGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCAc                                                                   >  1:2635775/1‑150 (MQ=255)

                                                                                         tGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAAt                                                              <  2:2419455/150‑1 (MQ=255)

                                                                                           cTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGc                                                            >  2:3046858/1‑150 (MQ=255)

                                                                                             aTGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGccc                                                          <  1:944895/150‑1 (MQ=255)

                                                                                             aTGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGccc                                                          <  1:5024683/150‑1 (MQ=255)

                                                                                                aCCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGAACTCAATGACCGGc                                                       >  2:3933874/1‑150 (MQ=255)

                                                                                                 ccGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCt                                                      <  2:1849188/150‑1 (MQ=255)

                                                                                                 ccGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCt                                                      >  2:731885/1‑150 (MQ=255)

                                                                                                    ttGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCg                                                   >  1:212573/1‑150 (MQ=255)

                                                                                                       gtggGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACACGCTCGTTCGTGAGGGAAAGCGTGGTGATTACCGACGAAACACAGCGCACCGAAATGATGCGCAACCTGTGTGGTTTCGCACAAGCACTCAATACCCGGCTGCGaaa                                                >  1:571610/1‑150 (MQ=255)

                                                                                                        tggGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGTGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAAc                                               <  1:3460856/150‑1 (MQ=255)

                                                                                                          ggAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGa                                             <  2:2590299/150‑1 (MQ=255)

                                                                                                            aGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAgg                                           >  1:4958495/1‑150 (MQ=255)

                                                                                                            aGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAgg                                           >  2:2893562/1‑150 (MQ=255)

                                                                                                            aGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAgg                                           >  1:5288800/1‑150 (MQ=255)

                                                                                                              gggCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGAc                                         >  1:5366860/1‑150 (MQ=255)

                                                                                                                gcgcAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCAATGCACTCAATGCCCGGCTGCGTAACGAGGACGa                                       >  1:3928859/1‑150 (MQ=255)

                                                                                                                  gcAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGacc                                     >  1:5514176/1‑148 (MQ=255)

                                                                                                                  gcAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAAc                                     >  2:1622672/1‑150 (MQ=255)

                                                                                                                   cAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACg                                    >  2:4751342/1‑150 (MQ=255)

                                                                                                                   cAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACg                                    >  1:689682/1‑150 (MQ=255)

                                                                                                                   cAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACg                                    <  2:3609993/150‑1 (MQ=255)

                                                                                                                     aGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGa                                             <  1:1152409/139‑1 (MQ=255)

                                                                                                                     aGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGa                                             >  2:1152409/1‑139 (MQ=255)

                                                                                                                          tGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCt                             <  2:4019930/150‑1 (MQ=255)

                                                                                                                                aaaCTGATTCTCGAAACCCGCTCGTTCGTGCGGGAACGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCGGGGTGGTGTCGCCCTTGCACTCGATGCCCGGGTGCGTACCGAGGGCGACCGGGCAGCTGCCCgg                       >  2:120457/1‑150 (MQ=255)

                                                                                                                                aaaCTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAgg                       >  2:3124398/1‑150 (MQ=255)

                                                                                                                                aaaCTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAgg                       >  1:1253944/1‑150 (MQ=255)

                                                                                                                                  aCTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGcc                     >  2:889805/1‑150 (MQ=255)

                                                                                                                                  aCTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGcc                     >  2:4167502/1‑150 (MQ=255)

                                                                                                                                  aCTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGcc                     >  2:2096581/1‑150 (MQ=255)

                                                                                                                                      aTTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGtgg                 <  2:2774099/150‑1 (MQ=255)

                                                                                                                                       ttATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGtggc                >  1:2252551/1‑150 (MQ=255)

                                                                                                                                        tATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGtggct               >  2:1035384/1‑150 (MQ=255)

                                                                                                                                        tATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGtggct               <  1:1516542/150‑1 (MQ=255)

                                                                                                                                          tCGAAACCCGCTCGTTCGTGCGGGAACGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGtggctgg             >  2:5074360/1‑150 (MQ=255)

                                                                                                                                          tCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGtggctgg             >  1:2656998/1‑150 (MQ=255)

                                                                                                                                             tttcccGCACGAACGAGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGt                   <  2:348302/138‑1 (MQ=255)

                                                                                                                                             aaaCCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCACGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGccc          >  2:1531720/1‑150 (MQ=255)

                                                                                                                                              ttcccGCACGAACGAGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAg                                            <  1:3456018/113‑1 (MQ=255)

                                                                                                                                              aaCCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCg         <  1:4000242/150‑1 (MQ=255)

                                                                                                                                               aCCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCgg        <  2:2325501/150‑1 (MQ=255)

                                                                                                                                                cccGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCggg       >  2:3698766/1‑150 (MQ=255)

                                                                                                                                                 ccGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCGGGc      >  1:4159394/1‑150 (MQ=255)

                                                                                                                                                  cGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCGGGcc     <  1:4017478/150‑1 (MQ=255)

                                                                                                                                                    cTCGTTCGTGCTGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATACCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCGGGCCGa   <  2:571610/150‑1 (MQ=255)

ACCCGGCGTATTTCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCGGGCCGAT  >  NC_002947/5223083‑5223381

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