BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,541,769	T→A	7.7%	E326D (GAA→GAT)	dnaB ←	replicative DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,541,769	0	T	A	7.7%	290.4 / 9.7	104	E326D (GAA→GAT)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand): ref base T (52/43); new base A (4/4); total (56/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 106 total aligned reads displayed.

AGCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCG  >  NC_002947/5541620‑5541916

aGCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGAtt                                                                                                                                                      >  2:1574034/1‑150 (MQ=255)

 gCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGAttt                                                                                                                                                     <  1:2230243/150‑1 (MQ=255)

  cTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTc                                                                                                                                                    >  1:2325752/1‑150 (MQ=255)

       aaaGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCAt                                                                                                                                               <  1:401324/150‑1 (MQ=255)

       aaaGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCAt                                                                                                                                               <  2:1624771/150‑1 (MQ=255)

         aGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGc                                                                                                                                             >  2:521/1‑150 (MQ=255)

          gggCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCt                                                                                                                                            >  2:1020076/1‑150 (MQ=255)

           ggCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTc                                                                                                                                           >  2:1739042/1‑150 (MQ=255)

           ggCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGAACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGAGCATGGCGAAATCGCCATGATc                                                                                                                                           <  2:484584/150‑1 (MQ=255)

             cAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAACTCTCGGAAAATTCGTTGGTACGGTTGTCACCGGCCGCACCGGGAATCTGCATCGGCTGCAGGTAGTCGACCATGATCATTGCGTTTTCGCCATACTCAc                                                                                                                                         >  1:2010148/1‑150 (MQ=255)

                   ccGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCa                                                                                                                                   >  1:2031352/1‑150 (MQ=255)

                    ctggCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAg                                                                                                                                  >  1:1759334/3‑150 (MQ=255)

                    cGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGAGCATGGCGAAATCGCCATGATc                                                                                                                                           <  2:759427/141‑1 (MQ=255)

                    cGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCCGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGAGCATGGCGAAATCGCCATGATc                                                                                                                                           >  1:759427/1‑141 (MQ=255)

                         gcttGAATTCTTGGGCCCGGGCCTTGCGAGCGCGGGCCATCTCGATGACGTCATTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAg‑gcggc                                                                                                                            <  2:2231424/148‑1 (MQ=255)

                         aGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAg‑gcggc                                                                                                                            >  2:740133/1‑150 (MQ=255)

                         aGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCAGGAATCTGCATCAGATGCAGGTAGTCAACCATGATCATGGCGAGTGCGCCATGCTCACTGGCCAg‑gcgga                                                                                                                            >  2:623580/1‑149 (MQ=255)

                           ttGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑Gcggcgg                                                                                                                          <  2:230540/150‑1 (MQ=255)

                           ttGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATATCGTAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑Gcggcgg                                                                                                                          <  1:1592595/150‑1 (MQ=255)

                           ttGAATTCTTTGGCAAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGAAATCGCCATga                                                                                                                                             <  1:423584/132‑2 (MQ=255)

                            tGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCggg                                                                                                                         >  1:997049/1‑150 (MQ=255)

                             gAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGAGCATGGCGATTTCGCCATGCTCACGGGCCAG‑gcggcgaga                                                                                                                        >  1:2231424/1‑147 (MQ=255)

                              aaTTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTa                                                                                                                       >  1:1608374/1‑150 (MQ=255)

                              aaTTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTa                                                                                                                       >  2:2373212/1‑150 (MQ=255)

                                ttCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTAcg                                                                                                                     >  2:1438310/1‑150 (MQ=255)

                                   tttGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTAcgcgc                                                                                                                  >  2:1810398/1‑150 (MQ=255)

                                   tttGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTAcgcgc                                                                                                                  >  1:1445726/1‑150 (MQ=255)

                                          agtgCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCAtt                                                                                                           >  1:898866/4‑150 (MQ=255)

                                            ggCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTc                                                                                                         >  2:372818/1‑150 (MQ=255)

                                            ggCCGTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTc                                                                                                         <  2:1131089/150‑1 (MQ=255)

                                               ttttagcgagcgTGAAATCTCGTAAATTTCGTTTGTACGGTTGTCACCGCCCCACCCCCGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGACTTCGCAATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGa                                                                                                      <  2:2010148/146‑1 (MQ=255)

                                               cTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCACCTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGa                                                                                                      >  1:144910/1‑150 (MQ=255)

                                                ttGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCCCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAg                                                                                                     <  2:1103938/150‑1 (MQ=255)

                                                   agcgagcgGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAgggg                                                                                                  <  2:898888/150‑1 (MQ=255)

                                                          gggAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATAc                                                                                           >  1:1020775/1‑150 (MQ=255)

                                                            gAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGCGTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATAccc                                                                                         >  2:1035174/1‑150 (MQ=255)

                                                                  tcGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTa                                                                                   >  1:124958/1‑150 (MQ=255)

                                                                     gAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑Gcggcgg                                                                                                                          <  2:1854804/108‑1 (MQ=255)

                                                                     gAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑Gcggcgg                                                                                                                          <  2:1857840/108‑1 (MQ=255)

                                                                     gAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑Gcggcgg                                                                                                                          >  1:1857840/1‑108 (MQ=255)

                                                                     gAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑Gcggcgg                                                                                                                          >  1:1854804/1‑108 (MQ=255)

                                                                       aaTGTCGTTGGTACGGTTGGCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGGGGGGCTGATACCCGCCGTAtcgtc                                                                              <  2:936091/150‑1 (MQ=255)

                                                                         tttCGTGGGTACGGTTGTCACCGGCCGCCCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGa                                                                            <  1:152260/150‑1 (MQ=255)

                                                                          ttCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAg‑gc                                                                                                                               >  1:1288883/1‑98 (MQ=255)

                                                                          ttCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAg‑gc                                                                                                                               <  2:1288883/98‑1 (MQ=255)

                                                                           tCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGAt                                                                           >  1:1389497/1‑149 (MQ=255)

                                                                           tCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGAt                                                                           <  2:1389497/149‑1 (MQ=255)

                                                                             gTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGaa                                                                        >  2:28543/1‑150 (MQ=255)

                                                                                 gTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAAt                                                                    >  2:957658/1‑150 (MQ=255)

                                                                                  tACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAAtt                                                                   <  1:1388497/150‑1 (MQ=255)

                                                                                  tACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAAtt                                                                   <  1:1387873/150‑1 (MQ=255)

                                                                                     ggTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTc                                                                >  2:2079957/1‑150 (MQ=255)

                                                                                     ggTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTc                                                                >  2:441885/1‑150 (MQ=255)

                                                                                      gtggTCACCGGCCGACCCCGGAATCTGCATCACCTGCAGGTAGTCGACCCTGATCACGGCGATTTGGCCATGCTCACGGGCCAG‑GCGGCGGGTACGAGCGCGCATTTCCGAGGGGCTGGTACCCGCCGTATCGTCGATGAGCAATTTTCg                                                               >  2:407995/4‑150 (MQ=255)

                                                                                       ttGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCgg                                                              >  2:599384/1‑150 (MQ=255)

                                                                                         gTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTc                                                            <  1:1020076/150‑1 (MQ=255)

                                                                                         gTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTc                                                            <  1:740133/150‑1 (MQ=255)

                                                                                                  cGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAg                                                   <  2:1445726/150‑1 (MQ=255)

                                                                                                   gACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAgg                                                  <  1:1438310/150‑1 (MQ=255)

                                                                                                    aCCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGt                                                 <  1:957658/150‑1 (MQ=255)

                                                                                                     ccccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGtt                                                <  2:124958/150‑1 (MQ=255)

                                                                                                      cccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAg                                                                                                     <  1:1019217/96‑1 (MQ=255)

                                                                                                       ccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCa                                              >  1:1541441/1‑150 (MQ=255)

                                                                                                       ccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCa                                              <  1:28543/150‑1 (MQ=255)

                                                                                                       ccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCa                                              >  2:1948803/1‑150 (MQ=255)

                                                                                                         ggAATCTGCATCAGCTGCAGGTAGCAGACCAAGATCCGGCCGATGTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACt                                            <  1:1739042/150‑1 (MQ=255)

                                                                                                              cTGCATCAGCTGCAGGTATTCGACCTTCATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGa                                       <  1:521/150‑1 (MQ=255)

                                                                                                               tGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAg                                      >  1:2087914/1‑150 (MQ=255)

                                                                                                               tGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAg                                      >  2:1392140/1‑150 (MQ=255)

                                                                                                                 cATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAgggg                                                                                                  >  2:770383/1‑88 (MQ=255)

                                                                                                                 cATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAgggg                                                                                                  <  1:770383/88‑1 (MQ=255)

                                                                                                                  aTCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAg                                                   >  2:822821/1‑134 (MQ=255)

                                                                                                                  aTCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAg                                                   <  1:822821/134‑1 (MQ=255)

                                                                                                                      gCTGCATGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCAC‑GGCAAGAGCGTCGGGTACGCGCGCGCATTTCCGGGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGgc                               <  1:2367749/150‑1 (MQ=255)

                                                                                                                      gCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGa                                                                            <  1:989259/105‑1 (MQ=255)

                                                                                                                      gCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGgc                               <  2:898866/150‑1 (MQ=255)

                                                                                                                      gCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGgc                               >  1:1000078/1‑150 (MQ=255)

                                                                                                                           aGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGc                          >  1:432639/1‑150 (MQ=255)

                                                                                                                           aGGTAGGCGACCATGAGCATGGCGATTTCCTCATGCTCACGGGCCGG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGTAGGTTCACTGCCGAGGTCAGGCGCGGc                          <  2:785746/150‑1 (MQ=255)

                                                                                                                            ggTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACccgcc                                                                                      <  1:1684325/89‑1 (MQ=255)

                                                                                                                             gTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCa                        >  1:2129469/1‑150 (MQ=255)

                                                                                                                                ggcGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtc                     <  2:1026757/148‑1 (MQ=255)

                                                                                                                                   gACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtc                  <  2:432639/150‑1 (MQ=255)

                                                                                                                                   gACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtc                  >  1:1482668/1‑150 (MQ=255)

                                                                                                                                   gACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtc                     <  1:1201114/147‑1 (MQ=255)

                                                                                                                                   gACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtc                     >  2:1201114/1‑147 (MQ=255)

                                                                                                                                      cATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGGCAGGCGCGGCCAgtcgtcgtc               >  2:1762021/1‑150 (MQ=255)

                                                                                                                                       aTGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCCATTCCGAGGGGATGATACCCGACCTATCGTGCATGAACAATTATCGGTCGTTGAGGAGGTTAAGTGCAGAGGTCAGGCGCGGCCAgtcgtagtcg              >  2:602089/1‑150 (MQ=255)

                                                                                                                                        tGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgtcgt             <  2:1541441/150‑1 (MQ=255)

                                                                                                                                         gagcATGGCGAAATCGCCATGATCACGGGTCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgtcgtc            >  1:822165/4‑150 (MQ=255)

                                                                                                                                          aTCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTcc           <  1:2373212/150‑1 (MQ=255)

                                                                                                                                            cATGGCGAAATCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGa                                       <  1:1498869/120‑1 (MQ=255)

                                                                                                                                            cATGGCGAAATCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGa                                       >  2:1498869/1‑120 (MQ=255)

                                                                                                                                             aTGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGc        >  2:1694790/1‑150 (MQ=255)

                                                                                                                                              tGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCg                                                               <  2:2448121/94‑1 (MQ=255)

                                                                                                                                              tGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCg                                                               >  1:2448121/1‑94 (MQ=255)

                                                                                                                                              tGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGt                                                          <  1:422395/99‑1 (MQ=255)

                                                                                                                                              tGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGt                                                          >  2:422395/1‑99 (MQ=255)

                                                                                                                                                   aTTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCg  <  1:599384/150‑1 (MQ=255)

                                                                                                                                                   aTTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCg  <  2:1482668/150‑1 (MQ=255)

AGCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAG‑GCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCG  >  NC_002947/5541620‑5541916

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