Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,353,873 A→C 5.1% intergenic (+66/+77) nucH → / ← USA300HOU_RS06735 Thermonuclease/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,353,8730AC 5.1% 368.2 / 11.7 119intergenic (+66/+77)nucH/USA300HOU_RS06735Thermonuclease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (49/64);  new base C (4/2);  total (53/66)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

Only 100 of 122 total aligned reads displayed.

CAAATATAGAAATGTATTTATAGAAGCACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCTTTATTTAAAGGGTGTTTAGATAGCTCT  >  USA300TCH1516_ALE/1353735‑1354013
                                                                                                                                          |                                                                                                                                                
cAAATATAGAAATGTATTTATAGAAGCACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAAtt                                                                                                                                                <  1:365458/141‑1 (MQ=255)
 aaaTATAGAAATGTATTTATAGAAGCACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAAttt                                                                                                                                               >  1:1009347/1‑141 (MQ=255)
 aaaTATAGAAATGTATTTATAGAAGCACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAAttt                                                                                                                                               >  1:504823/1‑141 (MQ=255)
  aaTATAGAAATGTATTTATAGAAGCACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTta                                                                                                                                              >  3:132319/1‑141 (MQ=255)
  aaTATAGAAATGTATTTATAGAAGCACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACAATAAATTta                                                                                                                                              <  2:223009/141‑1 (MQ=255)
    tataGAAATGTATTTATAGAAGCACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTtata                                                                                                                                            <  1:286217/141‑1 (MQ=255)
      taGAAATGTATTTATAGAAGCACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTtatata                                                                                                                                          >  1:176439/1‑141 (MQ=255)
          aaTGTATTTATAGAAGCACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑Ga                                                                                                                                  >  3:29204/1‑141 (MQ=255)
          aaTGTATTTATAGAAGCACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑Ga                                                                                                                                  >  1:659076/1‑141 (MQ=255)
            tGTATTTATAGAAGCACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑Gaaa                                                                                                                                >  2:1051979/1‑141 (MQ=255)
             gTATTTATAGAAGTACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAaat                                                                                                                               <  1:46158/141‑1 (MQ=255)
             gTATTTATAGAAGCACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAaat                                                                                                                               <  1:137512/141‑1 (MQ=255)
                 ttATAGAAGCACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAAtt                                                                                                                           <  2:1045281/141‑1 (MQ=255)
                 ttATAGAAGCACAAAATAAAGCTAAACAACAAAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAAtt                                                                                                                           <  3:329273/141‑1 (MQ=255)
                  tataGAAGCACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAACTTTATATATAAA‑‑GAaataat                                                                                                                            <  1:661770/141‑1 (MQ=255)
                    taGAAGCACAAAATAAAGCTAAAAAACGGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTaa                                                                                                                        >  2:1023827/1‑141 (MQ=255)
                        aGCACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGc                                                                                                                    <  2:117461/141‑1 (MQ=255)
                         gCATAAAATAAATCTAAACAACATAATTTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCa                                                                                                                   <  2:973169/141‑1 (MQ=255)
                            caAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTCAGTAGGTGATTATATTTATTAAAGTAAAATATTCTTGAAACATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCataa                                                                                                                >  1:541859/1‑141 (MQ=255)
                             aaaaTAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCataat                                                                                                               <  2:80399/141‑1 (MQ=255)
                             aaaaTAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCataat                                                                                                               <  1:731115/141‑1 (MQ=255)
                                aTAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAAc                                                                                                            >  2:540502/1‑141 (MQ=255)
                                 aaaaGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAAcc                                                                                                           <  2:602238/140‑1 (MQ=255)
                                  aaaGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCa                                                                                                          <  2:119071/141‑1 (MQ=255)
                                   aaGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCAt                                                                                                         <  3:50963/141‑1 (MQ=255)
                                    aGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTGTATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATa                                                                                                        >  1:904129/1‑141 (MQ=255)
                                    aGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATa                                                                                                        >  2:895686/1‑141 (MQ=255)
                                      cTAATCAACAGTAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATTATTTAAAGGCATACTAACCATAGt                                                                                                      <  2:463039/141‑1 (MQ=255)
                                      cTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGt                                                                                                      >  2:159398/1‑141 (MQ=255)
                                      cTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAAACATAGt                                                                                                      >  4:379175/1‑141 (MQ=255)
                                      atAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGt                                                                                                      <  2:541859/140‑1 (MQ=255)
                                        aaaCAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGc                                                                                                    >  2:44163/1‑141 (MQ=255)
                                              caGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCg                                                                                              <  2:906320/141‑1 (MQ=255)
                                               aGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGa                                                                                             >  1:975423/1‑141 (MQ=255)
                                                gAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAg                                                                                            <  4:40132/141‑1 (MQ=255)
                                                gAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAg                                                                                            <  3:153721/141‑1 (MQ=255)
                                                   aTTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGt                                                                                         <  3:88609/141‑1 (MQ=255)
                                                     tAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTaa                                                                                       <  2:823132/141‑1 (MQ=255)
                                                       aaTATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATg                                                                                     >  1:453621/1‑141 (MQ=255)
                                                        atatTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATgg                                                                                    >  3:419835/1‑141 (MQ=255)
                                                           tttGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATTTACTTGAATCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGatt                                                                                 <  1:701953/141‑1 (MQ=255)
                                                           tttGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGatt                                                                                 <  4:327277/141‑1 (MQ=255)
                                                           tttGGAGAAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATCTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGatt                                                                                 <  1:521158/141‑1 (MQ=255)
                                                             tGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGattat                                                                               <  2:266149/141‑1 (MQ=255)
                                                             tGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGattat                                                                               >  2:985984/1‑141 (MQ=255)
                                                              ggAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGattact                                                                              <  1:89733/141‑3 (MQ=255)
                                                               gAGTAAATAATTATATAAGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATGAGTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGAttatta                                                                             <  2:971863/141‑1 (MQ=255)
                                                               gAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGAttatta                                                                             >  1:825561/1‑141 (MQ=255)
                                                                 gTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATg                                                                           <  1:749907/141‑1 (MQ=255)
                                                                   aaaTAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAACTTTATATAT‑‑AACGACATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATg                                                                           >  1:48079/1‑141 (MQ=255)
                                                                   aaaTAATTAAATAGGTGATTATATATATTAAAGCAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGcc                                                                         <  1:483429/141‑1 (MQ=255)
                                                                     ataatTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCtt                                                                       <  1:59857/141‑1 (MQ=255)
                                                                       aatTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTg                                                                     >  2:471591/1‑141 (MQ=255)
                                                                       aagtAAATAGGTGATTATACATATTAAAGAAAAATATACTTGAATCATAGATTTATTATTAACATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTg                                                                     <  4:220519/138‑1 (MQ=255)
                                                                        atTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGa                                                                    >  2:550205/1‑141 (MQ=255)
                                                                          tAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGAtt                                                                  <  2:753641/141‑1 (MQ=255)
                                                                             aTAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTa                                                               >  2:122508/1‑141 (MQ=255)
                                                                               aGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACt                                                             >  3:202098/1‑141 (MQ=255)
                                                                                 gTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTg                                                           <  3:266142/141‑1 (MQ=255)
                                                                                     ttATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAata                                                       <  1:52872/141‑1 (MQ=255)
                                                                                       atatatatTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAAtata                                                     >  1:988495/1‑141 (MQ=255)
                                                                                        tatatatTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCAAACGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAg                                                    >  4:415012/1‑141 (MQ=255)
                                                                                         atatatTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGc                                                   <  4:23792/141‑1 (MQ=255)
                                                                                            tatTAAAGTAAAATATATTTAAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCtgt                                                <  2:127575/141‑1 (MQ=255)
                                                                                             atTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCtgtg                                               >  1:1054304/1‑141 (MQ=255)
                                                                                             atTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCtgtg                                               <  2:451716/141‑1 (MQ=255)
                                                                                               tAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGAt                                             <  4:253179/141‑1 (MQ=255)
                                                                                               tAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGAt                                             <  2:1020815/141‑1 (MQ=255)
                                                                                               tAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGAt                                             <  2:280022/141‑1 (MQ=255)
                                                                                                  aGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCATATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGAGTTg                                          <  4:263902/141‑1 (MQ=255)
                                                                                                    tAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGAt                                        >  2:718958/1‑141 (MQ=255)
                                                                                                    tAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGAt                                        >  2:719062/1‑141 (MQ=255)
                                                                                                    tAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGAt                                        >  1:895082/1‑141 (MQ=255)
                                                                                                     aaaaTATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATCATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGAACATTATGCCTCTGATTTCACTTGAATATAGCTGTGAGTTGATc                                       >  2:718907/1‑141 (MQ=255)
                                                                                                      aaaTATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCa                                      >  1:843101/1‑141 (MQ=255)
                                                                                                      aaaTATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCa                                      >  2:474832/1‑141 (MQ=255)
                                                                                                       aaTATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCaa                                     >  3:406377/1‑141 (MQ=255)
                                                                                                       aaTATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCaa                                     >  1:570951/1‑141 (MQ=255)
                                                                                                        atatACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATATCTGTGATTTGATCAAc                                    <  3:355653/141‑1 (MQ=255)
                                                                                                        atatACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAAc                                    <  2:285527/141‑1 (MQ=255)
                                                                                                          ataCTTGAAGCATCGATTTATTATTACCATAACTTTATATATAAA‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAAc                                    <  2:616737/141‑1 (MQ=255)
                                                                                                          ataCTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCt                                  <  3:305108/141‑1 (MQ=255)
                                                                                                          ataCTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGACTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCt                                  >  2:415313/1‑141 (MQ=255)
                                                                                                            aCTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGt                                <  1:619145/141‑1 (MQ=255)
                                                                                                             cTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGtt                               >  4:173131/1‑141 (MQ=255)
                                                                                                                  aGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCttta                          >  1:1022576/1‑141 (MQ=255)
                                                                                                                   gCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGATCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCtttat                         <  2:148935/141‑1 (MQ=255)
                                                                                                                   gCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCtttat                         >  4:159334/1‑141 (MQ=255)
                                                                                                                      tCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCtttattta                      <  4:266702/141‑1 (MQ=255)
                                                                                                                      tCGATTTATTATTACCATAAATTTATAAAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCtttattta                      <  2:937144/141‑1 (MQ=255)
                                                                                                                        gATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCTTTATTTaaa                    <  4:12037/141‑1 (MQ=255)
                                                                                                                          tttATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCTTTATTTAAAgg                  <  1:742880/141‑1 (MQ=255)
                                                                                                                          tttATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCTTTATTTAAAgg                  <  3:15475/141‑1 (MQ=255)
                                                                                                                           ttattattaCCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCTTTATTTAAAggg                 <  1:34716/141‑1 (MQ=255)
                                                                                                                           ttattattaCCATAAAATTATATAT‑‑A‑‑GAAATAAGTTAAAGACATAATAACCAAAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCTTTATTTAAAggg                 <  4:377335/141‑1 (MQ=255)
                                                                                                                                 ttaCCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCTTTATTTAAAGGGTGTTTa           >  1:340552/1‑141 (MQ=255)
                                                                                                                                   aCCATAACTTTATATATAAA‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCTTTATTTAAAGGGTGTTTa           >  2:959786/1‑141 (MQ=255)
                                                                                                                                     cATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCTTTATTTAAAGGGTGTTTAGATa       >  2:150561/1‑141 (MQ=255)
                                                                                                                                          aTTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCTTTATTTAAAGGGTGTTTAGATAGctct  >  2:905880/1‑141 (MQ=255)
                                                                                                                                          aTTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCTTTATTTAAAGGGTGTTTAGATAGctct  >  1:292975/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                                
CAAATATAGAAATGTATTTATAGAAGCACAAAATAAAGCTAAACAACAGAAATTAAATATTTGGAGTAAATAATTAAATAGGTGATTATATATATTAAAGTAAAATATACTTGAAGCATCGATTTATTATTACCATAAATTTATATAT‑‑A‑‑GAAATAATTTAAAGGCATAATAACCATAGTGCGTATCGAGCGTAATGGATTATTATGCCTTTGATTTTACTTGAATATAGCTGTGATTTGATCAACCTGTTCTTTATTTAAAGGGTGTTTAGATAGCTCT  >  USA300TCH1516_ALE/1353735‑1354013

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