BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,282,344	A→G	5.7%	L257L (TTG→CTG)	ampC ←	beta‑lactamase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,282,344	0	A	G	5.7%	842.8 / 23.2	262	L257L (TTG→CTG)	ampC	beta‑lactamase
Reads supporting (aligned to +/- strand): ref base A (110/137); new base G (9/6); total (119/143)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 270 total aligned reads displayed.

GCCATCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAAGGTAGG  >  NC_002947/3282195‑3282493

gCCATCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCaa                                                                                                                                                       <  2:2205810/150‑1 (MQ=255)

  cATCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGGTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAAcc                                                                                                                                                     <  2:4837082/150‑1 (MQ=255)

   aTCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCg                                                                                                                                                    >  2:2282782/1‑150 (MQ=255)

    tCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGc                                                                                                                                                   >  2:1136909/1‑150 (MQ=255)

    tCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGc                                                                                                                                                   <  2:4814970/150‑1 (MQ=255)

    tCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGGTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGc                                                                                                                                                   <  2:5277985/150‑1 (MQ=255)

     cGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGca                                                                                                                                                  <  2:4584749/150‑1 (MQ=255)

      gATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGTTGCACATGCAGCCGcag                                                                                                                                                 <  1:1000488/150‑2 (MQ=255)

      gATCAGTGGGCTCAGCGTGACCGGGTAGGGGTAGCGGTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGTTGCACATGCAGCCGcag                                                                                                                                                 <  1:3829295/150‑2 (MQ=255)

       aTCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGcaca                                                                                                                                                <  1:4367236/150‑1 (MQ=255)

       aTCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGcaca                                                                                                                                                <  1:2237862/150‑1 (MQ=255)

          aGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAg                                                                                                                                             >  1:2979693/1‑150 (MQ=255)

            tgtgCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAgcg                                                                                                                                           <  2:2803975/150‑1 (MQ=255)

               gCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCgcag                                                                                                                                        <  1:2837170/150‑1 (MQ=255)

                cTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCgcagc                                                                                                                                       <  2:3738838/150‑1 (MQ=255)

                cTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCgcagc                                                                                                                                       >  2:4293210/1‑150 (MQ=255)

                cTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCgcagc                                                                                                                                       <  1:3128886/150‑1 (MQ=255)

                cTCAGCGTGACAGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCgcagc                                                                                                                                       >  1:3970005/1‑150 (MQ=255)

                 tCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCgcagca                                                                                                                                      >  1:2418993/1‑150 (MQ=255)

                  cAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGcagcag                                                                                                                                     >  2:2323783/1‑150 (MQ=255)

                   aGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAgg                                                                                                                                    <  2:3490933/150‑1 (MQ=255)

                     cGTGACCGGGTAGGGGTAGCGTTCCCAGCCCGGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTc                                                                                                                                  <  2:5110642/150‑1 (MQ=255)

                       tGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTcgc                                                                                                                                >  2:5624285/1‑150 (MQ=255)

                       tGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTcgc                                                                                                                                >  1:2357606/1‑150 (MQ=255)

                         aCCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTAGCGTCATGGCGCCAACCTGGAAGAAGCCCTGCTGGGTAATGGCGGTGGCTTACGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTcgctg                                                                                                                              <  2:57186/150‑1 (MQ=255)

                           cGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTcgctggc                                                                                                                            <  1:1617319/150‑1 (MQ=255)

                             ggtaggggtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTcgctggcgc                                                                                                                          <  2:1704713/150‑1 (MQ=255)

                               taggggtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTcgctggcgctg                                                                                                                        <  1:4587277/150‑1 (MQ=255)

                                aggggtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTcgctggcgctgg                                                                                                                       <  1:1632006/150‑1 (MQ=255)

                                  gggtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGtt                                                                                                                     >  1:5392612/1‑150 (MQ=255)

                                   ggtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGttt                                                                                                                    <  2:4750005/150‑1 (MQ=255)

                                   ggtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGttt                                                                                                                    >  1:4991091/1‑150 (MQ=255)

                                    gtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGtttt                                                                                                                   <  2:1350262/150‑1 (MQ=255)

                                    gtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGtttt                                                                                                                   <  2:862199/150‑1 (MQ=255)

                                    gtagCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGtttt                                                                                                                   >  1:4688403/1‑150 (MQ=255)

                                       gCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGAt                                                                                                                <  2:4512463/150‑1 (MQ=255)

                                       gCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCGTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGAt                                                                                                                <  1:4813985/150‑1 (MQ=255)

                                        cGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATg                                                                                                               <  2:2904224/150‑1 (MQ=255)

                                          ttCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGcc                                                                                                             >  2:5004303/1‑150 (MQ=255)

                                          ttCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGcc                                                                                                             >  2:1062090/1‑150 (MQ=255)

                                           tCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCg                                                                                                            >  1:3278605/1‑150 (MQ=255)

                                            cccagcccagGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGt                                                                                                           >  1:1402852/1‑150 (MQ=255)

                                             ccagcccagGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTa                                                                                                          <  2:73759/150‑1 (MQ=255)

                                              cagcccagGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAg                                                                                                         >  2:5546527/1‑150 (MQ=255)

                                              cagcccagGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAg                                                                                                         <  1:5399361/150‑1 (MQ=255)

                                              cagcccagGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAg                                                                                                         >  1:4451669/1‑150 (MQ=255)

                                              cagcccagGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAg                                                                                                         <  2:1444365/150‑1 (MQ=255)

                                              cagcccagGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAg                                                                                                         >  2:1943003/1‑150 (MQ=255)

                                                  ccagGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCt                                                                                                     <  1:883118/150‑1 (MQ=255)

                                                       ccTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCg                                                                                                >  2:4095570/1‑150 (MQ=255)

                                                         ttGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGc                                                                                              <  1:4198362/150‑1 (MQ=255)

                                                                aTGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCACTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAggg                                                                                       >  1:5189607/1‑150 (MQ=255)

                                                                           cTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTac                                                                            >  1:1370121/1‑150 (MQ=255)

                                                                                aGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCt                                                                              <  2:1899359/143‑1 (MQ=255)

                                                                                    gCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCgg                                                                   <  1:2189686/150‑1 (MQ=255)

                                                                                     cccTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGc                                                                  >  1:5568329/1‑150 (MQ=255)

                                                                                     cccTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGc                                                                  <  1:1793482/150‑1 (MQ=255)

                                                                                       ctgctgGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCtt                                                                <  1:5075767/150‑1 (MQ=255)

                                                                                       ctgctgGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCtt                                                                >  2:3979310/1‑150 (MQ=255)

                                                                                               tAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCgg                                                        >  2:486006/1‑150 (MQ=255)

                                                                                               tAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCgg                                                        <  1:3664428/150‑1 (MQ=255)

                                                                                                 aTGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCg                                                      >  1:830714/1‑150 (MQ=255)

                                                                                                      ggTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATa                                                 <  2:531993/150‑1 (MQ=255)

                                                                                                            tttCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTg                                           >  1:1925606/1‑150 (MQ=255)

                                                                                                            tttCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTg                                           >  2:2553125/1‑150 (MQ=255)

                                                                                                            tttCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTg                                           >  2:4394260/1‑150 (MQ=255)

                                                                                                            tttCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTg                                           <  2:5425970/150‑1 (MQ=255)

                                                                                                            tttCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTg                                           >  2:3818974/1‑150 (MQ=255)

                                                                                                             ttCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTgg                                          <  1:4940296/150‑1 (MQ=255)

                                                                                                             ttCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTgg                                          <  2:3192457/150‑1 (MQ=255)

                                                                                                             ttCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTgg                                          <  2:2930086/150‑1 (MQ=255)

                                                                                                             ttCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTgg                                          <  1:3356556/150‑1 (MQ=255)

                                                                                                               cGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGc                                        >  1:5403778/1‑150 (MQ=255)

                                                                                                                       ccGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCTTACCCTTGGGCAGAAAGcc                                >  2:3903906/1‑150 (MQ=255)

                                                                                                                       ccGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGcc                                >  2:4665428/1‑150 (MQ=255)

                                                                                                                            ggCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTgg                           <  2:1566413/150‑1 (MQ=255)

                                                                                                                            ggCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTgg                           <  1:1252873/150‑1 (MQ=255)

                                                                                                                            ggCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTgg                           <  1:155566/150‑1 (MQ=255)

                                                                                                                             gCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTggg                          >  2:3934637/1‑150 (MQ=255)

                                                                                                                              cgctatgtGCGGTTGCACATGCAGCCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCACTGGgc                         >  2:268321/8‑150 (MQ=255)

                                                                                                                              cAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGgc                         <  1:3782623/150‑1 (MQ=255)

                                                                                                                                    cTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGc                   <  1:2323783/150‑1 (MQ=255)

                                                                                                                                     tGCGGTTGCACATGCAGCCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGcc                  >  1:2152455/1‑150 (MQ=255)

                                                                                                                                     tGCGGTTGCACATGCAGCCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCACTGGGCGCTGGcc                  <  1:268321/150‑1 (MQ=255)

                                                                                                                                      gCGGTTGCACATGCAGCCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCg                 >  1:4558980/1‑150 (MQ=255)

                                                                                                                                        tgcgctaTGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGa               >  2:2727254/7‑150 (MQ=255)

                                                                                                                                        ggCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGa               <  2:3057257/150‑1 (MQ=255)

                                                                                                                                              aTGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGa         <  2:3278605/150‑1 (MQ=255)

                                                                                                                                               tgtgCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATACAGCCCCTGGGCGCTGGCCGGTACCTGaa        <  1:644665/150‑1 (MQ=255)

                                                                                                                                               tgtgCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGaa        <  2:1557668/150‑1 (MQ=255)

                                                                                                                                               tgtgCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGaa        <  1:3930666/150‑1 (MQ=255)

                                                                                                                                               tgtgCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGaa        <  2:5426669/150‑1 (MQ=255)

                                                                                                                                                 tgCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAagg      <  1:5546527/150‑1 (MQ=255)

                                                                                                                                                  gCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAaggt     <  1:1136909/150‑1 (MQ=255)

                                                                                                                                                   cAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCTCCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGATCCTGAaggta    <  2:5189607/150‑1 (MQ=255)

                                                                                                                                                   cAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGaa        >  1:503650/1‑146 (MQ=255)

                                                                                                                                                   cAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAaggta    <  2:2365473/150‑1 (MQ=255)

                                                                                                                                                    aaCCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAaggtag   >  1:922180/1‑150 (MQ=255)

                                                                                                                                                    aaCCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAaggtag   <  1:4859077/150‑1 (MQ=255)

                                                                                                                                                     aCCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAaggtagg  >  2:3767129/1‑150 (MQ=255)

GCCATCGATCAGTGTGCTCAGCGTGACCGGGTAGGGGTAGCGTTCCCAGCCCAGGCCTTGCGTCATGGCGCCAACCTGGAAGTAGCCCTGCTGGGTAATGGCGGTGGCTTTCGCCAGGGCCGGTGGCAGGCCCTGCGGCTGCATGTGCAACCGCACATAGCGCAGCAGGTCGCTGGCGCTGGTTTTGATGCCGTAGGCCTCATCGGCGTAAGGGCCGGGCCCTACACGCACCGGCTTGTCGTCGGCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAAGGTAGG  >  NC_002947/3282195‑3282493

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