BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,120,685	G→A	6.2%	S30S (AGC→AGT)	USA300HOU_RS10785 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,120,685	0	G	A	6.2%	399.6 / 10.9	128	S30S (AGC→AGT)	USA300HOU_RS10785	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (69/51); new base A (5/3); total (74/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

Only 100 of 134 total aligned reads displayed.

TCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCATTTT  >  USA300TCH1516_ALE/2120545‑2120825

tCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACg                                                                                                                                               >  3:273068/1‑141 (MQ=255)

 ccTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGc                                                                                                                                              >  2:760113/1‑141 (MQ=255)

     ttCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCtata                                                                                                                                          >  2:538008/1‑141 (MQ=255)

     ttCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCtata                                                                                                                                          >  1:644608/1‑141 (MQ=255)

      tCATCTGTGTTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACtatat                                                                                                                                         >  1:526314/1‑140 (MQ=255)

      tCATCTGTGTTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACtatat                                                                                                                                         >  1:734959/1‑140 (MQ=255)

      tCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATAc                                                                                                                                         <  2:1003342/141‑1 (MQ=255)

          cTGTGTTCTTTTCCATCCCTTCATCATTACTTTGTTCCAATCTATCTAAATCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTc                                                                                                                                     <  3:199071/141‑1 (MQ=255)

          cTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTc                                                                                                                                     <  1:429292/141‑1 (MQ=255)

          cTGTGTTATTTTCAATCCAATCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTc                                                                                                                                     <  3:179267/141‑1 (MQ=255)

            gtgtTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCgc                                                                                                                                   <  1:455632/141‑2 (MQ=255)

            gtgtTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGt                                                                                                                                   >  2:635553/1‑141 (MQ=255)

               ttATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTc                                                                                                                                <  2:629719/141‑1 (MQ=255)

                 aTTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCtt                                                                                                                              <  1:299598/141‑1 (MQ=255)

                      cAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAActc                                                                                                                         >  1:372840/1‑141 (MQ=255)

                      cAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAActc                                                                                                                         <  1:309741/141‑1 (MQ=255)

                       aaTCCATTCATACACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAActct                                                                                                                        <  3:321195/141‑1 (MQ=255)

                       aaTCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCGCTTCTTAActct                                                                                                                        <  2:82870/141‑1 (MQ=255)

                         tCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTa                                                                                                                      >  4:131355/1‑141 (MQ=255)

                         tCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTa                                                                                                                      >  1:992887/1‑141 (MQ=255)

                             ttCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTg                                                                                                                  <  1:641314/141‑1 (MQ=255)

                              tCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑a                                                                                                                >  4:147125/1‑141 (MQ=255)

                                 tAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑aaaa                                                                                                             >  2:1046637/1‑141 (MQ=255)

                                     cacTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGGTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCaa                                                                                                         >  1:432525/1‑141 (MQ=255)

                                      acTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCaaa                                                                                                        >  2:188678/1‑141 (MQ=255)

                                         ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                                                                     >  1:940092/1‑141 (MQ=255)

                                         ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                                                                     >  2:944847/1‑141 (MQ=255)

                                         ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                                                                     >  1:1050983/1‑141 (MQ=255)

                                         ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                                                                     >  2:1054135/1‑141 (MQ=255)

                                           gTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCAt                                                                                                   <  1:929619/141‑1 (MQ=255)

                                           gTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCAt                                                                                                   >  3:271611/1‑141 (MQ=255)

                                                 aaaTATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAAATCAAAGCCATCTACa                                                                                              >  1:897473/1‑141 (MQ=255)

                                                   atatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACatta                                                                                           <  1:460266/141‑1 (MQ=255)

                                                   atatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACatta                                                                                           <  2:669992/141‑1 (MQ=255)

                                                    tatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaa                                                                                          <  1:797050/141‑1 (MQ=255)

                                                      tCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaatt                                                                                        <  1:631767/141‑1 (MQ=255)

                                                       cTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaatta                                                                                       <  4:228471/141‑1 (MQ=255)

                                                         aaaGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAattatt                                                                                     <  1:522513/141‑1 (MQ=255)

                                                          aaGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTAttt                                                                                    >  2:1058932/1‑141 (MQ=255)

                                                               gtgtATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGa                                                                               >  2:409349/1‑141 (MQ=255)

                                                                  tataGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACAt                                                                            <  1:439432/141‑1 (MQ=255)

                                                                      gATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGt                                                                        <  4:264319/141‑1 (MQ=255)

                                                                          gTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAcc                                                                    <  1:617052/141‑1 (MQ=255)

                                                                            tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCTTTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCACAGTCTTTTACATTAATTATTTTTCGACATACGGCACCCg                                                                  <  4:174304/141‑1 (MQ=255)

                                                                            tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                                                                  >  2:625438/1‑141 (MQ=255)

                                                                            tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                                                                  >  2:883004/1‑141 (MQ=255)

                                                                            tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                                                                  >  2:927714/1‑141 (MQ=255)

                                                                             ctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGt                                                                 >  2:723508/1‑141 (MQ=255)

                                                                             ctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGt                                                                 >  2:723512/1‑141 (MQ=255)

                                                                              tcgtcaGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCGCTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAAATATTTTTCGACATACGTCACCCGtt                                                                >  1:22892/4‑141 (MQ=255)

                                                                               catcaGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGttt                                                               >  4:421057/1‑141 (MQ=255)

                                                                                   aGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCTTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTa                                                           >  1:340122/1‑141 (MQ=255)

                                                                                      aaCTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATa                                                        <  4:248214/141‑1 (MQ=255)

                                                                                         tagtagCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCGCTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAAATATTTTTCGACATACGTCACCCGTTTTGTAATACAt                                                     >  1:715115/1‑141 (MQ=255)

                                                                                         tAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAcc                                                                    <  1:1074934/126‑1 (MQ=255)

                                                                                         tAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAcc                                                                    >  2:1074934/1‑126 (MQ=255)

                                                                                          aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt                                                    >  1:1014721/1‑141 (MQ=255)

                                                                                          aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt                                                    >  3:57508/1‑141 (MQ=255)

                                                                                          aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt                                                    >  2:993461/1‑141 (MQ=255)

                                                                                          aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt                                                    <  1:923004/141‑1 (MQ=255)

                                                                                           gTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTg                                                   >  2:176986/1‑141 (MQ=255)

                                                                                             aaCCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGtt                                                 >  1:26639/1‑141 (MQ=255)

                                                                                                  gTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGt                                            >  1:572862/1‑141 (MQ=255)

                                                                                                  gTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGt                                            >  1:789354/1‑141 (MQ=255)

                                                                                                    cAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGtcc                                          >  2:400888/1‑141 (MQ=255)

                                                                                                      aTTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGtcctc                                        <  3:403180/141‑1 (MQ=255)

                                                                                                      aTTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATAAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGtcctc                                        <  1:590673/141‑1 (MQ=255)

                                                                                                        tGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTcctcct                                      <  2:205492/141‑1 (MQ=255)

                                                                                                         gAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCtt                                     >  2:114196/1‑141 (MQ=255)

                                                                                                            tttCCTCCCTAGTCTAGAACAACCATTTCGACGCAATACGCTCGCTTCTTAACTCTTAATTG‑AAAAGCAAAGCCATGTACCTTAGTTCTTTTTCGCCATTCGTCACCCGTTTCGTAACGCATTGTTCTCGTCCTCCTTGTc                                  <  2:777534/141‑1 (MQ=255)

                                                                                                             ttCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGGCATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCa                                 <  1:842338/141‑1 (MQ=255)

                                                                                                              tCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCCTCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTATTACATTGTTCTCGTCCTCCTTGTCAt                                >  1:217886/1‑141 (MQ=255)

                                                                                                               ccTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTGATACATTGTTCTCGTCCTCCTTGTCAtt                               >  1:755978/1‑141 (MQ=255)

                                                                                                               ccTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCAtt                               >  1:996421/1‑141 (MQ=255)

                                                                                                                   cATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AACATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTAGTGTAGTACAGTGTTCCTGACGTCTTTGTCATTATCt                           >  2:271327/1‑141 (MQ=255)

                                                                                                                   cATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCt                           >  4:290744/1‑141 (MQ=255)

                                                                                                                    aTGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑CAAATCAAAGCCATCTACATTAATTATTTTTCGACATACTTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCta                          <  2:754503/141‑1 (MQ=255)

                                                                                                                     tGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCtat                         >  1:861688/1‑141 (MQ=255)

                                                                                                                     tGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCtat                         >  4:229346/1‑141 (MQ=255)

                                                                                                                         cTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCATCCTTGTCATTATCTATAGCt                     <  3:238133/141‑1 (MQ=255)

                                                                                                                           aGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAg                   <  1:249101/141‑1 (MQ=255)

                                                                                                                           aGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAg                   <  1:700825/141‑1 (MQ=255)

                                                                                                                               cgactaTTTCGACACTATATGCTCGCTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAAATATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTTTAGTCCTCCTTGTCGTCATCTATACCGAGAAtt               <  1:705909/136‑1 (MQ=255)

                                                                                                                               caacCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAAtt               >  1:21411/1‑141 (MQ=255)

                                                                                                                               caacCATNTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAAtt               <  2:816252/141‑1 (MQ=255)

                                                                                                                                  ccATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAAttttt            >  3:166615/1‑141 (MQ=255)

                                                                                                                                     tttCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGTCTTACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAAc         >  2:38862/1‑141 (MQ=255)

                                                                                                                                     tttCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAAc         >  1:66218/1‑141 (MQ=255)

                                                                                                                                      ttCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACg        <  1:189671/141‑1 (MQ=255)

                                                                                                                                       tCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGc       >  2:540546/1‑141 (MQ=255)

                                                                                                                                        cGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCa      >  2:600363/1‑141 (MQ=255)

                                                                                                                                        cGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCa      >  4:347669/1‑141 (MQ=255)

                                                                                                                                        cGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCa      >  2:366528/1‑141 (MQ=255)

                                                                                                                                        cGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTGCGACATACGTCACCCGTGTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCa      <  3:444248/141‑1 (MQ=255)

                                                                                                                                         gACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTcctcct                                      <  4:236595/108‑1 (MQ=255)

                                                                                                                                         gACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCAt     >  3:345848/1‑141 (MQ=255)

                                                                                                                                         gACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCAt     <  1:671547/141‑1 (MQ=255)

                                                                                                                                           cGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCAttt   >  3:9755/1‑141 (MQ=255)

                                                                                                                                            gCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTTTTCATACTTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCAtttt  >  2:316756/1‑141 (MQ=255)

                                                                                                                                            gCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCAtttt  >  1:949581/1‑141 (MQ=255)

TCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCATTTT  >  USA300TCH1516_ALE/2120545‑2120825

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