BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,745,828	G→T	7.5%	A672E (GCG→GAG)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,745,828	0	G	T	7.5%	479.9 / 19.7	147	A672E (GCG→GAG)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base G (54/82); new base T (6/5); total (60/87)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 151 total aligned reads displayed.

AAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAACACCACCTTCAGCAGGCCGCGC  >  NC_002947/5745679‑5745977

 aGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAgcgc                                                                                                                                                      <  2:1900520/150‑1 (MQ=255)

   cGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCAc                                                                                                                                                    <  1:3671916/150‑1 (MQ=255)

     gCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCa                                                                                                                                                  >  2:1523085/1‑150 (MQ=255)

      cTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGCCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCAt                                                                                                                                                 >  1:129589/1‑150 (MQ=255)

       tATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCCCGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCAtt                                                                                                                                                <  1:3890776/150‑1 (MQ=255)

       tATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGCGCTCACCAtt                                                                                                                                                <  1:2550665/150‑1 (MQ=255)

         tGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGa                                                                                                                                              <  2:266171/150‑1 (MQ=255)

         tGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGa                                                                                                                                              <  2:5207639/150‑1 (MQ=255)

         tGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGa                                                                                                                                              >  2:3148717/1‑150 (MQ=255)

          ggCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAg                                                                                                                                             <  1:3111158/150‑1 (MQ=255)

           gcagcaTTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGCGCTCACCATTggca                                                                                                                                            <  1:4762208/150‑4 (MQ=255)

           gcagcaTTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAgg                                                                                                                                            <  2:1776086/150‑1 (MQ=255)

           gcagcaTTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAgg                                                                                                                                            >  2:2245354/1‑150 (MQ=255)

           gcagcaTTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAgg                                                                                                                                            >  1:778609/1‑150 (MQ=255)

            cagcaTTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGc                                                                                                                                           >  1:2868150/1‑150 (MQ=255)

                aTTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGaa                                                                                                                                       >  2:753105/1‑150 (MQ=255)

                     cgcCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCa                                                                                                                                  <  2:594736/150‑1 (MQ=255)

                      gcCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCAt                                                                                                                                 >  1:5240655/1‑150 (MQ=255)

                      gcCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCAt                                                                                                                                 >  1:912890/1‑150 (MQ=255)

                      gcCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCAt                                                                                                                                 <  1:2556721/150‑1 (MQ=255)

                       ccAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATg                                                                                                                                <  1:3706561/150‑1 (MQ=255)

                           aaTCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGc                                                                                                                            <  2:1187779/150‑1 (MQ=255)

                           aaTCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCACGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGc                                                                                                                            >  2:3886581/1‑150 (MQ=255)

                             tCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTg                                                                                                                          <  1:723873/150‑1 (MQ=255)

                             tCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTg                                                                                                                          <  1:723956/150‑1 (MQ=255)

                               aCCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCCCCATTGAGGCCGAACTTCATGTTGCTGGt                                                                                                                        >  2:4901196/1‑150 (MQ=255)

                               aCCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGt                                                                                                                        <  2:4286192/150‑1 (MQ=255)

                                      cGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGaa                                                                                                                 <  1:1392129/150‑1 (MQ=255)

                                       gCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACTTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAg                                                                                                                <  2:2053336/150‑1 (MQ=255)

                                         ttGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCt                                                                                                              >  2:3366294/1‑150 (MQ=255)

                                          tGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCtt                                                                                                             <  2:4664927/150‑1 (MQ=255)

                                           gcgcgcTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTc                                                                                                            <  1:41160/150‑1 (MQ=255)

                                           gcgcgcTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCACACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTc                                                                                                            <  2:592091/150‑1 (MQ=255)

                                              cgcTTCAACCTGCTGCGAGGCCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTa                                                                                                         <  2:129589/150‑1 (MQ=255)

                                                cTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGTCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAAc                                                                                                       >  2:3876167/1‑150 (MQ=255)

                                                cTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCGGTTCGCACATCTCCCCGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAAc                                                                                                       <  1:4821659/150‑1 (MQ=255)

                                                 ttCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAAcc                                                                                                      >  2:3369729/1‑150 (MQ=255)

                                                 ttCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAAcc                                                                                                      >  1:3978685/1‑150 (MQ=255)

                                                 ttCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAAcc                                                                                                      >  2:2900096/1‑150 (MQ=255)

                                                 ttCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAAcc                                                                                                      >  2:2899459/1‑150 (MQ=255)

                                                     aCCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCcgcg                                                                                                  >  1:3015009/1‑150 (MQ=255)

                                                       ctgctgCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGt                                                                                                >  2:1794782/1‑150 (MQ=255)

                                                        tgctgcGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTg                                                                                               <  1:1477794/150‑1 (MQ=255)

                                                        tgctgcGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTg                                                                                               <  1:3149004/150‑1 (MQ=255)

                                                         gctgcGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTgg                                                                                              <  2:2972324/150‑1 (MQ=255)

                                                            gcGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGata                                                                                           <  2:4512497/150‑1 (MQ=255)

                                                             cGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGatat                                                                                          <  2:912701/150‑1 (MQ=255)

                                                                   cAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTc                                                                                    <  1:1712639/150‑1 (MQ=255)

                                                                      gCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTc                                                                                                            <  1:2040063/123‑1 (MQ=255)

                                                                      gCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGa                                                                                 >  1:4137766/1‑150 (MQ=255)

                                                                        cAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGaaa                                                                               >  1:3641037/1‑150 (MQ=255)

                                                                        cAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGaaa                                                                               >  1:2653103/1‑150 (MQ=255)

                                                                          aaGATGAACATGTTGTCGGACCCCACCTGTTCGCAAATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGAGAACCATTGAGGACGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACACGCGGTGGATATCTGCTCGGAAAgg                                                                             >  1:4692717/1‑150 (MQ=255)

                                                                           aGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGt                                                                            <  2:4303085/150‑1 (MQ=255)

                                                                              tGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTccg                                                                         >  1:1234747/1‑150 (MQ=255)

                                                                               gAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTccgc                                                                        >  2:794001/1‑150 (MQ=255)

                                                                               gAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTccgc                                                                        >  2:3535556/1‑150 (MQ=255)

                                                                               gAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTccgc                                                                        >  1:4226511/1‑150 (MQ=255)

                                                                                        gtCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGgat                                                               <  2:269418/150‑1 (MQ=255)

                                                                                         tCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGgatg                                                              <  1:915120/150‑1 (MQ=255)

                                                                                           ggCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGgatgat                                                            <  1:3362670/150‑1 (MQ=255)

                                                                                           ggCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGgatgat                                                            <  1:4418739/150‑1 (MQ=255)

                                                                                               cccACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGctc                                                        <  2:3732218/150‑1 (MQ=255)

                                                                                                 cACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCtctc                                                      >  2:3918445/1‑150 (MQ=255)

                                                                                                 cACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCtctc                                                      <  1:2054560/150‑1 (MQ=255)

                                                                                                 cACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCtctc                                                      >  1:4383534/1‑150 (MQ=255)

                                                                                                 cACCTGTTCGCACATCTCCACGTTGGCGCCAGCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGGTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGctctt                                                      <  2:3863291/150‑2 (MQ=255)

                                                                                                  aCCTGTTCGCTCATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCg                                                     <  2:891185/150‑1 (MQ=255)

                                                                                                  aCCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCg                                                     >  1:359249/1‑150 (MQ=255)

                                                                                                       ttCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAg                                                <  1:794001/150‑1 (MQ=255)

                                                                                                          gCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCt                                             <  1:5090038/150‑1 (MQ=255)

                                                                                                           cacaTCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTg                                            >  2:2604225/1‑150 (MQ=255)

                                                                                                                ctcCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGtt                                       >  1:5084008/1‑150 (MQ=255)

                                                                                                                ctcCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGtt                                       <  2:3641037/150‑1 (MQ=255)

                                                                                                                 tcCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTg                                      >  1:1413350/1‑150 (MQ=255)

                                                                                                                   cACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTa                                    >  2:4674713/1‑150 (MQ=255)

                                                                                                                       ttGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTg                                <  2:3776496/150‑1 (MQ=255)

                                                                                                                        tGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTgg                               <  2:1091926/150‑1 (MQ=255)

                                                                                                                        tGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTgg                               <  2:2152608/150‑1 (MQ=255)

                                                                                                                        tGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTgg                               <  2:257088/150‑1 (MQ=255)

                                                                                                                         ggCGCCATCCTGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTggg                              <  2:2368857/150‑1 (MQ=255)

                                                                                                                         ggCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTggg                              <  2:2369157/150‑1 (MQ=255)

                                                                                                                            gcCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAg                           >  1:1244301/1‑150 (MQ=255)

                                                                                                                            gcCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAg                           >  1:4008270/1‑150 (MQ=255)

                                                                                                                             ccATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAgg                          >  2:892198/1‑150 (MQ=255)

                                                                                                                             ccATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAgg                          <  1:3148717/150‑1 (MQ=255)

                                                                                                                               aTCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGaa                        <  1:3398555/150‑1 (MQ=255)

                                                                                                                               aTCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGaa                        <  1:2149349/150‑1 (MQ=255)

                                                                                                                                tCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCtctc                                                      <  1:1125373/119‑1 (MQ=255)

                                                                                                                                tCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAac                       <  2:3502384/150‑1 (MQ=255)

                                                                                                                                 ccAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAaca                      <  2:36713/150‑1 (MQ=255)

                                                                                                                                  cAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAtcac                     <  2:1413350/150‑1 (MQ=255)

                                                                                                                                   aGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAAcacc                    <  1:4662778/150‑1 (MQ=255)

                                                                                                                                        gCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAACACCACCtt               >  1:3889812/1‑150 (MQ=255)

                                                                                                                                        ctcAATGGTGCGCTCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGACCACCACCtg               >  2:1998343/3‑149 (MQ=255)

                                                                                                                                         tcAATGGTGCGCTCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAACACCACCTTc              >  1:2244637/2‑150 (MQ=255)

                                                                                                                                         tcAATGGTGCGCTCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAACACCACCTTc              >  2:2240158/2‑150 (MQ=255)

                                                                                                                                          cAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAACACCACCTTca             <  1:5391900/150‑1 (MQ=255)

                                                                                                                                           aaTGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAACACCACCTTcag            <  2:5125082/150‑1 (MQ=255)

                                                                                                                                                     gcACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAACACCACCTTCAGCAGGCCgcgc  >  1:4168106/1‑150 (MQ=255)

AAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCACCTGTTCGCACATCTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAACACCACCTTCAGCAGGCCGCGC  >  NC_002947/5745679‑5745977

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