BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,898,615	A→T	5.7%	L300* (TTA→TAA) ‡	PP_4307 ←	DMT superfamily permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,898,615	0	A	T	5.7%	314.4 / 8.3	106	L300* (TTA→TAA) ‡	PP_4307	DMT superfamily permease
Reads supporting (aligned to +/- strand): ref base A (48/51); new base T (2/4); total (50/56)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

Only 100 of 107 total aligned reads displayed.

CCGCGACCGTGGAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTCCACCAGAATTGCCGCAA  >  NC_002947/4898467‑4898762

ccGCGACCGTGGAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTaa                                                                                                                                                    <  1:2027153/150‑1 (MQ=255)

      ccGTGGAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCcgc                                                                                                                                              <  2:1533938/150‑1 (MQ=255)

         tGGAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTTTACCcgcgaa                                                                                                                                           <  1:810134/150‑3 (MQ=255)

         gggAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTTTACCcgcgaa                                                                                                                                           <  2:1669671/149‑3 (MQ=255)

          ggAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCt                                                                                                                                          >  1:1442691/1‑150 (MQ=255)

           gAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTa                                                                                                                                         <  1:2345938/150‑1 (MQ=255)

            aaCAGGCGGCGGGGTTCCTGTGGGGCAGGGTCGAGGGCTGCAAGCAGGTTTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTTTACCcgcgaagag                                                                                                                                        <  2:308013/150‑6 (MQ=255)

             aCAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACa                                                                                                                                       <  1:1036891/150‑1 (MQ=255)

              cAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTGAACCAGAGCCAATAGTAT‑GCTGATCGTAGTCTTTCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAg                                                                                                                                      <  1:2122432/150‑1 (MQ=255)

               aGGCAGAGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGAAAGCAGGTGTTGCTCAAGGGTG‑TAGATGGGCGAGGTCATGACGGATCGCAGGCTGGCAAAAAGCGGCATTCTAAATTGCGGGACTATTCGCGGGTAAACCCACCCCTACAgg                                                                                                                                     >  2:1698143/1‑150 (MQ=255)

                        gTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCgcgc                                                                                                                            <  2:1563626/150‑1 (MQ=255)

                         tttCCTGTGGCGCAGGGTCGAGGGCGGCAAGCAGTTGTTGCTCAAGGCTG‑TCGATGGGCTAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGcc                                                                                                                           <  1:708778/150‑1 (MQ=255)

                             cTGTGGGGCAGGGTAGAGGGCGGAAAGAAGATGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGATGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGc                                                                                                                       <  2:22085/150‑1 (MQ=255)

                              tgtgGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTTTACCcgcgaa                                                                                                                                           >  1:1465768/1‑127 (MQ=255)

                              tgtgGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTTTACCcgcgaa                                                                                                                                           <  2:1465768/129‑3 (MQ=255)

                              tgtgGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCg                                                                                                                      >  1:1019404/1‑150 (MQ=255)

                                  gggCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAac                                                                                                                  <  2:1025851/150‑1 (MQ=255)

                                  gggCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAac                                                                                                                  >  2:1389129/1‑150 (MQ=255)

                                  gggCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAac                                                                                                                  >  1:313256/1‑150 (MQ=255)

                                  gggCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAac                                                                                                                  <  1:1624396/150‑1 (MQ=255)

                                   ggCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAacc                                                                                                                 >  1:1474039/1‑150 (MQ=255)

                                   ggCAGGGTCGAGGGCGGAAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAacc                                                                                                                 <  2:688092/150‑1 (MQ=255)

                                        ggTCGATGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAaccaccac                                                                                                            >  2:906492/1‑150 (MQ=255)

                                        ggTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCcgc                                                                                                                                              >  2:356779/1‑116 (MQ=255)

                                        ggTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCcgc                                                                                                                                              <  1:356779/116‑1 (MQ=255)

                                         gTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCCGGCAAAAAGCGGCATTCTACCTTGCTGGCCTCTCCGCGGGCCATCCCGCGCCTACAGGATGCCGCGCCGAGCGGAacctccacc                                                                                                           <  2:338513/150‑1 (MQ=255)

                                          tCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCAccc                                                                                                          <  2:2184047/150‑1 (MQ=255)

                                          tCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCAccc                                                                                                          >  1:2214352/1‑150 (MQ=255)

                                           cGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAaccac                                                                                                               <  2:128118/144‑1 (MQ=255)

                                           cGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAaccac                                                                                                               >  1:128118/1‑144 (MQ=255)

                                           cGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCCACAGGATGCCGCGCCGAGCGGAACCACCTcccc                                                                                                         <  2:1273709/150‑1 (MQ=255)

                                              gggCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCctgc                                                                                                      <  1:1435003/150‑1 (MQ=255)

                                                  ggcAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCa                                                                                                  >  1:1965752/1‑150 (MQ=255)

                                                       gCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGccc                                                                                             <  1:570741/150‑1 (MQ=255)

                                                        cAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGAATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCg                                                                                            >  2:1326417/1‑150 (MQ=255)

                                                         aGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGc                                                                                           <  1:1931978/150‑1 (MQ=255)

                                                             gtTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCAt                                                                                                 >  1:1831920/1‑140 (MQ=255)

                                                             gtTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCAt                                                                                                 <  2:1831920/140‑1 (MQ=255)

                                                                  tCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCcag                                                                                  >  2:35275/1‑150 (MQ=255)

                                                                  tCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCcag                                                                                  <  1:503077/150‑1 (MQ=255)

                                                                    aaGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCcgc                                                                                                                              >  1:1791156/1‑104 (MQ=255)

                                                                    aaGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCcgc                                                                                                                              <  2:1791156/104‑1 (MQ=255)

                                                                    aaGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCcagca                                                                                >  2:480893/1‑150 (MQ=255)

                                                                       ggTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGa                                                                             >  2:2081759/1‑150 (MQ=255)

                                                                         tg‑tCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATc                                                                           >  2:655415/1‑150 (MQ=255)

                                                                             cGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAac                                                                        <  1:2223331/150‑1 (MQ=255)

                                                                              gATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAacc                                                                       >  2:859861/1‑150 (MQ=255)

                                                                               aTGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACa                                                                                                                                       >  2:2378575/1‑85 (MQ=255)

                                                                               aTGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACa                                                                                                                                       <  1:2378575/85‑1 (MQ=255)

                                                                                   gCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAaccccacc                                                                  <  1:480893/150‑1 (MQ=255)

                                                                                     gAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAaccccacccc                                                                >  2:1580455/1‑150 (MQ=255)

                                                                                           aTGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCa                                                          <  1:2173961/150‑1 (MQ=255)

                                                                                           aTGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCa                                                          >  1:1412062/1‑150 (MQ=255)

                                                                                           aTGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCa                                                          >  1:1414252/1‑150 (MQ=255)

                                                                                            tGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCAc                                                         <  2:1414252/150‑1 (MQ=255)

                                                                                            tGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCAc                                                         >  2:1045358/1‑150 (MQ=255)

                                                                                             ggCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACt                                                        <  1:906492/150‑1 (MQ=255)

                                                                                                ggATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTa                                                     >  1:1602194/1‑150 (MQ=255)

                                                                                                ggATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTa                                                     >  2:717256/1‑150 (MQ=255)

                                                                                                  aTCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAgc                                                   <  1:2080373/150‑1 (MQ=255)

                                                                                                     gCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCTAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAgcggc                                                <  1:35275/150‑1 (MQ=255)

                                                                                                      cAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGTCACTGTAgcggcg                                               >  2:550033/1‑150 (MQ=255)

                                                                                                      cAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAgcggcg                                               >  2:1686310/1‑150 (MQ=255)

                                                                                                       aGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCTCCCCTTCTGCATGCCCGCTGCCGCCAGCAAGATCAACCCCACCCCCAGCCACTGTAGCGgcgc                                              <  1:655415/150‑1 (MQ=255)

                                                                                                           tGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAgg                                          >  2:1451194/1‑150 (MQ=255)

                                                                                                           tGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGAGGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAgg                                          >  1:424484/1‑150 (MQ=255)

                                                                                                             gCAAAAAGCGGCATTCTACCTTGCGGGCGTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCg                                        <  2:177957/150‑1 (MQ=255)

                                                                                                               aacaaGCGGTATTCTACCTTGCGGGACTCCGCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGTATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGt                                      <  1:604305/147‑1 (MQ=255)

                                                                                                               aaaaaGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCggg                                      >  2:1352998/1‑149 (MQ=255)

                                                                                                               aaaaaGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGt                                      >  2:2298179/1‑150 (MQ=255)

                                                                                                               aaaaaGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGt                                      <  1:1045358/150‑1 (MQ=255)

                                                                                                               aaaaaGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGt                                      >  2:1347174/1‑150 (MQ=255)

                                                                                                                aaaaGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTg                                     <  1:1336290/150‑1 (MQ=255)

                                                                                                                 aaaGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGc                                    <  1:1580455/150‑1 (MQ=255)

                                                                                                                 aaaGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGc                                    <  1:1451194/150‑1 (MQ=255)

                                                                                                                      ggCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGaa                               <  1:253657/150‑1 (MQ=255)

                                                                                                                      ggCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGaa                               <  2:1474039/150‑1 (MQ=255)

                                                                                                                       gCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAAc                              <  1:633552/150‑1 (MQ=255)

                                                                                                                        cATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACACGTTCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACg                             >  1:624877/1‑150 (MQ=255)

                                                                                                                          ttCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCca                                                                                   <  1:2266326/94‑1 (MQ=255)

                                                                                                                          ttCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCca                                                                                   >  2:2266326/1‑94 (MQ=255)

                                                                                                                              aCCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCAccc                       >  1:1665550/1‑150 (MQ=255)

                                                                                                                               ccTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCa                      >  2:2253546/1‑150 (MQ=255)

                                                                                                                                cTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAg                     <  2:1965752/150‑1 (MQ=255)

                                                                                                                                 ttGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGt                    <  1:1326417/150‑1 (MQ=255)

                                                                                                                                  tGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTc                   <  1:717256/150‑1 (MQ=255)

                                                                                                                                     gggCCTGTGCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCTCCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTccac                <  1:716883/150‑1 (MQ=255)

                                                                                                                                     gggCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCca                                                                                   <  2:2160892/83‑1 (MQ=255)

                                                                                                                                     gggCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCca                                                                                   >  1:2160892/1‑83 (MQ=255)

                                                                                                                                     gggCCTCTTCGCCGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTccac                <  1:1376331/150‑1 (MQ=255)

                                                                                                                                       gCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTccacca              >  1:2230459/1‑150 (MQ=255)

                                                                                                                                       gCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGACCACACGCGACCCAGTccacca              >  2:1022258/1‑150 (MQ=255)

                                                                                                                                        ccTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGACCAGCCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTACCCGAACGCCACCCAGTCCACCAg             >  1:1240798/1‑150 (MQ=255)

                                                                                                                                         ctctTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGACGAGCGGCACCACCACCCCTGCTGCATGCCCGCTACCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTCCACCAGa            >  1:1616512/1‑150 (MQ=255)

                                                                                                                                          tctTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGAGGTGCCCGAACGCCACCCAGTCCACCAGaa           <  2:1412062/150‑1 (MQ=255)

                                                                                                                                            ttCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCAc                         <  1:598834/134‑1 (MQ=255)

                                                                                                                                            ttCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCAc                         >  2:598834/1‑134 (MQ=255)

                                                                                                                                              cgcgGGTAAACCCGCGCCTACCGGATGCCGCACAGAGCGGCCCCCCCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTCCACCAGAATTGc       <  2:2230459/150‑1 (MQ=255)

                                                                                                                                               gcgGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGTCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTCCACCAGAATTGcc      <  1:550033/150‑1 (MQ=255)

                                                                                                                                                   gTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCAccc                                                                                                          <  1:2269832/46‑1 (MQ=255)

CCGCGACCGTGGAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTG‑TCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTCCACCAGAATTGCCGCAA  >  NC_002947/4898467‑4898762

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