Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	882,089	G→A	8.0%	T11I (ACA→ATA)	USA300HOU_RS04300 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	882,089	0	G	A	8.0%	180.3 / 3.4	67	T11I (ACA→ATA)	USA300HOU_RS04300	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (38/22);  new base A (4/2);  total (42/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

AAACAAGTAAAGAATGATCGATGTTACTAATATTGCCATTCTCCAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACATTCATGTTCAACCAAAATTTTGTTCGAATTCAATAAATGTCTTGTTTAAAATA  >  USA300TCH1516_ALE/881952‑882227                                                                                                                                          |                                                                                                                                           aaaCAAGTAAAGAATGATCGATGTTACTAATAGTGCCATTCTCCAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTa                                                                                                                                          <  2:748492/140‑1 (MQ=255)    cAAGTAAAGAATGATCGATGTTACTAATATTGCCATTCTCCAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAgat                                                                                                                                       <  2:615003/140‑1 (MQ=255)       gTAAAGAATGATCGATGTTACTAATATTGCCATTCTCCAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAgcttga                                                                                                                                    >  2:728349/1‑140 (MQ=255)        tAAAGAATGATCGATGTTACTAATATTGCCATTCTCCAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAgattgat                                                                                                                                   <  2:675385/140‑1 (MQ=255)             aaTGATCGATGTTACTAATATTGCCATTCTCCAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATCGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGAttctta                                                                                                                              >  2:346540/1‑140 (MQ=255)             aaTGATCGATGTTACTAATATTGCCATTCTCCAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGAttatta                                                                                                                              >  1:228788/1‑140 (MQ=255)                gATCGATGTTACTAATATTGCCATTCTCCAAATTTATTTCTGTGAGTATTTGGACGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGCTGTATTAAAAAAGGTTTTAGAGGTTGATTCTTAata                                                                                                                           >  1:139052/1‑140 (MQ=255)                 aTCGATGTTACTAATATTGCCATTCTCCAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAatat                                                                                                                          <  1:84774/140‑1 (MQ=255)                   cGATGTTACTAATATTGCCATTCTCCAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGc                                                                                                                        >  2:581507/1‑140 (MQ=255)                           cTAATATTGCCATTCTCCAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACt                                                                                                                <  1:146650/140‑1 (MQ=255)                                    ccATTCTCCAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCatgtatgt                                                                                                       >  2:802023/1‑140 (MQ=255)                                     cATTCTCCAAATTTATTTCTGTGAGTATTTGGAAGCTACCAGTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGtt                                                                                                      <  1:410697/140‑1 (MQ=255)                                      aTTCTCCAAATTTATTTCTGTGAGTATTTGGAACCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGttt                                                                                                     >  2:206621/1‑140 (MQ=255)                                       ttCTCCAAATTTATTTCGGCGAGGATTTGGAAGCTCCCATTAGGCAACGGTTTAACAATCGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGACTGATTCTTAATATGCCATTCACTCAGGAATGcttt                                                                                                    >  2:99778/1‑140 (MQ=255)                                        tctcCAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGTCAATTGCTTTTCAGCTTGTTGTATTAAAAAAGGTTTTGTAGCTGGATTATTAATCTGCTATTCAATCATGTATGttttt                                                                                                   >  2:740853/1‑140 (MQ=255)                                        tctcCAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGttttt                                                                                                   >  1:201434/1‑140 (MQ=255)                                          tcCAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGTTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCa                                                                                                 >  1:719482/1‑140 (MQ=255)                                          tcCAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCACTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCa                                                                                                 >  1:127741/1‑140 (MQ=255)                                           ccAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTCTATTAAGTATTAATATGCCATTCACTCATGTATGTTTTTCAc                                                                                                <  1:308052/140‑1 (MQ=255)                                           ccAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCAc                                                                                                <  1:388455/140‑1 (MQ=255)                                            cAAATTTATTTCTGTGAGTATTTGGAAGCTACCCTTAGGCAACGGTTTAACAATAGACAATTGCTGTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCCTTCACTCGTGTATGTTTTCCACt                                                                                               >  2:790173/1‑140 (MQ=255)                                            cAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACt                                                                                               <  1:188826/140‑1 (MQ=255)                                              aaTTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTcc                                                                                             >  2:112332/1‑140 (MQ=255)                                                 ttatttCTGTAAGTATTTGAAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATACCACAATCGGGCGTTTATGTTTTTCACTCCTGc                                                                                          <  2:176060/140‑1 (MQ=255)                                                  tatttCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCt                                                                                         >  2:539865/1‑140 (MQ=255)                                                   atttCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCtt                                                                                        <  2:315647/140‑1 (MQ=255)                                                   atttCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCtt                                                                                        >  1:437421/1‑140 (MQ=255)                                                      tCTGTGAGTATTTGGAAGCTACCATTCGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTaa                                                                                     <  2:205216/140‑1 (MQ=255)                                                            aGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAACTAgg                                                                               >  1:442042/1‑140 (MQ=255)                                                            aGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAgg                                                                               >  2:108089/1‑140 (MQ=255)                                                             gTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTTCTTTAAAATAggg                                                                              >  2:41865/1‑140 (MQ=255)                                                             gTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAggg                                                                              >  1:726818/1‑140 (MQ=255)                                                               aTTTGGAAGCTACTATTAGGCAGCGGTTTAACAATAGAGAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTACTGTATATTTTTCACTTCTGCTTTAAAATAGGGtt                                                                            >  2:438364/1‑140 (MQ=255)                                                               aTTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGtt                                                                            >  1:663469/1‑140 (MQ=255)                                                               aTTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTATTCCGCGTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTCTGTTTTTCACTTCTGCTTTAAAATAGGGtt                                                                            >  2:742974/1‑140 (MQ=255)                                                                  tGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGa                                                                         >  2:670088/1‑140 (MQ=255)                                                                  tGGAAGCGACTATTGTGCATCCGTTTAACAATAGACCGTTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGa                                                                         <  2:222428/140‑1 (MQ=255)                                                                    gAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTATAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGaaa                                                                       >  1:622522/1‑140 (MQ=255)                                                                    gAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTATAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGaaa                                                                       >  1:27918/1‑140 (MQ=255)                                                                           ccATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTtata                                                                >  1:716760/1‑140 (MQ=255)                                                                             aTTATGCAACGTTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTCGAAAAGCATAGTTAGATGGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGt                                                              <  1:374618/140‑1 (MQ=255)                                                                             aTTAGGCAACGGTTTAACAATAGAGAAAGGCGTTCGAAGTTGTTGAACTAAAAAAGTTTTTATAGAGTGACTGTTAATATGCCATTCACTCAGTTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTCTAGt                                                              >  2:534686/1‑140 (MQ=255)                                                                                 gtcAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTgag                                                          <  2:39791/138‑1 (MQ=255)                                                                                    aaCGGTTTAACAATAGACAATTGCTTTTCTGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACa                                                       >  1:49963/1‑140 (MQ=255)                                                                                     aCGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGATTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATCGTTGAGCCAt                                                      >  1:782916/1‑140 (MQ=255)                                                                                      cGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACAtt                                                     <  2:654360/140‑1 (MQ=255)                                                                                         tttAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACATTAAt                                                  >  1:452489/1‑140 (MQ=255)                                                                                          ttAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTGTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATCGGGTTAGAACGCTTATACTTGAGACATTCATg                                                 >  2:466848/1‑140 (MQ=255)                                                                                            aaCAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGCCATTCATCtt                                               >  1:755062/1‑140 (MQ=255)                                                                                            aaCAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACATTCATGtt                                               >  1:90625/1‑140 (MQ=255)                                                                                            aaCAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACATTCATGtt                                               >  1:481054/1‑140 (MQ=255)                                                                                                 tAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACATTCATGTTCAAcc                                          <  2:376727/140‑1 (MQ=255)                                                                                                  cgACTATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATACGGTTAGAACGTTTATAGTTGAGACATTCATGTTCAACCa                                         <  2:669519/139‑1 (MQ=255)                                                                                                      aaTTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAACTAGGGTTAGAAAGTTTATAGTTGAGACATTCATGTTCAACCAAAAt                                     >  1:171672/1‑140 (MQ=255)                                                                                                      aaTTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACATTCATGTTCAACCAAAAt                                     >  1:440077/1‑140 (MQ=255)                                                                                                        ttGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACATTCATGTTCAACCAAAAttt                                   >  2:582589/1‑140 (MQ=255)                                                                                                            ttttCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACACATGTATGTTTTTCACTCCTGCTTTAACATAGGGTTAGTAAGTGTATAGTTGAGACATTCATGTTCAACCAAACCTTTggt                               >  2:53621/1‑138 (MQ=255)                                                                                                             tttCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACATTCATGTTCAACCAAAATTTTGTTc                              <  1:282099/140‑1 (MQ=255)                                                                                                                 cGCTTGTTGTATTAAAAAAGGTTTTATAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACATTCATGTTCAACCAAAATTTTGTTCGAAt                          <  2:27918/140‑1 (MQ=255)                                                                                                                 cGCTTGTTGTATTAAAAAAGGTTTTATAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACATCCATGTTCAACCAAAATTTTGTTCGAAt                          <  2:622522/140‑1 (MQ=255)                                                                                                                    ttgttgGATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCCCTCCTTCTTTAAAATAGGGTTAGAAAGTTTTTAGTTGAGACATTCATGTTCAACCAAAATTTTGTTCGAATTCa                       <  2:199931/140‑1 (MQ=255)                                                                                                                          tATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACATTCATGTTCAACCAAAATTTTGTTCGAATTCAATAAAt                 <  2:201434/140‑1 (MQ=255)                                                                                                                           aTTAAAAAAGGTTTTGTAGATTGAGTATGAATATGCCATTCACTCAGGTATGTTTTTCACTCCTGCTTGAAAATAGGGTTAGAAAGTTTATAGTTGAGACATTCATGTTCAACCAAAATTTTGTTCGAATTCAATAAATg                <  1:320010/140‑1 (MQ=255)                                                                                                                            ttAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACATTCATGTTCAACCAAAATTTTGTTCGAATTCAATAAATGt               >  2:211753/1‑140 (MQ=255)                                                                                                                                aaaaGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACATTCATGTTCAACCAAAATTTTGTTCGAATTCAATAAATGTCTTg           >  2:417428/1‑140 (MQ=255)                                                                                                                                   aGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACATTCATGTTCAACCAAAATTTTGTTCGAATGCAATAAATGTCTTGttt        >  2:684843/1‑140 (MQ=255)                                                                                                                                         tGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACATTCATGTTCAACCAAAATTTTGTTCGAATTCAATAAATGTCTTGTTTAAAATa  <  1:91724/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                           AAACAAGTAAAGAATGATCGATGTTACTAATATTGCCATTCTCCAAATTTATTTCTGTGAGTATTTGGAAGCTACCATTAGGCAACGGTTTAACAATAGACAATTGCTTTTCCGCTTGTTGTATTAAAAAAGGTTTTGTAGATTGATTATTAATATGCCATTCACTCATGTATGTTTTTCACTCCTGCTTTAAAATAGGGTTAGAAAGTTTATAGTTGAGACATTCATGTTCAACCAAAATTTTGTTCGAATTCAATAAATGTCTTGTTTAAAATA  >  USA300TCH1516_ALE/881952‑882227

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