BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,744,747	T→C	6.3%	S27G (AGC→GGC)	PP_4198 ←	OmpA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,744,747	0	T	C	6.3%	279.5 / 6.3	95	S27G (AGC→GGC)	PP_4198	OmpA family protein
Reads supporting (aligned to +/- strand): ref base T (48/41); new base C (4/2); total (52/43)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 103 total aligned reads displayed.

GAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTGCAGCACCTGCAGCTGCTCTTCGCTC  >  NC_002947/4744599‑4744894

gAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCActgct                                                                                                                                                     >  1:408640/1‑149 (MQ=255)

 aaCCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCActgctg                                                                                                                                                    <  4:179313/149‑1 (MQ=255)

 aaCCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCActgctg                                                                                                                                                    <  3:417257/149‑1 (MQ=255)

   ccGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTgcg                                                                                                                                                  <  3:198854/149‑1 (MQ=255)

     ggTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTACTCCAGCGCGGCACTGCTGCGCt                                                                                                                                                <  7:273515/149‑1 (MQ=255)

      gTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCt                                                                                                                                                >  6:445319/1‑148 (MQ=255)

      gTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCt                                                                                                                                                >  7:152009/1‑148 (MQ=255)

       tGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGa                                                                                                                                              >  3:419963/1‑149 (MQ=255)

         ccGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCCTCTTTGACCCCCTCTGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGAc                                                                                                                                             <  1:308510/149‑1 (MQ=255)

          cGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCACAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTg                                                                                                                                           <  6:326524/149‑1 (MQ=255)

           gATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCAGCAGTGCCGCGCTgg                                                                                                                                              >  4:100542/1‑144 (MQ=255)

           gATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCAGCAGTGCCGCGCTgg                                                                                                                                              <  3:100542/145‑2 (MQ=255)

             tGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGgcg                                                                                                                                        >  6:177584/1‑149 (MQ=255)

              gTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGgcg                                                                                                                                        >  3:316564/1‑148 (MQ=255)

                 aCTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTTCTGCGCTGACTGGCGCAAc                                                                                                                                    >  4:287843/1‑149 (MQ=255)

                  cTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAAcc                                                                                                                                   >  6:9590/1‑149 (MQ=255)

                   tGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCt                                                                                                                                  >  8:206102/1‑149 (MQ=255)

                   tGGACAGGCGCTCGCCGCGAGCCAGCGACTCACCCGCGCGAATCACGTCACGCCGCGCGCTGCGCAGCACATCGGAATCATCTTTGCCCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTTCGCTTACTGGCGCAACCt                                                                                                                                  >  3:63932/1‑149 (MQ=255)

                    ggACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCtt                                                                                                                                 >  5:119355/1‑149 (MQ=255)

                     gACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCt                                                                                                                                  <  3:386705/147‑1 (MQ=255)

                     gACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTg                                                                                                                                <  6:338886/149‑1 (MQ=255)

                         ggCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCActgct                                                                                                                                                     >  4:359254/1‑124 (MQ=255)

                          gcgcTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAAc                                                                                                                           <  8:152009/149‑1 (MQ=255)

                           cgcTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCt                                                                                                                                                <  6:242544/127‑1 (MQ=255)

                           cgcTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCt                                                                                                                                                >  5:242544/1‑127 (MQ=255)

                           cgcTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCACCCTTGCAAAcc                                                                                                                          >  6:132485/1‑149 (MQ=255)

                            gcTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCt                                                                                                                         >  4:216234/1‑149 (MQ=255)

                            gcTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCt                                                                                                                         <  5:9590/149‑1 (MQ=255)

                                 gcgcGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGgc                                                                                                                                         <  4:373443/128‑1 (MQ=255)

                                   gcgAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCa                                                                                                                  <  4:419963/149‑1 (MQ=255)

                                       gCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAgcgc                                                                                                              >  5:276907/1‑149 (MQ=255)

                                       gCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAgcgc                                                                                                              <  2:54941/149‑1 (MQ=255)

                                        ccAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAgcgc                                                                                                              >  3:409850/1‑148 (MQ=255)

                                          aGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTg                                                                                                                                >  8:348704/1‑128 (MQ=255)

                                          aGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTg                                                                                                                                <  7:348704/128‑1 (MQ=255)

                                           tcGACTCACCCGCCAGAATCACGTCACGCGGCGCGCTGCGCATCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGc                                                                                                          <  8:308142/148‑1 (MQ=255)

                                           gCGACTCACCCGCCCGAATCACGTCAAGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAg                                                                                                           >  1:292382/1‑148 (MQ=255)

                                                tCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGc                                                                                                     >  2:40998/1‑149 (MQ=255)

                                                  aCCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGc                                                                                                     <  4:409850/147‑1 (MQ=255)

                                                  aCCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGc                                                                                                     <  5:424853/147‑1 (MQ=255)

                                                  aCCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGc                                                                                                     <  2:408640/147‑1 (MQ=255)

                                                   cccgcccgAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGTATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAg                                                                                                  <  7:85699/149‑1 (MQ=255)

                                                   cccgcccgAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAg                                                                                                  >  3:21131/1‑149 (MQ=255)

                                                    ccgcccgAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGt                                                                                                 >  7:155038/1‑149 (MQ=255)

                                                     cgcccgAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTg                                                                                                >  3:307057/1‑149 (MQ=255)

                                                      gcccgAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGc                                                                                               <  4:307057/149‑1 (MQ=255)

                                                      gcccgAATCACCTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGCCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGc                                                                                               <  7:229583/149‑1 (MQ=255)

                                                          gAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGc                                                                                            >  2:74193/1‑148 (MQ=255)

                                                             tcacgtcacgCGGCGCGCTGCGCAGCACATCGGAATCATCTTAGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCa                                                                                          >  7:128795/1‑147 (MQ=255)

                                                                cgtcacgCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGc                                                                                     <  4:316564/149‑1 (MQ=255)

                                                                   cacgCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTa                                                                                  >  4:332161/1‑149 (MQ=255)

                                                                    acgCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTa                                                                                  <  1:296386/148‑1 (MQ=255)

                                                                      gcggcgCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGc                                                                                                     >  4:397776/1‑127 (MQ=255)

                                                                      gcggcgCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGc                                                                                                     <  3:397776/127‑1 (MQ=255)

                                                                       cggcgCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTacgac                                                                              <  3:216234/149‑1 (MQ=255)

                                                                             gcTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCa                                                                        >  5:422626/1‑149 (MQ=255)

                                                                              cTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAATGCAACGGTTGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCAt                                                                       <  1:274609/149‑1 (MQ=255)

                                                                                 cgcAGGACATGGGAATCATCTTTGCCCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGc                                                                    <  3:249571/149‑1 (MQ=255)

                                                                                  gcagcaCATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCt                                                                   <  1:74193/149‑1 (MQ=255)

                                                                                   cagcaCATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTg                                                                  <  4:368697/149‑1 (MQ=255)

                                                                                       acaTCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCAcc                                                              >  5:70237/1‑149 (MQ=255)

                                                                                       acaTCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCAc                                                               <  3:254981/148‑1 (MQ=255)

                                                                                         aTCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTc                                                            >  8:233098/1‑149 (MQ=255)

                                                                                            ggAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCa                                                          <  6:74510/148‑1 (MQ=255)

                                                                                            ggAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCa                                                          <  6:119355/148‑1 (MQ=255)

                                                                                            ggAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCa                                                          <  2:395303/148‑1 (MQ=255)

                                                                                             gAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGt                                                        >  7:330040/1‑149 (MQ=255)

                                                                                                  atcTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGttgctt                                                   >  2:197743/1‑149 (MQ=255)

                                                                                                  atcTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGttgctt                                                   >  1:196216/1‑149 (MQ=255)

                                                                                                       tGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGttgcttgc                                                 <  5:445319/146‑1 (MQ=255)

                                                                                                       tGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGttgcttgc                                                 <  4:131293/146‑1 (MQ=255)

                                                                                                       tGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGttgcttgc                                                 <  1:112684/146‑1 (MQ=255)

                                                                                                         aCCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAg                                            >  1:262398/1‑149 (MQ=255)

                                                                                                          cccccT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAgc                                           <  5:177584/149‑1 (MQ=255)

                                                                                                           ccccT‑GGAAGGCAACGGTGGCCTTATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAgcg                                          >  4:285768/1‑149 (MQ=255)

                                                                                                           ccccT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAgcg                                          >  8:329668/1‑149 (MQ=255)

                                                                                                             ccT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCt                                        <  8:128795/149‑1 (MQ=255)

                                                                                                             ccT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCt                                        >  1:33941/1‑149 (MQ=255)

                                                                                                              cT‑GGACGGCAACGGGGGCCTCATGCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCTTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCtt                                       <  4:63932/149‑1 (MQ=255)

                                                                                                              cT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCtt                                       >  1:120966/1‑149 (MQ=255)

                                                                                                              cT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCtt                                       >  3:358534/1‑149 (MQ=255)

                                                                                                                      gCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCcgc                                <  7:206102/149‑1 (MQ=255)

                                                                                                                       cAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCcgc                                >  2:433155/1‑148 (MQ=255)

                                                                                                                        aaCGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGa                              >  3:275328/1‑149 (MQ=255)

                                                                                                                                ccTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTgcagca                      >  3:145945/1‑149 (MQ=255)

                                                                                                                                ccTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTgcagca                      >  3:155388/1‑149 (MQ=255)

                                                                                                                                  tCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTGCAGCAcc                    >  3:29260/1‑149 (MQ=255)

                                                                                                                                   cagtcAGCGCAGCAGTGCCGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTGCAGCACCt                   >  6:225309/5‑149 (MQ=255)

                                                                                                                                     tCCGGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTGCAGCACCTg                  >  7:22622/1‑148 (MQ=255)

                                                                                                                                     tCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTGCAGCACCTg                  >  1:195185/1‑148 (MQ=255)

                                                                                                                                     gtcAGCGCAGCAGTGCCGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTGCAGCACCTg                  >  4:155374/3‑148 (MQ=255)

                                                                                                                                        aGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTGCAGCACCTGCAgc              >  8:408119/1‑149 (MQ=255)

                                                                                                                                        aGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTGCAGCACCTGCAgc              >  5:300829/1‑149 (MQ=255)

                                                                                                                                          cgcgGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCATGCGCAAGCGCTTGACCCGCGACTGCAGCACCTGCAgctg            <  6:422626/149‑1 (MQ=255)

                                                                                                                                          cgcaGCAGTGCCGCGCTGGCTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTGCAGCACCTGCAgctg            >  6:26265/1‑149 (MQ=255)

                                                                                                                                           gcagCAGTGCCGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTGCAGCACCTGCAgctgc           <  5:225309/146‑1 (MQ=255)

                                                                                                                                             ggCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTGCAGGTGCTGCAgtcgcgg         <  6:70237/149‑7 (MQ=255)

                                                                                                                                                 ctgctgCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTGCAGCACCTGCAGCTgctcttcg     <  7:233098/149‑1 (MQ=255)

                                                                                                                                                 ctgctgCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTGCAGCACCTGCAGCTgctcttcg     <  2:196216/149‑1 (MQ=255)

                                                                                                                                                    ctgcGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGATGCTTGCGCAAGCGCTTGACCCGCGACTGCAGCACCTGCAGCTgctcttcgctc  >  1:242344/1‑149 (MQ=255)

GAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTGCAGCACCTGCAGCTGCTCTTCGCTC  >  NC_002947/4744599‑4744894

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