Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	2,028,583	G→T	A21D (GCC→GAC)	USA300HOU_RS10230 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,028,583	0	G	T	98.7%	225.5 / NA	78	A21D (GCC→GAC)	USA300HOU_RS10230	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (33/43);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAATAGCGCGCCTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCATATAATTTATTAGATACGATGTTTTTTCTGTTGCTATTTTTTC  >  USA300TCH1516_ALE/2028446‑2028718                                                                                                                                          |                                                                                                                                        ttAATAGCGCGCCTATCTAACACCGCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCa                                                                                                                                       >  2:462211/1‑140 (MQ=255)  tAATAGCGCGAATATACAAAACGCCAAAAGAACACGAATTTTATCAAAATAAAGATATCCTACATAAACTAATACCCCACCAAATACTAGACCGACTAAAACGCCAATTCTGGCAAAAAATAGGAAGAAAAAATTGTCaa                                                                                                                                      <  2:586452/140‑1 (MQ=255)       gcgcgaCAATCAAAAACGCCAAAAGGAAACGCAGTTTAAGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAAAAAAAAAAAAAAATATGTAGTAAAAATTGTCAATAATc                                                                                                                                 <  2:800850/140‑1 (MQ=255)            ccTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATACCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTAAAAAAATATGTATTAAAAATTGTCAATAATCATTaa                                                                                                                            <  2:413946/140‑1 (MQ=255)             cTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCTATTTTTATGTAGTAAAAATTGGCAATCCTCATTAAt                                                                                                                           >  2:924681/1‑140 (MQ=255)             cTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAAtaattaat                                                                                                                           >  2:112341/1‑140 (MQ=255)               aTCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATgt                                                                                                                         >  1:597951/1‑139 (MQ=255)                tCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCa                                                                                                                        >  1:54090/1‑140 (MQ=255)                tCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCa                                                                                                                        >  2:30294/1‑140 (MQ=255)                   aaacCGCCAAGAGGCCACGCATGTTATGAAAATAAAGATATCCTACATAAAATAATACACATCCCAATAATAGAACGAGTAAACAGCCAATCTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTa                                                                                                                     <  1:290197/140‑1 (MQ=255)                   aaaaCGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTa                                                                                                                     <  1:42689/140‑1 (MQ=255)                         ccAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATg                                                                                                               >  1:695797/1‑140 (MQ=255)                            aaaGGCCATGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAAAAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATc                                                                                                            <  1:50172/140‑1 (MQ=255)                            aaaGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTCATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCCATAATCATTAATGCAGTAAAAATGATc                                                                                                            >  1:567834/1‑140 (MQ=255)                            aaaGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATTTGTAGTAAAAATTGTAAATAATCATTAATGAAGTAAAaataata                                                                                                            >  1:115222/1‑139 (MQ=255)                            aaaGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAGATTGTCAATAATCATTAATGCAGTAAAAATGAta                                                                                                            >  1:97371/1‑139 (MQ=255)                            aaaGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATc                                                                                                            <  1:878262/140‑1 (MQ=255)                            aaaGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATc                                                                                                            <  1:868651/140‑1 (MQ=255)                            aaaGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGACAATAATCATTAATGCAGTAAAAATGATc                                                                                                            >  2:539516/1‑140 (MQ=255)                            aaaGCCCACGTATTTTATGAAAAGAAAGATATCCATCATAAACTCAAACCCATCCAAATAATACACCAAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATc                                                                                                            <  2:859106/140‑1 (MQ=255)                                accACGTATTTTATAAAAATAAAGATATCCTACATAAACTAAGACAAATCCCAAAAATAGAACGAGTAAAAAGCAAATTTTTTCAAAAAATATTTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACa                                                                                                        <  2:450215/139‑1 (MQ=255)                                 ccACCTATTTGACGAAAATAAAGATATCCTACATAAACCAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACAt                                                                                                       <  2:75062/140‑1 (MQ=255)                                  cACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACAtt                                                                                                      <  1:1002151/140‑1 (MQ=255)                                  cACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACAtt                                                                                                      <  1:253125/140‑1 (MQ=255)                                     gTATTTTATGAAAATACAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCATAAAATAGGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAATTAATCAACATTTGc                                                                                                   >  1:558146/1‑140 (MQ=255)                                     gTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGc                                                                                                   >  1:309156/1‑140 (MQ=255)                                     gTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGc                                                                                                   >  1:593472/1‑140 (MQ=255)                                      tATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCg                                                                                                  <  1:306225/140‑1 (MQ=255)                                       aTTTTAGGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGCACGAGTAAAAAGCAAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGt                                                                                                 <  2:286859/140‑1 (MQ=255)                                           tATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGat                                                                                             >  1:977230/1‑140 (MQ=255)                                              gAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGatata                                                                                          <  2:878482/140‑1 (MQ=255)                                                aaaTAAAGATATCCTACATAAACCAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCCATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGatatata                                                                                        <  2:5729/140‑1 (MQ=255)                                                     aaGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTtgt                                                                                   <  2:696179/140‑1 (MQ=255)                                                     aaGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCAAAATATTAAAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTtgt                                                                                   <  2:723884/140‑1 (MQ=255)                                                        atatCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTtgtgtg                                                                                >  2:298208/1‑140 (MQ=255)                                                        atatCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTtgtgtg                                                                                >  1:780668/1‑140 (MQ=255)                                                          atCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTc                                                                              <  1:575376/140‑1 (MQ=255)                                                           accTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCa                                                                             <  1:875768/139‑1 (MQ=255)                                                               aCATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGtt                                                                         >  2:474144/1‑140 (MQ=255)                                                               aCATAAACTAATACACATCACAATAATAGAACGAATAATAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGtt                                                                         <  2:855101/140‑1 (MQ=255)                                                               aCATAAACTAATACACATCACAAAAATACAACGAGTATAAAGAAAATTAATTAAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGtt                                                                         <  2:855203/140‑1 (MQ=255)                                                                cATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGttt                                                                        >  2:954339/1‑140 (MQ=255)                                                                  tAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTcg                                                                      <  1:228307/140‑1 (MQ=255)                                                                  aaaaGTAATACACATCCCAATTATAGAACGAGTAAAAAGCCCATTTTTTGAAACAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTcg                                                                      <  1:77182/139‑1 (MQ=255)                                                                     tcTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTcaccg                                                                   <  2:886739/139‑1 (MQ=255)                                                                     actatgacacAGCCCAATAATAGAACGAGTAAAAAGCCAATTTGTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTcgccg                                                                   <  2:116770/134‑1 (MQ=255)                                                                        aaTACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGAAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCt                                                                >  2:683417/1‑140 (MQ=255)                                                                          tACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCttt                                                              >  1:215213/1‑140 (MQ=255)                                                                          tACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCACCttt                                                              >  1:216916/1‑140 (MQ=255)                                                                            cacaTCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCt                                                            >  1:785528/1‑140 (MQ=255)                                                                            cacaTCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCt                                                            >  1:655189/1‑140 (MQ=255)                                                                             acaTCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTTTTTTTTTTTGTAGAAAAAATTGTCAACAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATATTTCCCCGCCTTTCTa                                                           >  2:868856/1‑140 (MQ=255)                                                                               ataccAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGAGTTGATATATATTTGTGTGTCATGTTTCGCTGCCTTTCTAtt                                                         >  2:274891/4‑139 (MQ=255)                                                                                tCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGc                                                        <  1:148034/140‑1 (MQ=255)                                                                                 cctaataatAGAAAGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGcc                                                       <  1:224962/137‑1 (MQ=255)                                                                                 cccAATAATAGAACGAGTAAAAAGCGAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATCTTTCGCCGCCTTTCTATGcc                                                       <  1:610323/140‑1 (MQ=255)                                                                                  ccAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGGCGCCTTTCTATGCca                                                      >  1:55003/1‑140 (MQ=255)                                                                                  ccAATAAAAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGGAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCca                                                      <  1:360456/140‑1 (MQ=255)                                                                                    aataatAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCcaca                                                    >  1:25420/1‑140 (MQ=255)                                                                                       aataGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCg                                                 <  1:210388/140‑1 (MQ=255)                                                                                       aataGAACCAGTAAAAAGCAAATGTTTTCAAAAAACATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTACGTGGATATATATTTGTGTTTCATGTTTCGCCGCCTTTCTATGCCACAGCg                                                 <  1:363301/140‑1 (MQ=255)                                                                                         taGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTCTCTATGCCACAGCGtt                                               <  1:598975/140‑1 (MQ=255)                                                                                            aatgaatAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCat                                            <  2:22554/134‑1 (MQ=255)                                                                                            aaCGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCTTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCat                                            <  2:695797/140‑1 (MQ=255)                                                                                              cgtctacaaaaCCAATGTTTTCATATAATTTGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACTTTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCatat                                          <  2:529093/133‑1 (MQ=255)                                                                                              cGAGTAAAAAGCCAATTATTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCatat                                          <  1:46214/140‑1 (MQ=255)                                                                                              aaaaaaaaaaGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCatat                                          <  2:223172/135‑1 (MQ=255)                                                                                                  tAAAAAGCCAATTTTTTCAAATAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCATATAAtt                                      >  1:522976/1‑140 (MQ=255)                                                                                                     aaaGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCATATAATttat                                   >  2:889665/1‑140 (MQ=255)                                                                                                     aaaGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCATATAATttat                                   >  2:889624/1‑140 (MQ=255)                                                                                                       aGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCATATAATttatta                                 <  1:298208/140‑1 (MQ=255)                                                                                                              tttttCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTTTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCATATAATTTATTAGATACGa                          >  1:699986/1‑140 (MQ=255)                                                                                                              tttttCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCATATAATTTATTAGATACGa                          <  1:75924/140‑1 (MQ=255)                                                                                                                 ttCAAAAAATATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCATTTCTATGCCACAGCGTTCATATAATTTATTAGATACGATGt                       >  2:862242/1‑140 (MQ=255)                                                                                                                   cAAAAAATATGTAGTAAATATTGTGAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCATATAATTTATTAGATACGATGttt                     <  1:19290/140‑1 (MQ=255)                                                                                                                     aaaaaTATGTAGTAAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCATATAATTTATTAGATACGATGttttt                   >  1:394830/1‑140 (MQ=255)                                                                                                                             gtagtaAAAATTGTCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCATATAATTTATTAGATACGATGTTTTTTCTGTTGc           >  1:824916/1‑140 (MQ=255)                                                                                                                                      aTTGTCAATAATCATTAATGCAGTAAATATGATCAACAGTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCATATAATTTATTAGATACGATGTTTTTTCTGTTGCTATTTTTTc  <  1:620794/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                        TTAATAGCGCGCCTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCATATAATTTATTAGATACGATGTTTTTTCTGTTGCTATTTTTTC  >  USA300TCH1516_ALE/2028446‑2028718

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