Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,062,480	T→C	9.0%	E40G (GAA→GGA)	USA300HOU_RS05245 ←	Acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,062,480	0	T	C	9.0%	163.5 / 2.0	61	E40G (GAA→GGA)	USA300HOU_RS05245	Acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (20/35);  new base C (3/3);  total (23/38)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATTCTACCCATTCCTTGTCCACGATGTGATTTCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTGATTGTTTACTTTTGAAAACATAAAATTCCTCCACATCCTAAC  >  USA300TCH1516_ALE/1062343‑1062619                                                                                                                                          |                                                                                                                                            gCATTCTACCCATTCCTTGTCCACGATGTGATTTCATCACAGCTACTCGTTATATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAATGGCAACTGGCTGTCCAGTATCATATGCAATGAGTTGAATagattaga                                                                                                                                          >  1:215275/1‑138 (MQ=255)    ttCTACCCATTCCTTGTCCACGATGTGATTTCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATagattcagatt                                                                                                                                       >  1:43411/1‑141 (MQ=255)          cgcttCCTTGTCATCGCTGTGATTTCATCACAGCTACTCGTCATATATTGACAGTTGTGGCATTAATAGGGCGTATTCGCGCAGTGGCAACTGGCTGTGCATGATCATCTCCAATGAGGTGAATAgattcggatccacctt                                                                                                                                 >  2:87953/4‑137 (MQ=255)          ccATTCCTTGTCCACGATGTGATTTCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGattcatatt                                                                                                                                 <  1:337710/141‑1 (MQ=255)          ccATTCCTTGTCCACGATGTGATTTCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGattcatatt                                                                                                                                 >  2:573948/1‑141 (MQ=255)             ttCCTTGTCCACTATGTGCTTTCATCACAGCTACTCGTTCTAGTCTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTTAATAGATTCAGattcatattcat                                                                                                                              <  1:202288/141‑1 (MQ=255)             ttCCTTGTCCACGATGTGATTTCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGattcatattcat                                                                                                                              <  1:202320/141‑1 (MQ=255)                  tGTCCACGATGTGATTTCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGCACCAGATTCATATTCATCAAtt                                                                                                                         >  1:78241/1‑141 (MQ=255)                    tCCACGATGTGATTTCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTc                                                                                                                       >  1:305202/1‑141 (MQ=255)                       aCGATGTGATTTCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACt                                                                                                                    <  1:3593/141‑1 (MQ=255)                       aCGATGTGAGTTCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACt                                                                                                                    <  1:736397/141‑1 (MQ=255)                        cGATGTGATTTCATCGAAGCTACCCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGAATTCGAGAAGTGGCAACTGGCTGTCCATTATCTTATCCAATGAGGTGAATAGATTCAGATTCATATTCACCAATTTCCCtt                                                                                                                   >  2:20240/1‑141 (MQ=255)                        cGATGTGATTTCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACtt                                                                                                                   >  1:456634/1‑141 (MQ=255)                        cGATGTGATTTCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACtt                                                                                                                   >  2:209082/1‑141 (MQ=255)                        cGATGTGATTTCATCACAGCTACTCGTTCTATTTTGACAGTTCTTTCTTTAATAGGGCGTATTCGGCCAGCGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACtt                                                                                                                   <  2:544785/141‑1 (MQ=255)                         gATGTGATTTCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACttt                                                                                                                  <  1:420214/141‑1 (MQ=255)                         gATGTCATTTCATCACAGCTACTCTTTCTATTTTGACAATTGTTTCAGTAATAGGCCGTATGCGAGCAGTGGCAACTGGCTGGCCATTATCATATCCAATGAGGTGAATAGATCCAGATTCATATTCATCAATTTCACttt                                                                                                                  <  2:925/141‑1 (MQ=255)                             tgATTTCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTc                                                                                                              <  1:308316/141‑1 (MQ=255)                               aTTTCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAg                                                                                                            <  1:204177/141‑1 (MQ=255)                                tttCATCACAGCTACTGGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAgg                                                                                                           >  1:459621/1‑141 (MQ=255)                                tttCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATCCAGATTCATATTCATCAATTTCACTTTCCTCAgg                                                                                                           >  1:332667/1‑141 (MQ=255)                                tttCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATCCAGATTCATATTCATCAATTTCACTTTCCTCAgg                                                                                                           >  2:421638/1‑141 (MQ=255)                                   aatcaCAGCTACTCGTTCTCACTTAAAGCTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGAc                                                                                                        >  2:165919/2‑141 (MQ=255)                                     tcaCAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGc                                                                                                      <  2:565537/141‑1 (MQ=255)                                      cacaGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGcc                                                                                                     <  1:587946/141‑1 (MQ=255)                                         aGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTg                                                                                                  <  1:3682/141‑1 (MQ=255)                                             aCTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGttct                                                                                              <  2:195548/141‑1 (MQ=255)                                                  ttcttctttGACAGTTGTTTCAGTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAAGTGAATAGATTCAGATCCATATTCATCAATTTAACTTTCCTCAGGGACGCCTTGTTCTTCTAc                                                                                         <  2:183972/135‑1 (MQ=255)                                                   tCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACa                                                                                        >  1:96959/1‑141 (MQ=255)                                                   tCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACa                                                                                        <  1:727944/141‑1 (MQ=255)                                                    cTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGGTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGCTTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTATaa                                                                                       <  2:568429/141‑1 (MQ=255)                                                        ttttACAGTTTTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACCGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAacac                                                                                   <  2:325614/141‑1 (MQ=255)                                                        tttGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTTCAAacac                                                                                   <  1:486968/141‑1 (MQ=255)                                                        tttGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAacac                                                                                   <  1:406019/141‑1 (MQ=255)                                                         ttGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACAct                                                                                  <  2:734004/141‑1 (MQ=255)                                                         ctggcAGTTGTTTCATTAATAGGGCGTATTCGAGCGGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACAct                                                                                  <  2:379307/137‑1 (MQ=255)                                                           taCAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCCTATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAAGTGCACTTTCCTCAGGGACGCCTTGTTCTTCGACAAACActct                                                                                <  2:358284/140‑3 (MQ=255)                                                           gACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACActtt                                                                                <  2:695402/141‑1 (MQ=255)                                                            aCAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACActttc                                                                               >  1:585592/1‑141 (MQ=255)                                                               gttgttTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCCGGGACGCCTTGTTGTTCTACAAACCctttcttt                                                                            >  1:493092/1‑141 (MQ=255)                                                                       aTTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTtatata                                                                    <  1:776947/141‑1 (MQ=255)                                                                        ttAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAg                                                                   <  1:7200/141‑1 (MQ=255)                                                                        ttAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAg                                                                   <  1:19383/141‑1 (MQ=255)                                                                        ttAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAg                                                                   >  1:199005/1‑141 (MQ=255)                                                                          aaTAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGaa                                                                 <  2:673125/141‑1 (MQ=255)                                                                          aaTAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGaa                                                                 <  1:397004/141‑1 (MQ=255)                                                                             aGGGCGTAGTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCa                                                              <  1:379636/141‑1 (MQ=255)                                                                                        gAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACa                                                   >  1:458367/1‑141 (MQ=255)                                                                                              tGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCtttt                                             >  1:358475/1‑141 (MQ=255)                                                                                               ggCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTg                                            <  1:573948/141‑1 (MQ=255)                                                                                                gCAACTGGCTGTCCATTATAATATCCAATGAGGTGAATAGATTCAGATTCATGTTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTTTATAGAAGCAACCTTCTAACACCttttct                                           >  2:534322/1‑139 (MQ=255)                                                                                                  aaCTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTGAtt                                         >  1:521983/1‑141 (MQ=255)                                                                                                   aCTGGCTGTCCATTATCATATCCAATGAGGTGCATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTCCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTGATTg                                        <  1:519255/141‑1 (MQ=255)                                                                                                                tATCATATCCAATGAGGTGAATAGATCCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTGATTGTTTACTTTTGaaa                           <  2:332667/141‑1 (MQ=255)                                                                                                                tATCATATCCAATGAGGTGAATAGATCCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTGATGGTTTACTTTTGaaa                           <  1:421638/141‑1 (MQ=255)                                                                                                                    atatCCAATGAGGTGAATAGATTCAGATTCGTATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTTTATAGTAGCAATCTTCTCACATCTTTTTATTGTTTATTTTTAaaaaaaa                       >  2:690109/1‑138 (MQ=255)                                                                                                                      atCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTGATTGTTTACTTTTGAAAACATaa                     <  2:78241/141‑1 (MQ=255)                                                                                                                      atCCAATCAGGTGAATAGATTCAGCTTCCTATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTGATTGTTTACTTTTGAAAACATaa                     >  1:492035/1‑141 (MQ=255)                                                                                                                       tCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTGATTGTTTACTTTTGAAAACATaaa                    >  2:470225/1‑141 (MQ=255)                                                                                                                            tGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTGATTGTTTACTTTTGAAAACATAAAATtcc               <  1:549839/141‑1 (MQ=255)                                                                                                                             gcggTGAATCGATTCAGACTCATATTCATCCATTTCACTTTCCGCACGGACGCCTTGTTCTTCTTCAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTGTTTGTGTACTCTAGAAAACATAAAATtcct              <  2:263615/139‑1 (MQ=255)                                                                                                                              aGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTGATTGTTTACTTTTGAAAACATAAAATtcctc             >  1:735479/1‑141 (MQ=255)                                                                                                                               ggTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTGATTGTTTACTTTTGAAAAAATAAAATtcctcc            >  1:114384/1‑141 (MQ=255)                                                                                                                                   aaTCGATTCAGATTCATCTTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATATAAGCAATCTTCTAACATCTTTTGATTGTTTACTTTTGAAAACATAAAATTCCTCCACAt        <  2:517609/141‑1 (MQ=255)                                                                                                                                      agattcagattcATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTGATTGTTTACTTTTGAAAACCTAAAATTCCTCCACATCCt     >  2:397759/1‑141 (MQ=255)                                                                                                                                       gatgcaggttcaTATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTGTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTGATTGTTTACTTTTGAAATCTTAAAATTCCTCCATATCCTa    >  2:99380/1‑141 (MQ=255)                                                                                                                                        attcagattcaTATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACAGCTTTTGACTGTTTACTTTTGAAAACATAAAATTCCTCCACATCCTaa   <  2:484775/141‑1 (MQ=255)                                                                                                                                         ttcagattcaTATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTGATTGTTTACTTTTGAAAACATAAAATTCCTCCACATCCTAAc  >  1:172168/1‑141 (MQ=255)                                                                                                                                          |                                                                                                                                            GCATTCTACCCATTCCTTGTCCACGATGTGATTTCATCACAGCTACTCGTTCTATTTTGACAGTTGTTTCATTAATAGGGCGTATTCGAGCAGTGGCAACTGGCTGTCCATTATCATATCCAATGAGGTGAATAGATTCAGATTCATATTCATCAATTTCACTTTCCTCAGGGACGCCTTGTTCTTCTACAAACACTTTCTTTCTTATATAGAAGCAATCTTCTAACATCTTTTGATTGTTTACTTTTGAAAACATAAAATTCCTCCACATCCTAAC  >  USA300TCH1516_ALE/1062343‑1062619

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