Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	110,660	T→A	10.4%	I173I (ATT→ATA)	USA300HOU_RS00515 →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	110,660	0	T	A	10.4%	179.3 / 10.1	66	I173I (ATT→ATA)	USA300HOU_RS00515	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (24/35);  new base A (3/4);  total (27/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAAT  >  USA300TCH1516_ALE/110522‑110798                                                                                                                                           |                                                                                                                                           ggTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAAAAATGGCTATTATAAAGACGGAGATAttt                                                                                                                                           <  1:660476/140‑1 (MQ=14) ggTTAAAATAGTCTATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATAttt                                                                                                                                           <  2:633028/140‑1 (MQ=255)  gTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTc                                                                                                                                          <  1:627446/140‑1 (MQ=255)    tAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCat                                                                                                                                        <  2:607879/140‑1 (MQ=255)      aaaaGGTCAATAATAAAATCATTCCAACAGAATAAATTACAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCatat                                                                                                                                      <  2:371993/140‑1 (MQ=17)         tGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTAAAAAATTATAAAGACAGAGATATTTAATATAAt                                                                                                                                   >  1:336940/1‑140 (MQ=255)              gATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATTATTCATTTAATCCaat                                                                                                                              >  2:359454/1‑138 (MQ=255)                  acaaaaGCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTg                                                                                                                          <  2:256923/138‑1 (MQ=255)                    aaaaataTTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                                                                        <  2:274670/134‑1 (MQ=255)                    aaaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                                                                        <  1:267580/140‑1 (MQ=255)                    aaaaTCACTCCAACAAAAGTAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                                                                        <  1:278493/140‑1 (MQ=255)                     aaaTTATTCTATCAAAAGAAATCAACGATCAAAAAATCAAAAAAGGATTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGATTTTATAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCg                                                                                                                       <  1:15838/140‑1 (MQ=17)                     aaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCg                                                                                                                       >  1:124396/1‑140 (MQ=255)                        tCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGt                                                                                                                    <  1:306927/140‑1 (MQ=255)                                   aaaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATATGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaa                                                                                                         >  1:485955/1‑140 (MQ=255)                                    aaGAAATTAAAGATGAAGAAATAAAAAAAGAAATCGAAAACCTTAATTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaaa                                                                                                        <  2:328966/140‑1 (MQ=255)                                    aaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaaa                                                                                                        <  2:601839/140‑1 (MQ=255)                                    aaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCAGAGTTATTCGGCTaaa                                                                                                        >  2:182931/1‑140 (MQ=255)                                      gAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATGATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                                                                      <  1:1754/140‑1 (MQ=255)                                            aaaGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                                >  2:97442/1‑140 (MQ=255)                                            aaaGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                                >  2:460542/1‑140 (MQ=255)                                            aaaGATAAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                                <  2:180169/140‑1 (MQ=255)                                             aagaataaaaaaTAAAAAAAGAAATCTAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                                                                               <  2:321233/134‑1 (MQ=255)                                              ataTGAAAAAATAAAAAAAGAAATCGATAACTTTAAGTTCTTTGTTCAATATTGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTGATTTTGCTAAATATCAATTaa                                                                                              <  2:144190/138‑1 (MQ=17)                                              aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTCGAAAAATCATCCCGACCGCGCTATTCCATATCCTCCAGCGCTCCCCAGTCACTCGGCTAACTCTCaaccaa                                                                                              >  2:139520/1‑136 (MQ=12)                                                aTGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                            <  1:106827/140‑1 (MQ=255)                                                 agAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTGTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTCCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                           <  2:168978/139‑1 (MQ=255)                                                   aaaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGAGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                                         <  1:550340/140‑1 (MQ=255)                                                      ataTAAAAAAAGAAATCGAAAACTTTAAGTTCTGTGTTCAATATGGCGTCTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                                      <  1:317793/138‑1 (MQ=255)                                                       aaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgatg                                                                                     <  2:234261/140‑1 (MQ=255)                                                        aTAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATTTTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAtgatga                                                                                    <  1:200981/140‑1 (MQ=255)                                                           aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATta                                                                                 >  2:580704/1‑140 (MQ=255)                                                              aaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGCAAAATTTCATATAATACAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTaaaa                                                                              >  1:155071/1‑140 (MQ=17)                                                              aaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATaa                                                                              >  1:542904/1‑140 (MQ=255)                                                               aaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAAt                                                                             >  1:641631/1‑140 (MQ=255)                                                                aGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATg                                                                            <  1:643105/140‑1 (MQ=255)                                                                  aaaTCGAAAACTTTCAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTGTAATGTa                                                                          <  1:416035/140‑1 (MQ=255)                                                                  aaaTCGAAAACTTTAAGTTCTTTGTTCAATACGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                                                                          <  2:16744/140‑1 (MQ=255)                                                                             tttAAGTTCTTTGTTCAATATGGTGACTTTAAAAATTTTAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGTTGCCGAGTTATTCGGCTAAATATCAATTAACTAAGGATGATTATAATGTAAAACAATTACg                                                               <  1:205390/140‑1 (MQ=255)                                                                                 atttCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGATGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa                                                           <  1:78548/138‑1 (MQ=255)                                                                                         gTTCAATATGGCGGCTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatga                                                   >  1:662986/1‑140 (MQ=255)                                                                                                  ggCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgac                                          >  2:591227/1‑140 (MQ=255)                                                                                                   gCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacg                                         <  2:382504/140‑1 (MQ=255)                                                                                                        tttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACATCTATCCAAAAGATATGATATACCGACGAGTaa                                    <  1:159177/140‑1 (MQ=255)                                                                                                          tAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAGcgacgaacaaca                                  >  2:327056/1‑134 (MQ=255)                                                                                                          tAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAg                                  >  2:228524/1‑140 (MQ=255)                                                                                                           aaaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGc                                 >  2:644096/1‑140 (MQ=255)                                                                                                           aaaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGc                                 <  1:255915/140‑1 (MQ=255)                                                                                                             aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTc                               >  2:82103/1‑140 (MQ=255)                                                                                                              aaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGATTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                              <  1:244433/140‑1 (MQ=255)                                                                                                              aaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATAAACCGACGAGTAAAGCTcc                              >  1:285008/1‑140 (MQ=255)                                                                                                               aTTTGAAAAATCATAAAGACGTGTGTGCTTTATATAATCCAGAAGCGCCTAATTATTCGGCTAAACATCAATTAACTAATGATGATTACAATGTAAAACAATTACGCAAAAGATATGAGATACCGACGAGTAAAGCTCCa                             <  2:513341/140‑1 (MQ=17)                                                                                                                tttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaa                            >  2:509017/8‑140 (MQ=14)                                                                                                                tttGAAAATTTATTAAGACGGAGATATTTCATCTAATCCAGTGATGCCGATTTATTCGTCTAAATTTCAATTACCCTATTACGATTATGATGTAAAACAACTCCGCATTATGTATGATATAACGATGAGTAGAGCTCCat                            <  2:207072/140‑2 (MQ=12)                                                                                                                     tgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGGCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa                       <  1:655658/138‑1 (MQ=255)                                                                                                                      gaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt                      <  1:60097/139‑1 (MQ=255)                                                                                                                      aaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATAGCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt                      >  2:335105/1‑140 (MQ=255)                                                                                                                         ttATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGt                   <  1:527809/140‑1 (MQ=255)                                                                                                                         ttATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGt                   >  1:548564/1‑140 (MQ=255)                                                                                                                              aaGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCTAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaaa              <  1:534995/140‑1 (MQ=255)                                                                                                                               aGACGGAGATATTTCATATAATCCAGATGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAg             <  1:510600/140‑1 (MQ=255)                                                                                                                               aGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGGGTAAAGCTCCAAAGTTATTGTTAAAAg             >  2:445946/1‑140 (MQ=255)                                                                                                                                  cGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGtt          <  2:124396/140‑1 (MQ=255)                                                                                                                                   ggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTc         <  1:92868/140‑1 (MQ=255)                                                                                                                                    gagaTATATCATATAATCCAGAGGTGCCGAGTTATTTGTCTAAGTCTAAAATGACTAGATATTGGTATGATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTAATTGTTAAAAGGTTCa        >  2:303962/1‑140 (MQ=17)                                                                                                                                      gaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAgg      >  1:459526/1‑140 (MQ=255)                                                                                                                                       atatTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAggg     >  2:330405/1‑140 (MQ=255)                                                                                                                                        tatTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGa    >  2:334682/1‑140 (MQ=255)                                                                                                                                          tttCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATTTACCGACGAGTAAAGCTCCAACGTTATTGTTAATAGGTTGGGGGaag  >  1:594006/1‑139 (MQ=18)                                                                                                                                           |                                                                                                                                           GGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAAT  >  USA300TCH1516_ALE/110522‑110798

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