Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 171,052 A→T 5.9% intergenic (+401/‑49) USA300HOU_RS07235 → / → adhE hypothetical protein/Aldehyde‑alcohol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE171,0520AT 5.9% 199.7 / 3.0 68intergenic (+401/‑49)USA300HOU_RS07235/adhEhypothetical protein/Aldehyde‑alcohol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (26/38);  new base T (2/2);  total (28/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

GATTATAAAATCATTCTAAAAACACCGAAATAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTATGATTGATGCTCTAGCTGACAAAGGGAAAAA  >  USA300TCH1516_ALE/170913‑171189
                                                                                                                                           |                                                                                                                                         
gATTATAAAATCATTCTAAAAACACCGAAATAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTAc                                                                                                                                          >  1:566792/1‑141 (MQ=255)
aaTTCTAAAATCATTCTAAAAACACCGAAATAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTAc                                                                                                                                          <  1:260916/140‑1 (MQ=255)
   tataAAATCATTCTAAAAACACCGAAATAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACtat                                                                                                                                       <  1:577922/141‑1 (MQ=255)
      aaaaTCATTCTAAAAACACCGAAATAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAAt                                                                                                                                    <  1:605369/141‑1 (MQ=255)
       aaaTCATTCTAAAAACACCGAAATAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTTCTATAATg                                                                                                                                   >  1:423011/1‑141 (MQ=255)
       aaaTCATTCTAAAAACACCGAAATAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTTCTATAATg                                                                                                                                   >  1:70533/1‑141 (MQ=255)
       aaaTCATGCTAAAAACACCGAAATAACAATGATTTCATGAAAACATATATTTTAAAATTTGATATTTGTTCAAATAATATTCGAACTTAAACTTTTTTGTATAGAATTTTCTTTATATCTTTAGAGACATGTACtatactt                                                                                                                                   >  2:151205/1‑138 (MQ=255)
           cATTCTAAAAACACCGAAATAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTtg                                                                                                                               >  2:728923/1‑141 (MQ=255)
              tCTAAAAACACAGAAATAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGa                                                                                                                            <  2:702974/141‑1 (MQ=255)
               cTAAAAAGAACAAAAGAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGaa                                                                                                                           <  2:21910/141‑1 (MQ=255)
                   aaaCACCGAAATAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTTCTATAATGTTTGTGAaataa                                                                                                                       <  2:70533/141‑1 (MQ=255)
                   aaaCACCGAAATAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTTCTATAATGTTTGTGAaataa                                                                                                                       <  2:423011/141‑1 (MQ=255)
                   aaaCACCGAAATAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAaataa                                                                                                                       <  2:60042/141‑1 (MQ=255)
                      tacCGAAATAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTc                                                                                                                    <  2:708811/140‑1 (MQ=255)
                      catcGAAATAATAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTc                                                                                                                    <  1:722052/138‑1 (MQ=255)
                      cacCGAAATAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATCCACATGGAAGCGAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTc                                                                                                                    <  1:775808/141‑1 (MQ=255)
                        ccGAAATAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTcac                                                                                                                  <  1:183690/141‑1 (MQ=255)
                           aaaTAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTTAAAATTAAACTTTTTTTTATAGAATTTTTTTTATATAATTAGAGATATGTACTATAATATTTCTGAAATAATTCACaaa                                                                                                               >  2:670239/1‑141 (MQ=255)
                             aTAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGt                                                                                                             >  1:628246/1‑141 (MQ=255)
                             aTAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGt                                                                                                             >  1:131110/1‑141 (MQ=255)
                             aTAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGt                                                                                                             >  2:289431/1‑141 (MQ=255)
                                aCAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATCATGTTTGTGACATAATTCACCAAGtatt                                                                                                          >  1:712414/1‑140 (MQ=255)
                                aCAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGtata                                                                                                          <  2:475725/141‑1 (MQ=255)
                                    tGATTTCATGAACACCTTTATTGTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAgg                                                                                                      <  1:104796/141‑1 (MQ=255)
                                    tGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAACgg                                                                                                      >  1:119640/1‑141 (MQ=255)
                                    tGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAgg                                                                                                      >  1:500573/1‑141 (MQ=255)
                                    tGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAgg                                                                                                      <  1:444358/141‑1 (MQ=255)
                                    tGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAgg                                                                                                      <  1:59580/141‑1 (MQ=255)
                                    tGATTTCAGTAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCTTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAgg                                                                                                      <  2:329087/141‑1 (MQ=255)
                                    tGAGTTTATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAgg                                                                                                      <  1:272787/141‑1 (MQ=255)
                                     taTTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGa                                                                                                     <  2:225952/140‑1 (MQ=255)
                                      aTTTCATGAAAACATTTATTTTAAAATTTGATAGTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAg                                                                                                    <  1:341059/141‑1 (MQ=255)
                                          taTGAAAATATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGt                                                                                                <  2:193915/140‑1 (MQ=255)
                                          cATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGt                                                                                                <  1:335449/141‑1 (MQ=255)
                                           aTGAAAGCATGTATTTTAAAATTTGGTATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGtt                                                                                               <  2:224371/141‑1 (MQ=255)
                                                      tatttTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTaa                                                                                    >  1:154058/1‑141 (MQ=255)
                                                                     tatTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGCAGTGGTTGTATATGTTCACTATACCTGAAAAAg                                                                     >  1:341836/1‑141 (MQ=255)
                                                                         tGTTCAAATAATATTCGAAATTAAACTTTCTTGTATATAATTTTCTTAATAACTTGACAATCATGTACTCTAATCTTTCTTAAATCTTTCACAACGTATAAACGACCACTTGTATATGTTAACTATCCCTGaaaaaccaaa                                                                 >  2:320543/1‑136 (MQ=255)
                                                                                ataataTTAGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGa                                                          <  2:771772/141‑1 (MQ=255)
                                                                                 taataTTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATACTGAGAGATATGTACTATAATGTTTGTGAAATCATTCGCAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGAt                                                         <  1:112737/141‑1 (MQ=255)
                                                                                 taataTTCGAAAGTAAACTTTTTTGGATAGAATTTTCTTTATAACCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGAt                                                         <  1:309345/141‑1 (MQ=255)
                                                                                 taataTACGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGAt                                                         <  2:233816/141‑1 (MQ=255)
                                                                                    tatTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAACAAAGTGGAGCGa                                                      >  2:271364/1‑141 (MQ=255)
                                                                                     atTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGaa                                                     >  1:593931/1‑141 (MQ=255)
                                                                                       tCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGAATATGTAAACTATACCTGAAAAAGAAAATCGTGGATCGgaaa                                                   >  2:120622/1‑140 (MQ=255)
                                                                                         gAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGaa                                                 <  1:712229/141‑1 (MQ=255)
                                                                                         gAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGaa                                                 <  1:53720/141‑1 (MQ=255)
                                                                                           aagtAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACa                                               <  2:206880/138‑1 (MQ=255)
                                                                                           aaTTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACa                                               >  1:640740/1‑141 (MQ=255)
                                                                                               aaaCTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACaagaa                                           >  1:59343/1‑141 (MQ=255)
                                                                                                aaCTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACaagaag                                          <  1:128552/141‑1 (MQ=255)
                                                                                                 aCTTTGTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGt                                         <  1:393195/141‑1 (MQ=255)
                                                                                                    tttttGTATAGAATTGTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGc                                      <  2:628122/141‑1 (MQ=255)
                                                                                                      tttGTCTATAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCaa                                    <  1:289431/141‑1 (MQ=255)
                                                                                                      tttGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCaa                                    >  1:235557/1‑141 (MQ=255)
                                                                                                        tGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAAtt                                  <  1:320756/141‑1 (MQ=255)
                                                                                                         gTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTa                                 <  1:728923/141‑1 (MQ=255)
                                                                                                            taGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTAtga                              >  1:277692/1‑141 (MQ=255)
                                                                                                            taGAATTTTCTTTATAGCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTAtga                              <  2:696554/141‑1 (MQ=255)
                                                                                                               aaTTTTCTTTATATACTGAGAGAGATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTAtgattg                           <  2:423728/141‑1 (MQ=255)
                                                                                                                 ttttCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTAtgattgat                         <  1:623983/141‑1 (MQ=255)
                                                                                                                          tatCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTATGATTGATGCTCTAGCt                >  2:226976/1‑141 (MQ=255)
                                                                                                                            tCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTATGATTGATGCTCTAGCTGa              >  1:685903/1‑141 (MQ=255)
                                                                                                                               tGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTATGATTGATGCTCTAGCTGACaa           >  1:386568/1‑141 (MQ=255)
                                                                                                                                   agaCATGGAATATAATGTTTGTGAAATAGTTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTATGATTGATGCTCTAGCTGGCAAAggg       <  2:326164/141‑1 (MQ=255)
                                                                                                                                    gaCATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTATGATTGATGCTCTAGCTGACAAAGGGa      >  1:159222/1‑141 (MQ=255)
                                                                                                                                        tGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTATGATTGATGCTCTAGCTGACAAAGGGaaaaa  >  1:341289/1‑141 (MQ=255)
                                                                                                                                        tGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTATGATTGATGCTCTAGCTGACAAAGGGaaaaa  >  2:151830/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                         
GATTATAAAATCATTCTAAAAACACCGAAATAACAATGATTTCATGAAAACATTTATTTTAAAATTTGATATTTGTTCAAATAATATTCGAAATTAAACTTTTTTGTATAGAATTTTCTTTATATCCTGAGAGACATGTACTATAATGTTTGTGAAATAATTCACAAAGTATAAAGGAGTGGTTGTATATGTTAACTATACCTGAAAAAGAAAATCGTGGATCGAAAGAACAAGAAGTGGCAATTATGATTGATGCTCTAGCTGACAAAGGGAAAAA  >  USA300TCH1516_ALE/170913‑171189

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: