BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	552,256	T→A	6.8%	L186M (TTG→ATG)	rplD →	50S ribosomal protein L4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	552,256	0	T	A	6.8%	297.2 / 7.0	101	L186M (TTG→ATG)	rplD	50S ribosomal protein L4
Reads supporting (aligned to +/- strand): ref base T (51/43); new base A (2/5); total (53/48)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

GCTGAACGGCATGGGTCTGAGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCAGTTCGTATTCAAGGTTGCTACT  >  NC_002947/552108‑552404

gCTGAACGGCATGGGTCTGAGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAgtgt                                                                                                                                                      >  5:41640/1‑149 (MQ=255)

gCTGAACGGCATGGGTCTGAGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAgtgt                                                                                                                                                      <  6:241610/149‑1 (MQ=255)

    aaCGGCATGGGTCTGAGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTg                                                                                                                                                    <  2:218148/147‑1 (MQ=255)

         cATGGGTCTGAGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCAGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTg                                                                                                                                             >  5:405063/1‑149 (MQ=255)

          aTGGGTCTGAGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTg                                                                                                                                             >  8:139960/1‑148 (MQ=255)

             ggTCTGAGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGtctc                                                                                                                                         >  1:2268/1‑149 (MQ=255)

              gTCTGAGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCg                                                                                                                                        <  8:71279/149‑1 (MQ=255)

               tCTGAGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGACAGTGATCAACACTTtct                                                                                                                                          <  7:186517/146‑1 (MQ=255)

               tCTGAGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGACAGTGATCAACACTTtct                                                                                                                                          >  8:186517/1‑146 (MQ=255)

                cTGAGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGc                                                                                                                                      >  8:247008/1‑149 (MQ=255)

                 tGAGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGcc                                                                                                                                     >  8:5663/1‑149 (MQ=255)

                   aGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGACAGTGATCAACACTTTCTCgtatg                                                                                                                                    <  2:66611/148‑5 (MQ=255)

                   aGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGACAGTGATCAACACTTTCTCgtat                                                                                                                                     <  1:413913/147‑4 (MQ=255)

                   aGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGACAGTGATCAACACTTTCTCgtat                                                                                                                                     <  5:261510/147‑4 (MQ=255)

                   aGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGACAGTGATCAACACTTTCTCgtat                                                                                                                                     <  7:14757/147‑4 (MQ=255)

                   aGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGcc                                                                                                                                     <  5:188276/147‑1 (MQ=255)

                     cGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTga                                                                                                                                 <  2:370631/149‑1 (MQ=255)

                           aCTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGaaa                                                                                                                            >  1:183869/1‑148 (MQ=255)

                             tGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTc                                                                                                                         >  1:160904/1‑149 (MQ=255)

                               aTCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCga                                                                                                                       <  2:156207/149‑1 (MQ=255)

                                  gTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCgagga                                                                                                                    <  5:364283/149‑1 (MQ=255)

                                  gTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCgag                                                                                                                      >  7:51457/1‑147 (MQ=255)

                                  gTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCgag                                                                                                                      >  8:194249/1‑147 (MQ=255)

                                  gTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCgag                                                                                                                      >  1:205293/1‑147 (MQ=255)

                                        gACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAgctgct                                                                                                              >  6:67008/1‑149 (MQ=255)

                                           gCTGTTGATCAGATCGTGTAGCTGGCTGCTCGTAACGTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACGCGGT‑CAGTCTGTTCGCATTGGTGATAGTGTTGATCAGTGTCTCGTCCGTGAAGAAATTCGAGGAGCTGCTggg                                                                                                           <  6:82360/149‑1 (MQ=255)

                                            cTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGa                                                                                                          <  7:194249/149‑1 (MQ=255)

                                              gTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTCTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGAt                                                                                                         >  5:440905/1‑148 (MQ=255)

                                               ttGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCAGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTAGGACGTGAAGAAATTCGAGGAGCTGCTGGGATg                                                                                                        >  6:328074/1‑148 (MQ=255)

                                               ttGATCAGAAACTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATg                                                                                                        >  6:123703/1‑148 (MQ=255)

                                                 gATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGACGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGaa                                                                                                      <  8:159640/148‑1 (MQ=255)

                                                  aTCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAAcc                                                                                                    <  8:245222/149‑1 (MQ=255)

                                                  aTCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAAc                                                                                                     >  6:172412/1‑148 (MQ=255)

                                                  aTCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAAc                                                                                                     <  4:221413/148‑1 (MQ=255)

                                                   tCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAAc                                                                                                     >  2:62693/1‑147 (MQ=255)

                                                   tCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTACGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAAc                                                                                                     >  8:230036/1‑147 (MQ=255)

                                                      gAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGctgctg                                                                                                             <  5:9986/136‑1 (MQ=255)

                                                      gAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGctgctg                                                                                                             >  6:9986/1‑136 (MQ=255)

                                                         ccTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAAcg                                                                                             >  1:8759/1‑149 (MQ=255)

                                                          cTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAAcgc                                                                                            >  5:391486/1‑149 (MQ=255)

                                                           tGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAAcgcg                                                                                           >  8:368410/1‑149 (MQ=255)

                                                               ccTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTAt                                                                                        >  6:5312/1‑148 (MQ=255)

                                                                cTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTAtt                                                                                       >  1:11243/1‑148 (MQ=255)

                                                                   gctgctCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCaaa                                                                                   <  8:51457/149‑1 (MQ=255)

                                                                    ctgctCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAg                                                                                  <  3:43780/149‑1 (MQ=255)

                                                                    ctgctCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAg                                                                                  <  2:183869/149‑1 (MQ=255)

                                                                     tgctCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGt                                                                                 >  8:392328/1‑149 (MQ=255)

                                                                     tgctCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGATGGGATGAACCAGGAACGCGTATTCAAAGt                                                                                 >  6:277516/1‑149 (MQ=255)

                                                                      gctCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGtc                                                                                >  4:89781/1‑149 (MQ=255)

                                                                       ctCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGtcc                                                                               <  5:405698/149‑1 (MQ=255)

                                                                        tCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGtcct                                                                              >  4:58779/1‑149 (MQ=255)

                                                                          gTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGtcctc                                                                             <  2:205293/148‑1 (MQ=255)

                                                                               cTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACTAGAAAGTGTTGATCACTGTCTCGGCCTTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGc                                                                       <  5:328074/149‑1 (MQ=255)

                                                                                  ccGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGCTCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCg                                                                    <  8:410536/149‑1 (MQ=255)

                                                                                    gCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCa                                                                  >  8:39594/1‑149 (MQ=255)

                                                                                     cACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCAc                                                                 <  1:62693/149‑1 (MQ=255)

                                                                                     cACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCAc                                                                 <  5:67741/149‑1 (MQ=255)

                                                                                        gTCGACGTACGTGACGTCCAAGGTTCCGACCCGGTACAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGt                                                               >  1:284262/1‑149 (MQ=255)

                                                                                             cgtacgtGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCga                                                         <  7:5663/149‑1 (MQ=255)

                                                                                             cgtacgtGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCga                                                         <  2:284262/149‑1 (MQ=255)

                                                                                               tacgtGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCgaga                                                       <  8:47973/149‑1 (MQ=255)

                                                                                                acgtGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACAGTCTCGGCCGTGAAGAAATTCGATGAGCTGCTGGGATGAAACAGGAACGCGTATACAAAGTCCTCCTTGGCCCGCACGTATCCgag                                                        >  3:304295/1‑147 (MQ=255)

                                                                                                        tCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCg                                              >  2:251641/1‑149 (MQ=255)

                                                                                                        tCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCg                                              >  5:380562/1‑149 (MQ=255)

                                                                                                          cAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGtt                                            >  5:46465/1‑149 (MQ=255)

                                                                                                          cAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGt                                             >  5:264438/1‑148 (MQ=255)

                                                                                                              gTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTgg                                        <  3:58779/149‑1 (MQ=255)

                                                                                                               ttCCGACCCGGT‑CAGTCTGATCGCATCCGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAAGGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTgg                                        <  5:123703/148‑1 (MQ=255)

                                                                                                               ttCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGc                                       >  3:420060/1‑149 (MQ=255)

                                                                                                               ttCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGc                                       >  7:161198/1‑149 (MQ=255)

                                                                                                                 ccGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCt                                      <  2:2268/148‑1 (MQ=255)

                                                                                                                  cGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTga                                    <  4:420060/149‑1 (MQ=255)

                                                                                                                  cGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTga                                    <  7:247008/149‑1 (MQ=255)

                                                                                                                   gACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCt                                      >  4:408666/1‑146 (MQ=255)

                                                                                                                      ccGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGaa                                 >  2:317382/1‑148 (MQ=255)

                                                                                                                      ccGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGaa                                 >  4:98573/1‑148 (MQ=255)

                                                                                                                      ccGGT‑CAGTCTGATCGCATACGAGAAAGAGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGaa                                 >  1:242708/1‑148 (MQ=255)

                                                                                                                        ggT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGaaaaa                              <  6:391486/149‑1 (MQ=255)

                                                                                                                          t‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAg                            >  5:161620/1‑149 (MQ=255)

                                                                                                                          t‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAg                            >  4:314255/1‑149 (MQ=255)

                                                                                                                             aGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGc                          <  1:317382/149‑1 (MQ=255)

                                                                                                                              gTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGcc                         <  3:89781/149‑1 (MQ=255)

                                                                                                                               tCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTAGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGAACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCa                        >  1:248850/1‑149 (MQ=255)

                                                                                                                               tCTGATCGCATACGAGAAAGTGATGATCACTGTCTCGGCCGTGAAGAAAATCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTAAAAGTCCTCCTTGGCCCGCACGTTTCAGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGc                          >  4:68153/1‑147 (MQ=255)

                                                                                                                                cTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCa                        <  5:172412/148‑1 (MQ=255)

                                                                                                                                 tGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCAGt                      <  7:368410/149‑1 (MQ=255)

                                                                                                                                 tGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCAGt                      <  2:242708/149‑1 (MQ=255)

                                                                                                                                    tCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCAGTTc                    <  2:160904/148‑1 (MQ=255)

                                                                                                                                      gCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCAGTTCGTa                 <  5:5312/149‑1 (MQ=255)

                                                                                                                                        aTACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCAGTTCGTAtt               <  6:440905/149‑1 (MQ=255)

                                                                                                                                        aTACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCAGTTCGTAt                >  1:49695/1‑148 (MQ=255)

                                                                                                                                          aCGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTga                                    <  1:229237/126‑1 (MQ=255)

                                                                                                                                          aCGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTga                                    >  2:229237/1‑126 (MQ=255)

                                                                                                                                            gagaAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCAGTTCGTATTCAAg           >  1:111672/1‑149 (MQ=255)

                                                                                                                                            gagaAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCAGTTCGTATTCAAg           >  6:334109/1‑149 (MQ=255)

                                                                                                                                            gagaAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCAGTTCGTATTCAAg           >  4:391519/1‑149 (MQ=255)

                                                                                                                                             agaAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCAGTTCGTATTCAAgg          <  2:8759/149‑1 (MQ=255)

                                                                                                                                              gaAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGTCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCAGTTCGTATTCAAGGt         >  7:83294/1‑149 (MQ=255)

                                                                                                                                                 aGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCAGTTCGTATTCAAGGTTGc      <  7:39594/149‑1 (MQ=255)

                                                                                                                                                     ttGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCAGTTCGTATTCAAGGTTGCTACt  >  6:43297/1‑149 (MQ=255)

GCTGAACGGCATGGGTCTGAGCGATGTACTGATCGTTTCCGACGCTGTTGATCAGAACCTGTACCTGGCTGCTCGTAACCTGCCGCACGTCGACGTACGTGACGTCCAAGGTTCCGACCCGGT‑CAGTCTGATCGCATACGAGAAAGTGTTGATCACTGTCTCGGCCGTGAAGAAATTCGAGGAGCTGCTGGGATGAACCAGGAACGCGTATTCAAAGTCCTCCTTGGCCCGCACGTTTCCGAGAAGGCTACCGTTCTGGCTGAGAAAAAAGGCCAGTTCGTATTCAAGGTTGCTACT  >  NC_002947/552108‑552404

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