BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,744,743	C→G	6.3%	S28T (AGT→ACT)	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,743	0	C	G	6.3%	312.0 / 9.3	95	S28T (AGT→ACT)	PP_4198	OmpA family protein
Reads supporting (aligned to +/- strand): ref base C (48/41); new base G (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 = 4.42e-01

Only 100 of 104 total aligned reads displayed.

ATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTGCAGCACCTGCAGCTGCTCTTCG  >  NC_002947/4744595‑4744891

aTCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCAc                                                                                                                                                      >  6:424853/1‑149 (MQ=255)

  cGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCActg                                                                                                                                                    >  3:131293/1‑149 (MQ=255)

   ggAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCActg                                                                                                                                                    <  5:102054/148‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

               gATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCAGCAGTGCCGCGCTgg                                                                                                                                           >  4:100542/1‑144 (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)

                              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                                                                                                                                      >  3:373443/1‑128 (MQ=255)

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

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

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

                                           gCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAgcgc                                                                                                           >  5:276907/1‑149 (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                                                                                                  <  2:408640/147‑1 (MQ=255)

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

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

                                                       cccgcccgAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGTATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAg                                                                                               <  7:85699/149‑1 (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)

                                                                                         gcaCATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCa                                                             >  5:74510/1‑149 (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:119355/148‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                                                           tGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGttgcttgc                                              <  4:131293/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                                     >  1:33941/1‑149 (MQ=255)

                                                                                                                 ccT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCt                                     <  8:128795/149‑1 (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           >  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  <  2:196216/149‑1 (MQ=255)

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

ATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCT‑GGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGACTGGCGCAACCTTGCAAACCTGCCGTCAGCGCCAGCAATGCCAGTGCAGCCAAAGGCGTACGACGGATCATGGCTGCACCTCCAGTTGCTTGCGCAAGCGCTTGACCCGCGACTGCAGCACCTGCAGCTGCTCTTCG  >  NC_002947/4744595‑4744891

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