Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	110,660	T→A	21.2%	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	21.2%	138.3 / 30.8	68	I173I (ATT→ATA)	USA300HOU_RS00515	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (24/28);  new base A (6/8);  total (31/37)
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.99e-01

CGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGG  >  USA300TCH1516_ALE/110521‑110795                                                                                                                                            |                                                                                                                                        cGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATAtt                                                                                                                                         <  1:351749/140‑1 (MQ=255)              cGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                                                                            <  1:382931/140‑1 (MQ=255)              cGATAATAAAAACATACCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                                                                            <  1:560057/140‑1 (MQ=255)              cGAGAATAAAATCATACCAACAGAAGAAATGAGACATGAAAAAATACAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                                                                            <  2:200410/140‑1 (MQ=17)                  aataaAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGAAATTTCATATAATCCagagat                                                                                                                        >  2:557388/1‑138 (MQ=255)                    taaAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATGAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGACAAATTATCAAGACGGAGATATTTCATATAATCCAGAGGTGc                                                                                                                      >  1:18133/1‑140 (MQ=255)                    taaAATCATTCCAACAAAAGAAATTAAAGATGAAAAAAGAAAAAAAGAAATCGAAAACTTTAAGGTCTTGGGTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGc                                                                                                                      <  1:574892/140‑1 (MQ=255)                      aaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCg                                                                                                                    <  1:589426/140‑1 (MQ=255)                            ttCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGAAGGAGATATCTAATATAATCCAAAGGTACCGaataat                                                                                                              >  2:59002/1‑135 (MQ=17)                              ccAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTc                                                                                                            <  2:297916/140‑1 (MQ=255)                                       gAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCAACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGATGCCGAGTTATTCGGCTAAata                                                                                                   >  2:685917/1‑140 (MQ=255)                                          aTTAGAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCa                                                                                                <  2:20635/140‑1 (MQ=255)                                               aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATCGCGACAAACAGGATGTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                           <  1:697691/140‑1 (MQ=14)                                                    aaaaaaTAAAAAAAGAAATCGAAAACTTTGAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                                      <  2:450327/140‑1 (MQ=255)                                                      aaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                                    <  2:331310/140‑1 (MQ=255)                                                      aaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                                    <  2:633657/140‑1 (MQ=255)                                                      aaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACAGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                                    >  2:247526/1‑140 (MQ=255)                                                          aaaaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                                                                <  2:123092/139‑1 (MQ=255)                                                            aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGAttatta                                                                              >  1:480160/1‑140 (MQ=255)                                                            aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATta                                                                              >  2:166390/1‑140 (MQ=255)                                                            aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATta                                                                              >  1:603936/1‑140 (MQ=255)                                                             aaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCTGCTAAATATCAATTAACTAATGATGATtat                                                                             <  2:216822/140‑1 (MQ=255)                                                              aaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCATAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtata                                                                            >  2:512805/1‑140 (MQ=255)                                                                     aTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                                                                     >  2:494772/1‑140 (MQ=255)                                                                      tCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                                                                    >  1:472098/1‑140 (MQ=255)                                                                      tCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                                                                    >  2:551879/1‑140 (MQ=255)                                                                          aaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                <  1:112266/140‑1 (MQ=255)                                                                          aaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATACTTCATATAATCCAGAGGTGCCGAGCTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                <  2:351482/140‑1 (MQ=255)                                                                           acatttAAGTTATTTGTTCAATATGGGGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAGCAAtt                                                               <  2:229559/137‑1 (MQ=255)                                                                              tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGCTATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGACGATTATACTGTAAAACAATTACg                                                            >  2:272063/1‑140 (MQ=255)                                                                              tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAGAACAATTACg                                                            >  2:635719/1‑140 (MQ=255)                                                                               ttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                           <  2:367149/140‑1 (MQ=255)                                                                                 aaGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTACGTATCAATTAACTAATGATGATTATAATGTATCACAATTTCGCat                                                         >  1:360326/1‑139 (MQ=18)                                                                                          gTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatga                                                >  1:176486/1‑140 (MQ=255)                                                                                                      gACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAg                                    <  1:581843/140‑1 (MQ=255)                                                                                                         tttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                                 <  2:326765/140‑1 (MQ=11)                                                                                                         tttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTACATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                                 <  2:47169/140‑1 (MQ=255)                                                                                                         tttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAACATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                                 <  2:673157/140‑1 (MQ=255)                                                                                                            aaaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGc                              <  2:65547/140‑1 (MQ=255)                                                                                                             aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAGAGCt                             >  1:121044/1‑140 (MQ=255)                                                                                                              aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTc                            >  2:671701/1‑140 (MQ=255)                                                                                                               aaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCACAAGATATGATATAACGACGAGTAAAGcccc                           >  2:139458/1‑140 (MQ=255)                                                                                                                aTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCa                          <  1:618312/140‑1 (MQ=255)                                                                                                                  ttGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAAATATGATATACCGACGAGTAAAGCTCCaaa                        >  2:262967/1‑140 (MQ=255)                                                                                                                   tGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTACCGGTATCGCCGTCTTAATATAAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCCAAg                       >  1:194600/1‑140 (MQ=18)                                                                                                                    gatgaagtataAAGATGGAGATATATCATATAAGCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCCCCAAAGt                      <  2:305259/133‑1 (MQ=14)                                                                                                                    aatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGt                      <  2:556564/136‑1 (MQ=255)                                                                                                                    aatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGt                      <  2:106459/136‑1 (MQ=255)                                                                                                                     aaaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGtt                     <  1:531901/140‑1 (MQ=255)                                                                                                                      aaaaTTATAAAGACGGAGATATTTCATATAATCCAGAAGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATCTGATATACCGACGAGTAAAGCTCCACAGTTa                    >  2:702216/1‑140 (MQ=255)                                                                                                                      aaaaTTATAAAGACGGAGATATTTCATATAATACAGAGGTGCCAAGCTATTCGGCTGCATATCTATTAACTCATGATGATTATAATGTAAAACAATTACGCAACAGCTATGCTATAACGCCGATTACAGATCCAAAgctg                    >  2:318181/1‑137 (MQ=17)                                                                                                                       aaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAAGAAAGCTCCAAACTTAt                   >  2:144720/1‑140 (MQ=255)                                                                                                                         aTTATAGAGATGGAGATATATCATATAATCCAGAGGTGTCGAGCTAGTCGTCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTg                 >  1:102559/1‑140 (MQ=17)                                                                                                                         aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTg                 >  1:539115/1‑140 (MQ=255)                                                                                                                           tgtcaaGACGGAGATATTTCATAGAATCAAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGtt               <  2:409377/136‑1 (MQ=255)                                                                                                                             tagaGATGGAGATATGTCATATAATCAAGAGGTGCCGAGTGATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaa             <  2:480160/137‑1 (MQ=18)                                                                                                                             taAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaa             <  1:533963/140‑1 (MQ=255)                                                                                                                             taAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaa             <  1:143340/140‑1 (MQ=255)                                                                                                                             taAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaa             <  1:166390/140‑1 (MQ=255)                                                                                                                              aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaa            <  2:538461/140‑1 (MQ=255)                                                                                                                                agatggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAg          >  2:567316/5‑140 (MQ=255)                                                                                                                                  atggAGATATATCATATAATCCAGAGGTGCCGATTTATTCGGATAAATATCAATTAACTAATGATGATTATAATGTAACACAATTACGCAAAAGATATGATATACTGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGt        >  1:200030/3‑140 (MQ=18)                                                                                                                                  atggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGt        >  1:525948/3‑140 (MQ=255)                                                                                                                                  atggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGt        >  1:147143/3‑140 (MQ=255)                                                                                                                                   cGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAGGCTCCAAAGTTATTGGTAAAAGGtg       >  1:182051/1‑139 (MQ=255)                                                                                                                                    ggAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTc      <  2:190846/140‑1 (MQ=255)                                                                                                                                    ggAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTc      <  1:557388/140‑1 (MQ=255)                                                                                                                                       gaTATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAGTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAgg   <  1:166166/140‑1 (MQ=255)                                                                                                                                        agatTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAggg  <  2:601418/138‑1 (MQ=255)                                                                                                                                            |                                                                                                                                        CGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGG  >  USA300TCH1516_ALE/110521‑110795

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