Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	1,872,184	G→C	Q381E (CAA→GAA)	USA300HOU_RS09300 ←	Putative dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,872,184	0	G	C	98.9%	308.6 / NA	91	Q381E (CAA→GAA)	USA300HOU_RS09300	Putative dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (55/34);  minor base T (1/0);  total (57/34)
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

CAAATGCTACACCCTTGTCTAAGTTTCTCGCATAAGTACCGCCACCTATAGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTGGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACA  >  USA300TCH1516_ALE/1872040‑1872326                                                                                                                                                 |                                                                                                                                               cAAATGCTACACCCTTGTCTAAGTTTCTCGCATAAGTACCGCCACCTATAGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGAc                                                                                                                                             >  2:817425/1‑148 (MQ=255)  acaTGCTACACCCTTGTCTAAGTTTCTCGCATAAGTACCGCCACCTATAGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAGCATAATGTGGTGGTTCGACt                                                                                                                                            <  1:223745/146‑1 (MQ=255)   aaTGCCACCCCCTTTTCTAAGTTTCTCGCATAAGTACCGCCACCTATAGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACtt                                                                                                                                           <  2:901880/148‑1 (MQ=255)               ttGTCTAAGTTTCTCGCATAAGTACCGCCACCTATAGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCa                                                                                                                               >  2:423932/1‑148 (MQ=255)                tGTCTAAGTTTCTCGCATAAGTACCGCCACCTATAGTATAAGGTTCAGTCATCTCAAGTGTTGGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCAc                                                                                                                              <  1:472956/148‑1 (MQ=255)                      agattCACGCATAAGTACCGCCACCTATAGTAGAAAGTCCAGTCATATCATTTGTTTGAGTTCTATATGCAGTAACTAACTTTTGTACTAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCaaa                                                                                                                        <  1:897110/145‑1 (MQ=255)                      aGTTTCTCGCATAAGTACCGCCACCTATCGTATTAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCaaa                                                                                                                        <  1:922660/148‑1 (MQ=255)                      aGTTTCTCGCATAAGTACCGCCACCTATAGTATAAGGTTCAGTGCTATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATGTTTCTCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCaaa                                                                                                                        <  2:343166/148‑1 (MQ=255)                      aGTTTCTCGCATAAGTACCGCCACCTATAGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCaaa                                                                                                                        <  2:940017/148‑1 (MQ=255)                      aGTTTCTCGCATAAGTACCGCCACCTATAGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCaaa                                                                                                                        <  1:652657/148‑1 (MQ=255)                      aGTTTCTCGCATAAGTACCGCCACCTATAGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCaaa                                                                                                                        <  1:710235/148‑1 (MQ=255)                      aGTTTCTCGCATAAGTACCGCCACCTATAGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCaaa                                                                                                                        <  1:881685/148‑1 (MQ=255)                       gTTTCTCGCATAAGTACCGCCACCTATAGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCaacc                                                                                                                       >  2:225329/1‑146 (MQ=255)                       gTTTCTCGCATAAGTACCGCCACCTATAGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTAGACTTTAACTAATTTCACTTCaaac                                                                                                                       >  2:154925/1‑147 (MQ=255)                        tttCTCGCATAAGTACCGCCACCTATAGTATAAGGGTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTACTTTCACTTCAAAGc                                                                                                                      >  2:308729/1‑148 (MQ=255)                         ttCTCGCATAAGTACCGCCACCTATAGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGcc                                                                                                                     >  1:855641/1‑148 (MQ=255)                          tctgGCCTAAGTACCGCCACCTATACTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTCATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACGTCAAAGCCa                                                                                                                    <  2:856407/148‑1 (MQ=255)                          tctcGCATAAGTACCGCCACCTATAGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACGAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTTGACTTTACCTAATTTCACTTCAAAGCCa                                                                                                                    >  1:689466/1‑148 (MQ=255)                               cATAAGTACCGCCACCTATAGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATAttg                                                                                                               >  2:400445/1‑148 (MQ=255)                                                 taGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCaaa                                                                                             >  1:411896/1‑148 (MQ=255)                                                  aGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAGACGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTTCCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCaaaa                                                                                            >  1:433829/1‑148 (MQ=255)                                                  aGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTGTTGTACAAAAGGATCATTTTTATCAACATAATTTGTTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTAAATCTCCTTTCCaaaa                                                                                            >  2:482213/1‑148 (MQ=255)                                                    tataAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACg                                                                                          >  1:72764/1‑148 (MQ=255)                                                    tataAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAATAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTACTTCCACTTCAAAGCCATATTGTTGAACACAATTTGCAAAACg                                                                                          >  2:762908/1‑148 (MQ=255)                                                     ataAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCGAAGCCATATTGTTGAATCTCATTTGCAAAACGa                                                                                         >  1:673544/1‑148 (MQ=255)                                                     ataAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAACACGa                                                                                         >  2:654954/1‑148 (MQ=255)                                                      taAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGaa                                                                                        >  2:109530/1‑147 (MQ=255)                                                      taAGGTTCAGTCATATCATTTGTTTGAGTTCTATATGGAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACTAt                                                                                        >  1:37957/1‑148 (MQ=255)                                                       aaGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATc                                                                                       <  2:573731/148‑1 (MQ=255)                                                       aaGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATc                                                                                       >  1:915629/1‑148 (MQ=255)                                                       aaGCTTCACTCATTACCTTTGGTTGATTTCTATATGCAGGAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATc                                                                                       <  2:932566/148‑1 (MQ=255)                                                         ggTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCa                                                                                     >  2:74249/1‑148 (MQ=255)                                                         ggTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAAAGATCCa                                                                                     >  2:832708/1‑148 (MQ=255)                                                          gTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCAt                                                                                    >  2:755765/1‑148 (MQ=255)                                                          gTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACGTTTGTACAAAAGGGTCATTTGTATCAACATAATGTGGTGGTTCGAGTTTACCTAATTTCACTGCAAAGCTATATTCTTGGATCTGATTTGCAAGCCGATCCaa                                                                                    >  2:25362/1‑147 (MQ=255)                                                            tcggTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCCTTTGCAAAACGATCCATAg                                                                                  >  1:890299/4‑148 (MQ=255)                                                             cAGCCATATCATTTGTTTGATTTCTATATGCAGGAACTAACTTTTGTACAAAAGGATCATTTTTACCAACATAATGTGGGGGCTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGc                                                                                 <  1:146951/148‑1 (MQ=255)                                                              aGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCt                                                                                >  2:563767/1‑148 (MQ=255)                                                                  atatCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCttttt                                                                            >  1:757483/1‑148 (MQ=255)                                                                  atatCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTGATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCtgttt                                                                            >  1:438515/1‑148 (MQ=255)                                                                    atCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTttca                                                                          >  1:614837/1‑148 (MQ=255)                                                                    atCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTttca                                                                          >  1:148036/1‑148 (MQ=255)                                                                     tCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTttcaa                                                                         >  2:835784/1‑148 (MQ=255)                                                                      cATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTttcaaa                                                                        >  2:714902/1‑148 (MQ=255)                                                                         ttgtttgATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTtttatattt                                                                     >  2:832197/1‑147 (MQ=255)                                                                         ttgtttgATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTttcaaattc                                                                     >  1:909791/1‑148 (MQ=255)                                                                         ttgtttgATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTttcaaattc                                                                     >  2:859007/1‑148 (MQ=255)                                                                            tttgATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTttcaaattcaaa                                                                  <  2:767774/148‑1 (MQ=255)                                                                            tttgATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTttcaaattaaaa                                                                  >  2:188634/1‑148 (MQ=255)                                                                             ttgATTTCTATATGCAGTAACTAACTTTTGTACACAAGTATCATTTTCATCAGCATAATGGGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTtcaaattcaaat                                                                 <  2:349928/148‑1 (MQ=255)                                                                              tgATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCCTATTGTTGAATCTCATTTGCAAAACGTTCCATATCTTTTTCAAATTCAAATc                                                                >  2:783207/1‑148 (MQ=255)                                                                                aTTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCt                                                              >  1:920329/1‑148 (MQ=255)                                                                                aTTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCt                                                              >  1:435816/1‑148 (MQ=255)                                                                                 tttCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCtt                                                             <  1:155569/148‑1 (MQ=255)                                                                                  ttCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTc                                                            >  2:907704/1‑148 (MQ=255)                                                                                   tCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCt                                                           >  2:161428/1‑148 (MQ=255)                                                                                       tatGCAGTAACTAACCTTTGTACAAAAGGATCATTTTTATCAACCTAATTTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGt                                                       <  1:633554/148‑1 (MQ=255)                                                                                         tGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTGAAATCCTTCTGGGTAg                                                     >  1:128032/1‑148 (MQ=255)                                                                                           cAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGCATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCg                                                   >  2:792994/1‑148 (MQ=255)                                                                                            aGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGt                                                  >  2:313630/1‑148 (MQ=255)                                                                                             gTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTTCCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAACCGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGATAGCGTa                                                 >  2:701099/1‑148 (MQ=255)                                                                                             gTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGCTTCAGGCCGTTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTa                                                 >  1:721362/1‑148 (MQ=255)                                                                                                 ctaactTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGtt                                             >  2:13284/1‑148 (MQ=255)                                                                                                 ctaactTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGtt                                             >  1:11723/1‑148 (MQ=255)                                                                                                 ctaactTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTCCTGGGAAGCGTAAGtt                                             >  1:419020/1‑148 (MQ=255)                                                                                                    actTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGAt                                          >  1:652416/1‑148 (MQ=255)                                                                                                            aCAAAAGGATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGaaaa                                  <  1:381392/148‑1 (MQ=255)                                                                                                                  ggATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTg                            >  2:829572/1‑148 (MQ=255)                                                                                                                   gATCATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGc                           >  2:116035/1‑148 (MQ=255)                                                                                                                      cATTTTTATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGtt                        <  2:855641/148‑1 (MQ=255)                                                                                                                         ttttATCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTc                     <  2:689466/148‑1 (MQ=255)                                                                                                                             aTCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTa                 <  1:400445/148‑1 (MQ=255)                                                                                                                              tCAACATAATGTGGTGGTTCGTCTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTAt                >  1:548480/1‑148 (MQ=255)                                                                                                                              tCAACATAATGTGGTGGTTCGACTTTACCTAGTTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAATAGACCTGCGTTTTCATTAt                >  1:634123/1‑148 (MQ=255)                                                                                                                              tCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTGAGTTGATACCGAAAAGACCTGCGTTTTCATTAt                <  2:614837/148‑1 (MQ=255)                                                                                                                              tCAACATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTAt                <  2:148036/148‑1 (MQ=255)                                                                                                                                 aCATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTTATTTGCAAAACGATCCATAGCTTTTTCAAATTCAGATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCat             >  1:986211/1‑148 (MQ=255)                                                                                                                                 aCATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTTATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTTTGGGTAGCGTAAGTTGTTACCGAAAAGACCTGCGTTTTCATTATCat             >  1:302793/1‑148 (MQ=255)                                                                                                                                  cATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAACCGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCata            <  1:330099/148‑1 (MQ=255)                                                                                                                                  cATAATGTGGTGGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAACACCTGCGTTTTCATTCTCata            <  2:673544/148‑1 (MQ=255)                                                                                                                                        gtggtgGTTCGACTTTACCTTATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGtaat      <  2:107873/148‑1 (MQ=255)                                                                                                                                        gtggtgGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGtaat      <  2:172276/148‑1 (MQ=255)                                                                                                                                        gtggtgGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGtaat      >  1:833837/1‑148 (MQ=255)                                                                                                                                        gtggtgGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGATTTCATTATCATATGtaat      >  1:366230/1‑148 (MQ=255)                                                                                                                                        gtggtgGTTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGGAGCGTAAGTTGATACCGCAAAGACCTGAGTTTTCATTATCATATGtaat      >  1:621606/1‑148 (MQ=255)                                                                                                                                        ggggtggtTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGACTCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGtaat      <  2:201234/146‑1 (MQ=255)                                                                                                                                           gtggtTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAac   <  1:829572/148‑1 (MQ=255)                                                                                                                                            tggtTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCTTTTGCAAAACGATCCATAGACTTGTCAAATTCACATCCTTCTTTGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAaca  <  1:975891/148‑1 (MQ=255)                                                                                                                                            tggtTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAAACCTGCGTTTTCATTATCATATGTAATAaca  <  1:631998/148‑1 (MQ=255)                                                                                                                                            gggtTCGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAACTCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAaca  <  2:949351/147‑1 (MQ=255)                                                                                                                                            gggtTCGACTTTACCGAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAaca  <  1:936404/147‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                               CAAATGCTACACCCTTGTCTAAGTTTCTCGCATAAGTACCGCCACCTATAGTATAAGGTTCAGTCATATCATTTGTTTGATTTCTATATGCAGTAACTAACTTTTGTACAAAAGGATCATTTTTATCAACATAATGTGGTGGTTGGACTTTACCTAATTTCACTTCAAAGCCATATTGTTGAATCTCATTTGCAAAACGATCCATAGCTTTTTCAAATTCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACA  >  USA300TCH1516_ALE/1872040‑1872326

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