BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,501	G→A	5.2%	A17A (GCC→GCT)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,501	0	G	A	5.2%	410.2 / 9.9	175	A17A (GCC→GCT)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base G (76/88); new base A (5/4); total (82/93)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 175 total aligned reads displayed.

CTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATT  >  CP000730/36363‑36637

cttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                                        <  2:378390/141‑1 (MQ=35)

cttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                                        <  2:536554/141‑1 (MQ=35)

cttcCAAAGTTGATACAAAAGTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                                        <  1:277769/141‑1 (MQ=25)

            aTACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATGTCTCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAGCGGCTACAGTGATaaa                                                                                                                            >  1:26528/1‑140 (MQ=255)

                  aaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCtt                                                                                                                      <  1:883063/141‑1 (MQ=32)

                    ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                                                    <  2:714914/141‑1 (MQ=32)

                     tGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCCCGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                                                   <  2:288636/141‑1 (MQ=17)

                      ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATCTTTCAGATATATCACGATCACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTCCAGTGATAACATCCTTGTTa                                                                                                                  >  2:48257/1‑141 (MQ=12)

                      ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATGACATCCTTGTTa                                                                                                                  >  2:134463/1‑141 (MQ=17)

                      ggggCATATTCTTGAACCCAACGGTAGACCGCTGAATGATGAAAGTTTACGCAACGCAACTTTAATACTTTTGTTATATGCTGATAAATTG‑ATCTAAAGCTTAAGGGGTAGCAAACTGCTACAGTTTTAACATCCTTGTTa                                                                                                                  >  2:785553/1‑141 (MQ=2)

                         gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                                               >  2:653764/1‑141 (MQ=32)

                         gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                                               >  2:880897/1‑141 (MQ=32)

                           atatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGCGATAACATCCTTGTTAAATTg                                                                                                             <  1:286313/141‑1 (MQ=17)

                            tatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGt                                                                                                            <  2:27748/141‑1 (MQ=32)

                                 ttGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTtata                                                                                                       >  2:571846/1‑141 (MQ=32)

                                  tGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATCTTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTatat                                                                                                      >  1:153761/1‑141 (MQ=17)

                                   gAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTGTATc                                                                                                     >  2:894903/1‑141 (MQ=17)

                                    aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                                                    <  1:512041/141‑1 (MQ=32)

                                    aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                                                    <  2:879713/141‑1 (MQ=32)

                                            ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                                            <  1:530140/141‑1 (MQ=32)

                                            ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                                            <  1:462828/141‑1 (MQ=32)

                                            ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                                            >  2:291693/1‑141 (MQ=32)

                                            ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAAATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAAAGGCTACAGTGATAACATACTTGTTAAATTGTTTATATCTGAAATAGt                                                                                            >  2:224637/1‑141 (MQ=12)

                                              tAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTc                                                                                          >  1:806273/1‑141 (MQ=32)

                                               aGACCGTTGAATGAAGAACGTTTACACCACGTTCCCTTAATATTTCAGATATAACACGGTAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAGTGTTTATATCTGAAATAGTTCa                                                                                         <  1:659888/141‑1 (MQ=255)

                                                gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                                                        >  1:755353/1‑141 (MQ=32)

                                                  ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                                                      >  2:65562/1‑141 (MQ=32)

                                                  ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGATACAGTCATACCATCCTTGTTAAATTTGTTATATTTTAAATAGTTCATAc                                                                                      >  2:267128/1‑141 (MQ=11)

                                                      tGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaa                                                                                  >  2:100273/1‑141 (MQ=32)

                                                               taCGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                                                         <  2:411153/140‑1 (MQ=17)

                                                               aaCTTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                                                         <  2:80749/141‑1 (MQ=17)

                                                               aaCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                                                         <  2:819449/141‑1 (MQ=32)

                                                                 cGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGATGTTGTTCAGACAGAAGACTCCttttt                                                                       <  2:669330/141‑1 (MQ=255)

                                                                 cGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttttt                                                                       <  2:73262/141‑1 (MQ=32)

                                                                 cGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATAATTTAGATAGTAGCCAACGGCTACAGTGATAACATCATTGTTAAATTGTTTATATCTGAAATAGTTCATACAAAAGCATCCttttt                                                                       >  2:124551/1‑141 (MQ=11)

                                                                  gTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTg                                                                      >  1:18119/1‑141 (MQ=32)

                                                                     tACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACAACCAAGTTAAATTGTTTATATCTGAACGAGTTCGTACAGAAAACTCCTTTTTGTTa                                                                   >  2:150523/1‑141 (MQ=11)

                                                                      acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                                                  >  1:858728/1‑141 (MQ=32)

                                                                            cgttatcTTAATATTACGGATATATCACGAAAGATCAGTTGTA‑AACTTAGATTGTGGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTAAAATtat                                                            >  1:452373/7‑141 (MQ=11)

                                                                                 ccTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtat                                                       >  1:588038/1‑141 (MQ=32)

                                                                                 ccTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtat                                                       >  1:377630/1‑141 (MQ=32)

                                                                                   ttAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAGGACTCCTTTTTGTTAACATTATACTCTaa                                                     >  1:74860/1‑141 (MQ=12)

                                                                                      atatTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAtt                                                  <  2:367354/141‑1 (MQ=35)

                                                                                      atatTTCAGATATAACACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAtt                                                  <  2:120770/141‑1 (MQ=25)

                                                                                       aatTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAGTATACTATAAATTc                                                 <  2:641276/140‑1 (MQ=25)

                                                                                           actgATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGTAGACTCCTGTTTGTTAAAATTATACTATAAATTCAACt                                             <  1:187929/138‑1 (MQ=255)

                                                                                              gATATATCACGATAGCTCA‑ATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAATTTCAACTTTg                                          <  1:683885/141‑1 (MQ=255)

                                                                                               atatatCACGATAGCTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGCGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGc                                         <  1:587874/141‑1 (MQ=255)

                                                                                               atatatCACGATAGCTCA‑ATGGATCTCATGGAGAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGc                                         >  1:456209/1‑141 (MQ=11)

                                                                                                 atatCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATATCATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCaa                                       <  1:538986/141‑1 (MQ=25)

                                                                                                  tatCACGATAGCTCA‑ATGTATATATTAGTTAGTAGCCAACGGCTACGGTGATAACATCCTTGTTAAATTGTTGATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCGACTTTGCAAc                                      <  2:184198/141‑1 (MQ=11)

                                                                                                  tatCACGATAACTCA‑ATGTAAATTGTAGCTAGTATCCAACGGCTACAGTGATAACATCCTCGTTAAATTGTTTATACCTGAAATATTTCATACAGAAGACTCCTTTTTGTTAAAATTCTACTATACATTCAACTTTGCaca                                      >  2:435309/1‑139 (MQ=11)

                                                                                                     cACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCATTTTTGTTAAAATTATACTATAAATTCAACTTTCCacccga                                   >  1:116292/1‑136 (MQ=255)

                                                                                                          tagcTCA‑ATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                              <  1:868682/138‑1 (MQ=9)

                                                                                                          tagcTCA‑ATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                              <  2:727138/138‑1 (MQ=12)

                                                                                                          aagccta‑aTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                              <  2:324286/135‑1 (MQ=255)

                                                                                                           agcTCA‑ATTCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta                             <  2:598123/139‑1 (MQ=255)

                                                                                                            gcTCA‑ATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat                            <  2:770435/140‑1 (MQ=16)

                                                                                                            aCTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat                            <  1:817193/141‑1 (MQ=38)

                                                                                                             cTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatt                           >  2:578432/1‑141 (MQ=39)

                                                                                                             cTCA‑ATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt                           >  1:788610/1‑141 (MQ=12)

                                                                                                             aatcaacttataCCTCAGATAGTAGCCAACGCCTACAGCTATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTCCACACCGAAGACCCCTGTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat                            <  1:431876/133‑1 (MQ=11)

                                                                                                              tCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta                          <  2:502065/141‑1 (MQ=255)

                                                                                                              tCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta                          >  1:378210/1‑141 (MQ=40)

                                                                                                              tCA‑ATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTTTTTTTTTTTTATTAATTCAACTTTGCAACAGAACCGttttt                          >  2:611834/1‑140 (MQ=11)

                                                                                                               cA‑ATGTATAGCGTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtattat                         <  1:823997/141‑1 (MQ=255)

                                                                                                                a‑aTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATACAACTTTGCAACAGAACCGtatttta                        >  2:583813/1‑138 (MQ=255)

                                                                                                                 aattgatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATg                        <  2:589180/136‑1 (MQ=255)

                                                                                                                 aaaccatatCTTAAATATTAGCCAACGACTACAGTGATAACATCCTTGTTAAATTGTTTATATATGAAATAGTTCATACAGAATACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaact                        <  2:857539/136‑9 (MQ=11)

                                                                                                                  atgcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAACGtattttat                       >  1:372918/5‑137 (MQ=255)

                                                                                                                   tgcatatCTTATATATTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatg                      <  2:230690/138‑6 (MQ=11)

                                                                                                                     catatatTAGATAGTAGACAACGGCTACAGTGATAACATCCTCGTTAAATTGTTTACCTCTGAAATAGTTCATACAGACGACTATTTTTTGTTAAAATTATATTATAAATGCAACTTTGCATCAGAACCGtattttatgag                    <  2:818040/140‑8 (MQ=11)

                                                                                                                     catatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgag                    <  1:483229/140‑8 (MQ=255)

                                                                                                                     catatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgag                    <  1:745899/140‑8 (MQ=255)

                                                                                                                      atatCTTAGATAGTAGCCAACGGCTACAGTTATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATa                   <  1:124551/141‑1 (MQ=255)

                                                                                                                      atatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTTAAATTATATTATAAATTCAACTTTGCAACAGAAcatcctaactgtatt                   >  2:586003/1‑127 (MQ=11)

                                                                                                                      atatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtatt                   <  2:309778/141‑14 (MQ=11)

                                                                                                                       tatCTTAGCTGTTATCCAACAGCTACAGTGATTATTCCCGTGTTCAATTGTTTCTTGCACAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactatatg                  >  2:429497/1‑127 (MQ=2)

                                                                                                                       tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATTTGAAATAGTGCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagta                  <  1:220023/141‑10 (MQ=11)

                                                                                                                       tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCACCAGAAcctcctaactgtatta                  >  2:828973/1‑127 (MQ=11)

                                                                                                                       tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtatta                  >  2:162247/1‑127 (MQ=11)

                                                                                                                       tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagta                  <  1:235621/141‑10 (MQ=255)

                                                                                                                       tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAAAGAAccacctaactttatta                  >  2:496240/1‑127 (MQ=11)

                                                                                                                       tatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactatatg                  >  1:298913/1‑127 (MQ=11)

                                                                                                                       tatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactatatg                  <  1:774077/141‑15 (MQ=11)

                                                                                                                       tatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactatatg                  >  1:299376/1‑127 (MQ=11)

                                                                                                                         tCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtaca                >  1:929303/1‑130 (MQ=11)

                                                                                                                         tCTTACATAGTTGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACATAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtaca                <  2:771740/141‑12 (MQ=11)

                                                                                                                           ttAGATAGTAGCCAACGGCTAGAGTGATAACATCCTTGTTAATTTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacact              <  2:4924/141‑14 (MQ=11)

                                                                                                                           ttAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacact              >  1:300190/1‑128 (MQ=11)

                                                                                                                           ttAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacact              >  1:792446/1‑128 (MQ=11)

                                                                                                                            tAGATAGTAGCCAACGGCTAGAGTGATAACATCCTTGTTAAATGGTTTATATCTGAAATAGTTCATACAGAAGGCTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacactc             <  2:343271/141‑15 (MQ=11)

                                                                                                                            tAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacactc             <  2:924246/141‑15 (MQ=11)

                                                                                                                            tAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACTGTATTATGGAATAGAGATg             <  1:519723/141‑1 (MQ=255)

                                                                                                                            tAGATAGTAGCCAACGGCGACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGCACCGtattttatgagtacacgc             <  1:132767/141‑15 (MQ=11)

                                                                                                                              gatatgcgCCAACGGTCACAGTGATAACATCCTGGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGCCTCCTTTTTGTTAACCTTATACACCAAATCCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGtt           <  2:387371/134‑1 (MQ=17)

                                                                                                                                   tagCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTaa      >  2:761238/1‑141 (MQ=255)

                                                                                                                                   tagCCAACCGCTACAGTTATAACATCTTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTaa      <  1:417290/141‑1 (MQ=255)

                                                                                                                                    agCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAAc     >  2:753232/1‑141 (MQ=255)

                                                                                                                                    agCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAAc     >  1:331469/1‑141 (MQ=255)

CTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATT  >  CP000730/36363‑36637

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