Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	1,659,541	A→T	Y42N (TAT→AAT)	glcK ←	Glucokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,659,541	0	A	T	98.7%	257.0 / ‑3.4	75	Y42N (TAT→AAT)	glcK	Glucokinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (38/36);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

GCTCCATTTACTGTACCTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATATCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAGCTGCTAAAATAATTTTTGTCATTTCTGTT  >  USA300TCH1516_ALE/1659405‑1659671                                                                                                                                         |                                                                                                                                   gCTCCATTTACTGTACATTTTGCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAATTTCCTTTCAAAAGTGTATTTCCt                                                                                                                                >  1:916234/1‑141 (MQ=255) gCTCCATTTACTGTAACTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCt                                                                                                                                >  1:916237/1‑141 (MQ=255)    ccATTTACTGTACCATTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTa                                                                                                                             >  1:492625/1‑141 (MQ=255)    ccATTCACTGTACCTTGTTAAGAGTCAACAGGACCTGGTACACCTATACCTACGCCAAGTACATTTGAAACATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTa                                                                                                                             <  2:444269/141‑1 (MQ=255)      aTTTACTGTACCTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACt                                                                                                                           <  2:876863/141‑1 (MQ=255)      aTTTACTGTACCTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACt                                                                                                                           <  2:821761/141‑1 (MQ=255)       tttACTGTACCTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTa                                                                                                                          >  1:371569/1‑141 (MQ=255)        ttACTGTACCTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTAt                                                                                                                         <  1:162792/141‑1 (MQ=255)        ttACTGAACCTTTTTCAAAGTCAATATGACCTGGTACACCAATAACTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTAt                                                                                                                         <  2:575037/141‑1 (MQ=255)          aCTGTACCTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCa                                                                                                                       <  1:214012/141‑1 (MQ=255)           cTGTACCTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCACCAAAGGAAACATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAg                                                                                                                      >  1:860296/1‑141 (MQ=255)                 cTTTTTCAAAGTCAACAGGAGCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCTTGTACTATCAGATGTAt                                                                                                                >  1:647789/1‑141 (MQ=255)                 cTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTAt                                                                                                                <  1:644580/141‑1 (MQ=255)                     tttaaaTTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAgt                                                                                                            <  1:164191/138‑1 (MQ=255)                     ttCAAAGTCAACATGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTACCAgt                                                                                                            <  2:701281/141‑1 (MQ=255)                     ttCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTGCTATCAGATGTATCCgt                                                                                                            >  1:470185/1‑141 (MQ=255)                       cAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAAATGTAACAgatt                                                                                                          >  2:45777/1‑138 (MQ=255)                         aaGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGa                                                                                                        <  2:587413/141‑1 (MQ=255)                         aaGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGa                                                                                                        <  2:918167/141‑1 (MQ=255)                         aaGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGa                                                                                                        >  1:576390/1‑141 (MQ=255)                         aaGTCAACAGGACCTGGTACACCAATACCTACGCCAAATACATTTGAAAAATTATAATTTTTTTCATTTACTTTTTCAACAAACGAATCATAACTTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGa                                                                                                        <  1:54630/141‑1 (MQ=255)                          aGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTATTATTAGATGTATCAGTTTGat                                                                                                       >  2:314755/1‑140 (MQ=255)                          aGTCAACAGGACCCGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGaa                                                                                                       <  1:110306/141‑1 (MQ=255)                          aGTCAACAGGAACTGATACACCCATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAAGAAACGAATCATAAATTTCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGtgtgca                                                                                                       <  2:239916/141‑3 (MQ=255)                           gTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCCGTGTGAAt                                                                                                      >  1:473012/1‑141 (MQ=255)                           gTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTACAAGTGTGAAt                                                                                                      >  1:167155/1‑141 (MQ=255)                               aCAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTTTACTATCAGATGTATCAGTGTGAATAGaa                                                                                                  >  2:532470/1‑140 (MQ=255)                               aCAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTCTCAGTGTGAATAGAc                                                                                                  >  2:119885/1‑141 (MQ=255)                               aCAGGACCTGGTACACCAATACCTACGCCAAGTACATGTGAAAAATTATAATTATTTTCATTTACTTATTCAACAAACGAATGATAAATTCCTTTCAAAAGTATAATTCCTGTACTATCACATGTATCAGTTTGACTAGAc                                                                                                  >  2:532397/1‑141 (MQ=255)                                cAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCTTCGTTTGCTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGAcc                                                                                                 <  1:437800/141‑1 (MQ=255)                                cAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGAcc                                                                                                 <  1:141296/141‑1 (MQ=255)                                cAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGAcc                                                                                                 <  1:414393/141‑1 (MQ=255)                                cAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGAcc                                                                                                 <  1:467918/141‑1 (MQ=255)                                   gACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACAAtt                                                                                              >  2:370520/1‑141 (MQ=255)                                    aacTGGTACACCAAAACCTGCGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTACTGTCTGATCAGATGTATCAGTGTGAATAGACCAttt                                                                                             >  2:190423/3‑141 (MQ=255)                                        ggTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAAAAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTAACGGTGTGAATAAACCATTTATGt                                                                                         >  1:153874/1‑141 (MQ=255)                                         gTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAGATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATGGACCATTTATGTa                                                                                        >  2:665871/1‑141 (MQ=255)                                         gTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTTTACTATCACGTGTATGCGTGTAAATGGACCATTTATCTa                                                                                        >  1:196464/1‑141 (MQ=255)                                         gTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTCTCAGTGTGAATTGACCAtttattta                                                                                        >  1:196324/1‑141 (MQ=255)                                           acacCAATTCCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAAt                                                                                      <  1:326550/141‑1 (MQ=255)                                           acacCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAAt                                                                                      <  2:48914/141‑1 (MQ=255)                                             acCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTg                                                                                    <  2:216461/141‑1 (MQ=255)                                             acCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTg                                                                                    <  1:650992/141‑1 (MQ=255)                                                        ctctaaGTACACTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAg                                                                         <  2:684655/137‑1 (MQ=255)                                                        acaaaaGTACATGTCAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTCTTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAg                                                                         <  2:351009/137‑1 (MQ=255)                                                           cAAGTACATTTGAAAAATTATCATTATTTTCTTTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGgtg                                                                      >  1:445369/1‑141 (MQ=255)                                                           cAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGgtg                                                                      >  1:51586/1‑141 (MQ=255)                                                             aGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCCGGtgtg                                                                    >  1:430621/1‑141 (MQ=255)                                                             aGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGtgtg                                                                    <  2:110751/141‑1 (MQ=255)                                                                  aTTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAAt                                                               >  2:241663/1‑141 (MQ=255)                                                                   tttGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATa                                                              <  1:578386/141‑1 (MQ=255)                                                                       aaaaaTTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTa                                                          >  2:566692/1‑141 (MQ=255)                                                                            ttataattatTTTCATTTACTTTTTTAACATACTAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAttttt                                                     >  1:313992/1‑140 (MQ=255)                                                                            ttataattatTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCAGTTATGTAATGGTTCTATTGCAGGTGTGAAAATTCGTATTTTa                                                     >  1:853692/1‑141 (MQ=255)                                                                              ataattatTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACa                                                   <  2:97866/141‑1 (MQ=255)                                                                                    attTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAgtcgt                                             <  2:18877/141‑1 (MQ=255)                                                                                      tttCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTcc                                           <  1:149854/141‑1 (MQ=255)                                                                                           tttACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAGTTTACAAGTCTTCCCGCCt                                      >  2:637874/1‑141 (MQ=255)                                                                                             tACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTCCAAGTAGTCCAGCCTAc                                    >  2:867716/1‑141 (MQ=255)                                                                                               cttcttCAACAACCGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGCCCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACAt                                  <  1:247793/141‑1 (MQ=255)                                                                                               cttcttCAACAAACGCATCATACATTCCTTGCAAAAGTGTATTTCCTTTACTATCAGATGTATCAGTGTGAATAGACCATTTATGCAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACAt                                  <  2:415283/141‑1 (MQ=255)                                                                                               cTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTATAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCTCGCCTACAt                                  <  2:592579/141‑1 (MQ=255)                                                                                                   ttCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAgc                              >  2:901019/1‑141 (MQ=255)                                                                                                    tCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATATCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAgat                             >  2:768848/1‑139 (MQ=255)                                                                                                       ttaaaCTAATCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAgctgct                          >  2:524794/3‑141 (MQ=255)                                                                                                              aaTCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCTGGAGAGAGAATACCTAAACTAAAAGTAGTCCCGCCTACTTCAGCTGCTAAaatat                   >  1:901874/1‑140 (MQ=255)                                                                                                              aaTCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGCAAATACCTAATTTACAAGTCGTCCCGCCTACATCAGCTGCTAAaataa                   >  1:113459/1‑141 (MQ=255)                                                                                                               aTCATAAATTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAGCTGCTAAaataat                  >  1:24367/1‑141 (MQ=255)                                                                                                                     aaTTCCTTTCGAAAGTGTTTTTCCGGTACTATCATATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAGCTGCTAAAATAATTTTTGt            >  2:914637/1‑141 (MQ=255)                                                                                                                     aaTTCCTTTCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCCCCTACATCAGCTGCTAAAATAATTTTTGt            >  2:409317/1‑141 (MQ=255)                                                                                                                           tttCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAGCTGCTAAAATAATTTTTGTCATTTc      <  1:454481/141‑1 (MQ=255)                                                                                                                            ttCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGTATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAGCTGCTAAAATAATTTTTGTCATTTCt     <  2:918046/141‑1 (MQ=255)                                                                                                                            ttCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAGCTGCTAAAATAATTTTTGTCATTTCt     <  1:241663/141‑1 (MQ=255)                                                                                                                            ttCAAAAGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAGCTGCTAAAATAATTTTTGTCATTTCt     <  2:671634/141‑1 (MQ=255)                                                                                                                               aaaaGTGTATTTCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGATCTAATTCAGGTGTGAAAATACCTAATTTACAAGTGGTCCCGCCTACCTCAGCTGCTATAATAATTTTTGGCCTTTTTGtt  >  2:513491/1‑141 (MQ=255)                                                                                                                                         |                                                                                                                                   GCTCCATTTACTGTACCTTTTTCAAAGTCAACAGGACCTGGTACACCAATACCTACGCCAAGTACATTTGAAAAATTATAATTATTTTCATTTACTTTTTCAACAAACGAATCATAAATTCCTTTCAAAAGTGTATATCCTGTACTATCAGATGTATCAGTGTGAATAGACCATTTATGTAATTGTTCTAATTCAGGTGTGAAAATACCTAATTTACAAGTCGTCCCGCCTACATCAGCTGCTAAAATAATTTTTGTCATTTCTGTT  >  USA300TCH1516_ALE/1659405‑1659671

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