BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	1,006,124	A→T	S169R (AGT→AGA)	USA300HOU_RS04965 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,006,124	0	A	T	99.4%	564.9 / NA	173	S169R (AGT→AGA)	USA300HOU_RS04965	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base T (110/62); minor base G (0/1); total (110/63)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

AAGTTGTTGTTGCTGTATATAAGTTTCTAATTTAGGAGGAAGATTCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCACTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCTAAATCAATACCTGCATTTTGAAT  >  USA300TCH1516_ALE/1005980‑1006261

aaGTTGTTGTTGCTGTATATAAGTTTCTAATTTAGGAGGAAGATTCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTaa                                                                                                                                       <  2:767021/149‑1 (MQ=255)

   ttgttgttgCTGTATATAAGTTTCTAATTTAGGAGGAAGATTCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTGTAGATGTGATGTAGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGa                                                                                                                                    >  2:1246418/1‑149 (MQ=255)

     tttgttgCTGTATATAAGTTTCTAATTTAGGAGGAAGATTCTTTTCGATAGGTTTACCCATCAGTTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGaaa                                                                                                                                  <  2:1176046/148‑1 (MQ=255)

         ttgCTGTATATAAGTTTCTAATTTAGGAGGAAGATTCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACga                                                                                                                              <  2:1685286/149‑1 (MQ=255)

             tGTATATAAGTTTCTAATTTAGGAGGAAGATTCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACgagaga                                                                                                                          >  1:1263408/1‑149 (MQ=255)

               tatataAGTTTCTAATTTAGGAGGAAGATTCCTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAg                                                                                                                        <  1:1581936/149‑1 (MQ=255)

               tatataAGTGTCTAATTTAGAAGGAAGATTCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATGATCCTTAGACTTTTAAACCTTCTTCATGACCATCTTCTATCAAGGAAACgaaggaag                                                                                                                        >  2:561171/1‑149 (MQ=255)

                 tataAGTTTCTAATTTAGGAGGAAGATTCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAggg                                                                                                                      >  1:618935/1‑149 (MQ=255)

                 tataAGTTTCTAATTTAGGAGGAAGATTCTTTTCGAGAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAggg                                                                                                                      >  1:1092886/1‑149 (MQ=255)

                    aaGTTTCTAATTTAGGAGGAAGATTCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAACCGAGAGAAGGGGCt                                                                                                                   >  1:806783/1‑149 (MQ=255)

                    aaGTTTCTAATTTAGGAGGAAGATTCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCt                                                                                                                   >  2:1706749/1‑149 (MQ=255)

                    aaGTTTCTAATTTAGGAGGAAGATTCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAAAGAGAGAAGGGGCt                                                                                                                   >  2:157815/1‑149 (MQ=255)

                            aaTTTAGGAGGAAGATTCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGAt                                                                                                           >  2:929589/1‑149 (MQ=255)

                                tAGGAGGAAGATTCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTcc                                                                                                       >  1:172210/1‑149 (MQ=255)

                                  ggaggaAGATTCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCAt                                                                                                     >  2:1017219/1‑149 (MQ=255)

                                  ggaggaAGATTCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCAt                                                                                                     >  1:1462450/1‑149 (MQ=255)

                                          cttCTTTTCTATAGGTNTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATTGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCGCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTg                                                                                             <  1:1493850/148‑1 (MQ=255)

                                          aTTCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTg                                                                                             <  1:605213/149‑1 (MQ=255)

                                           ttCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGc                                                                                            <  2:1676018/149‑1 (MQ=255)

                                           ttCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCAGGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGc                                                                                            <  2:892893/149‑1 (MQ=255)

                                            tCTTTTCGATAGGTTTAACCATCCCTTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCa                                                                                           <  2:1122475/149‑1 (MQ=255)

                                                    aTAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGTTTTCCAGTTCTCTTGCAATATggca                                                                                   >  1:212240/1‑146 (MQ=255)

                                                    aTAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGcca                                                                                   <  2:1462450/149‑3 (MQ=255)

                                                    aTAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCaa                                                                                   <  2:654025/149‑1 (MQ=255)

                                                    aTAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGCGGCTTGTTCGATTTCCATTTCTCTTGCAACATGCaa                                                                                   >  1:44234/1‑149 (MQ=255)

                                                    aTAGGTTTACACAGCAATTCATTAATTATATAAGTATAGTTGAGCCACAGCTATTAACGTTCGTCTTATAGTCCTTGCTCATTAACATCTAATCTAAACAGCACGTGGGAAGGGGCTTGTTCGATTTCCATTTCTCATGCATTATGCaa                                                                                   <  2:226189/149‑1 (MQ=255)

                                                         tttACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATcca                                                                              >  1:177830/1‑147 (MQ=255)

                                                         tttACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCaa                                                                              >  1:240472/1‑149 (MQ=255)

                                                                     ttCATTAATTATATAAGTATAGATGTGATATGGGTATAATACTTCAACGTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTACAATATGAAAATCAATTTTCAAGCttt                                                                  >  2:32481/1‑149 (MQ=255)

                                                                       cATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCg                                                                <  2:703990/149‑1 (MQ=255)

                                                                            attaTATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCACGCTTTCGGTTAg                                                           <  1:1170704/149‑1 (MQ=255)

                                                                             ttaTATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGt                                                          <  1:361503/149‑1 (MQ=255)

                                                                             ttaTATAAGGATAGATGTGATATGGGAATAATACTTCGACTATTAAACCTTCATCATGAACATCTTCTCTCAAGAAAACGAGAGAAGGGGCATGTTCGATTTCCATTAATCTTGCACGGTGCAAATATATTTTCAAGATTTCGGATatt                                                          >  2:429185/1‑147 (MQ=255)

                                                                                tataAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTa                                                       <  2:851970/149‑1 (MQ=255)

                                                                                tataAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTa                                                       <  2:1721285/149‑1 (MQ=255)

                                                                                tataAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTa                                                       <  2:1227322/149‑1 (MQ=255)

                                                                                  taAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGACAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACt                                                     >  2:1559976/1‑149 (MQ=255)

                                                                                  taAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACt                                                     >  1:523225/1‑149 (MQ=255)

                                                                                  taAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACt                                                     >  2:1135443/1‑149 (MQ=255)

                                                                                  taAGTATAGATGTAATATAGGTACGATGCTTGGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTGGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACt                                                     >  2:870556/1‑149 (MQ=255)

                                                                                   aaGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGCTTAGTTTACtt                                                    >  2:173583/1‑149 (MQ=255)

                                                                                    aGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACttt                                                   >  1:790229/1‑149 (MQ=255)

                                                                                      tcacactgtgCTATCGGTATAATCCTTCGACGTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACttttt                                                 <  2:333258/143‑1 (MQ=255)

                                                                                      tataGATGTTATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACttttt                                                 <  2:562132/149‑1 (MQ=255)

                                                                                      tataGATGTGATATGGGTATAATCCTTCGACTTTTACACCTTCTTCCTGCCAATCTTCTCTTAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACttttt                                                 <  2:1648252/149‑1 (MQ=255)

                                                                                      tataGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAACCCGATTGTCAAGCTTTCGGTTAGTTTACttttt                                                 <  2:1223266/149‑1 (MQ=255)

                                                                                      tataGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACttttt                                                 >  1:1069828/1‑149 (MQ=255)

                                                                                      tataGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAAAGAGGGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATAAATTTTCAAGCTTTCGGTTAGGTTACttttt                                                 >  2:314079/1‑149 (MQ=255)

                                                                                         aGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGt                                              >  2:867736/1‑149 (MQ=255)

                                                                                          gATGTGGTATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCATTCGGTTAGTTTGCTTTTTTGTa                                             >  1:6103/1‑149 (MQ=255)

                                                                                          gATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCGAATCAATTTTCAAGCTTTCGGTTAGTTTACttttttttt                                             >  2:404537/1‑146 (MQ=255)

                                                                                          gATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACttttttttt                                             >  2:1686756/1‑146 (MQ=255)

                                                                                          gATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTa                                             <  1:995934/149‑1 (MQ=255)

                                                                                          gATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTa                                             >  1:705328/1‑149 (MQ=255)

                                                                                           aTGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTaa                                            >  2:136325/1‑149 (MQ=255)

                                                                                           aTGTGATATGGGTATAATCCTTCGACTTTGAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGACGGGGCTTGTTCGATTTCCATTTCTCTTGCAATGTGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTaa                                            >  1:897593/1‑149 (MQ=255)

                                                                                             gtgATATGGGTATAATTCTTGGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAgt                                          >  2:203121/1‑149 (MQ=255)

                                                                                             gtgATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAgt                                          >  1:1716452/1‑149 (MQ=255)

                                                                                             gtgATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAgt                                          >  2:1262313/1‑149 (MQ=255)

                                                                                             gtgATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAgt                                          >  2:375634/1‑149 (MQ=255)

                                                                                               gATATGGGTATAATCCTTCGGCTTTTAATCCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAgtcg                                        <  1:1157388/149‑1 (MQ=255)

                                                                                               gATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAgtcg                                        >  1:1574382/1‑149 (MQ=255)

                                                                                               gATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAgtcg                                        <  1:877001/149‑1 (MQ=255)

                                                                                               gATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTATCTAAAGAAAACGAGGGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGAAAATCAATTTTCAAGATTTCGGTTAGTTTACTTTTTTGTaatgct                                        >  2:620877/1‑145 (MQ=255)

                                                                                                atatGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAgtcgt                                       <  1:23782/149‑1 (MQ=255)

                                                                                                atatGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAgtcgt                                       <  1:1476498/149‑1 (MQ=255)

                                                                                                 tatGGGTATAATCCTTCGACTTTTAAACCTTCTTCATTAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAgtcgtc                                      <  2:1390701/149‑1 (MQ=255)

                                                                                                 tatGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAgtcgtc                                      <  1:1419810/149‑1 (MQ=255)

                                                                                                 tatGCGTATAATCCTTCGGCGTGTAACCCGCCTTCATGAACATTTTCTCTAAAGAAACCGAGAGAAGGGGCTTGTTCTATTTCCATTTCTCCTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAgtcgtc                                      <  1:202393/149‑1 (MQ=255)

                                                                                                  atGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAACATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTAGTCt                                     >  1:317447/1‑149 (MQ=255)

                                                                                                  atGGGTATAATCCTTCGACTTTGAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCt                                     >  1:958214/1‑149 (MQ=255)

                                                                                                    gggTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCttt                                   <  2:1338910/149‑1 (MQ=255)

                                                                                                    gggTATAATCCCTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCttt                                   <  2:1424876/149‑1 (MQ=255)

                                                                                                     ggTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACttttt                                                 <  1:502847/134‑1 (MQ=255)

                                                                                                      gTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTaa                                 >  1:1177056/1‑149 (MQ=255)

                                                                                                       tataATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTACGTCGTCTTTaaa                                >  2:602115/1‑149 (MQ=255)

                                                                                                         taATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATa                              <  1:1606693/149‑1 (MQ=255)

                                                                                                           aTCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGCTTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTTTATTTTGTCTTtacttact                            >  2:1609258/1‑149 (MQ=255)

                                                                                                            tCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACtt                           >  1:225451/1‑149 (MQ=255)

                                                                                                              cTTCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCt                         <  2:1069560/149‑1 (MQ=255)

                                                                                                                tCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGGGAAGGGGCTTGTTCGATTTCCATTTCTCTTCCAATATGCAAATCAATTTTCAAGCTTTCGGTCAGTTTACTTTTTTGTACGTCGTCTTTAACTACTTCtca                       >  2:1295651/1‑147 (MQ=255)

                                                                                                                tCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCTaa                       >  1:957909/1‑149 (MQ=255)

                                                                                                                tCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCTaa                       >  2:1369984/1‑149 (MQ=255)

                                                                                                                tCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCTaa                       >  1:1332821/1‑149 (MQ=255)

                                                                                                                tCGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCTaa                       >  1:1437725/1‑149 (MQ=255)

                                                                                                                 cGACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCTaaa                      >  1:192572/1‑149 (MQ=255)

                                                                                                                  gACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCTAAAt                     >  1:89684/1‑149 (MQ=255)

                                                                                                                  gACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCTAAAt                     >  1:1646108/1‑149 (MQ=255)

                                                                                                                  gACTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAACCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTCATTCATCTTTAAATACTTCTAAAt                     >  2:1234605/1‑149 (MQ=255)

                                                                                                                   aCTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCTAAATc                    >  2:1598787/1‑149 (MQ=255)

                                                                                                                   aCTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCTAAATc                    >  2:1227705/1‑149 (MQ=255)

                                                                                                                   aCTTTTAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCAGTTCTCTTACAATATGCAAACCAATTTTCAAGCTTACGGTTAGTTTACTTTTTTTTAAGTCGTCTTTATATATGTGTAAATc                    >  2:1456513/1‑149 (MQ=255)

                                                                                                                      tttAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGGCGTCTTTAAATACTTCTAAATCatt                 >  2:500484/1‑147 (MQ=255)

                                                                                                                       ttAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAGGTCGTCTTTAAATACTTCTAAATCAATa                <  1:1051220/149‑1 (MQ=255)

                                                                                                                       ttAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCTAAATCAATa                <  2:1627622/149‑1 (MQ=255)

                                                                                                                       ttAAACCTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCTAAATCAATa                <  2:1590661/149‑1 (MQ=255)

                                                                                                                            ccTTCTTCATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCTAAATCAATACCTGc           >  1:995330/1‑149 (MQ=255)

                                                                                                                                   cATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTTTATTTCTTCTTTAAATACTTCTAAATCAATACTTGCATTTTGa    >  2:1622031/1‑149 (MQ=255)

                                                                                                                                   cATGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCTAAATCAATACCTGCATTTTGa    <  1:1369494/149‑1 (MQ=255)

                                                                                                                                    aTGAACATCTTCTCTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCTAAATCAATACCTGCATTTTGaa   >  1:549738/1‑149 (MQ=255)

AAGTTGTTGTTGCTGTATATAAGTTTCTAATTTAGGAGGAAGATTCTTTTCGATAGGTTTACCCATCAATTCATTAATTATATAAGTATAGATGTGATATGGGTATAATCCTTCGACTTTTAAACCTTCTTCATGAACATCTTCACTAAAGAAAACGAGAGAAGGGGCTTGTTCGATTTCCATTTCTCTTGCAATATGCAAATCAATTTTCAAGCTTTCGGTTAGTTTACTTTTTTGTAAGTCGTCTTTAAATACTTCTAAATCAATACCTGCATTTTGAAT  >  USA300TCH1516_ALE/1005980‑1006261

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