Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 253,991 A→C 7.4% intergenic (+204/+158) asbF → / ← USA300HOU_RS01125 3‑dehydroshikimate dehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE253,9910AC 7.4% 742.8 / 37.5 227intergenic (+204/+158)asbF/USA300HOU_RS011253‑dehydroshikimate dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (126/84);  new base C (9/8);  total (135/92)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 231 total aligned reads displayed.

GTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATAT  >  USA300TCH1516_ALE/253902‑254080
                                                                                         |                                                                                          
 tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATAACCATTCACAGTAACAATCCTCACCATTTAAAAGAGTATATAAcc                                                                                          <  1:3837431/91‑1 (MQ=255)
  ttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCt                                                                                         >  2:1387571/1‑91 (MQ=255)
  ttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCt                                                                                         <  1:296192/91‑1 (MQ=255)
   tgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtt                                                                                        <  2:2894416/91‑1 (MQ=255)
    gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCttt                                                                                       >  2:20883/1‑91 (MQ=255)
    gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCttt                                                                                       >  1:1703507/1‑91 (MQ=255)
    gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCttt                                                                                       >  1:169661/1‑91 (MQ=255)
    gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCttt                                                                                       >  1:1438027/1‑91 (MQ=255)
     tAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtttt                                                                                      <  1:2833612/91‑1 (MQ=255)
      aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTc                                                                                     <  1:1020057/91‑1 (MQ=37)
      aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                     >  1:13176/1‑91 (MQ=255)
      aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                     >  1:13012/1‑91 (MQ=255)
      aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                     >  1:2630029/1‑91 (MQ=255)
        agaatGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCaa                                                                                   <  1:1332870/91‑1 (MQ=25)
        agaagGTGTACTTTATATAAGTATATAGCGATGGTCGTACCCATTCACACTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCaa                                                                                   <  1:3355507/91‑1 (MQ=255)
         gaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAAt                                                                                  <  2:1525163/91‑1 (MQ=25)
         gaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAAt                                                                                  <  2:1677554/91‑1 (MQ=255)
          aagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATa                                                                                 >  1:648220/1‑91 (MQ=255)
          aagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATa                                                                                 >  1:2298354/1‑91 (MQ=255)
            ggTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGt                                                                               <  1:2663702/91‑1 (MQ=255)
             gtgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAg                                                                                <  2:2547144/89‑1 (MQ=255)
             gtgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAg                                                                                >  1:2547144/1‑89 (MQ=255)
              tgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGa                                                                             >  1:2583219/1‑91 (MQ=255)
                tACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAgg                                                                           >  2:1772591/1‑91 (MQ=255)
                 aCTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGt                                                                          >  2:2165677/1‑91 (MQ=255)
                 aCTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGt                                                                          >  1:3667798/1‑91 (MQ=255)
                  cTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGta                                                                         >  1:2656230/1‑91 (MQ=255)
                      atataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATg                                                                     >  1:2287357/1‑91 (MQ=255)
                        ataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGat                                                                   >  2:2548520/1‑91 (MQ=255)
                         taAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGata                                                                  >  2:313482/1‑91 (MQ=255)
                          aaGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataa                                                                 >  2:1762456/1‑91 (MQ=255)
                          aaGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataa                                                                 <  2:2726656/91‑1 (MQ=255)
                           aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                                >  2:1408830/1‑91 (MQ=255)
                           aGTATATAGCGATGGTCATACCCATACACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                                <  1:3606606/91‑1 (MQ=255)
                            gTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataata                                                              <  1:2143309/91‑1 (MQ=255)
                            gTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataata                                                              >  1:3318406/1‑91 (MQ=255)
                              atataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaa                                                            <  2:2947629/91‑1 (MQ=255)
                                 taGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaa                                                         >  1:2937046/1‑91 (MQ=255)
                                  aGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaa                                                        >  2:2424190/1‑91 (MQ=255)
                                   gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                       >  2:428940/1‑91 (MQ=255)
                                   gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                       >  2:480543/1‑91 (MQ=255)
                                   gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                       >  1:2832194/1‑91 (MQ=255)
                                    cGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAg                                                      <  1:1202849/91‑1 (MQ=255)
                                     gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGc                                                     <  1:1016109/91‑1 (MQ=255)
                                     gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGc                                                     >  1:2417075/1‑91 (MQ=255)
                                       tGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGCTAATAAAAAAAAGCCt                                                   >  2:3711530/1‑91 (MQ=255)
                                         gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                                 <  2:301600/91‑1 (MQ=255)
                                          tCATACCCATTCCCAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgtt                                                >  1:2994273/1‑91 (MQ=255)
                                          tCATACCCATTCACAGTAACAATCCTCACCATTGAACAGACTATATACCCTTTTCAATAGTGAGGTACATGATAATAAAAAAAAGCCtgtt                                                >  2:3608081/1‑91 (MQ=255)
                                          tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaa                                                        >  1:1926428/1‑83 (MQ=255)
                                          tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaa                                                        <  2:1926428/83‑1 (MQ=255)
                                           cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttg                                               >  2:3205153/1‑91 (MQ=255)
                                            aTACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTAAGGTATATGATAATAAAAAAAAGCCtgttgt                                              <  2:2499087/91‑1 (MQ=255)
                                              aCCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTca                                            >  1:1617674/1‑91 (MQ=255)
                                              aCCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTca                                            <  2:3398171/91‑1 (MQ=255)
                                                ccATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcaca                                          <  2:3871123/91‑1 (MQ=255)
                                                   ttCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATg                                       >  2:1789791/1‑91 (MQ=255)
                                                    tCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATgg                                      >  1:909059/1‑91 (MQ=255)
                                                     cacaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGt                                     >  1:2354023/1‑91 (MQ=255)
                                                     cacaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAATGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGt                                     >  2:329562/1‑91 (MQ=255)
                                                      acaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTc                                    <  1:1575459/91‑1 (MQ=255)
                                                       caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                   >  2:2262420/1‑91 (MQ=255)
                                                       caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                   <  1:1772591/91‑1 (MQ=255)
                                                       caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                   <  1:3452079/91‑1 (MQ=255)
                                                       caGTAACAAGCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                   <  2:3522648/91‑1 (MQ=255)
                                                        aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                   <  1:3021137/90‑1 (MQ=255)
                                                        aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                                  >  2:604055/1‑91 (MQ=255)
                                                        aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                                  >  1:744395/1‑91 (MQ=255)
                                                         gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                                 >  1:2227700/1‑91 (MQ=255)
                                                         gTAACAATCCTCACCATTGAAAAGAGTATATAACCCTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCCTa                                 >  2:357335/1‑91 (MQ=255)
                                                           aaCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGa                               >  1:1894371/1‑91 (MQ=255)
                                                             cAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGaca                             <  2:3404819/91‑1 (MQ=255)
                                                             cAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAGCCTGTTGTCACAATGGTCATAGacac                            <  1:3850542/91‑1 (MQ=255)
                                                              aaTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGacac                            <  2:1124685/91‑1 (MQ=255)
                                                               aTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACg                           >  1:3416523/1‑91 (MQ=255)
                                                                  cTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTAAAAAAAAAAAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACa                        <  2:2720530/91‑1 (MQ=255)
                                                                   tCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACAt                       >  2:572064/1‑91 (MQ=255)
                                                                    cactaTTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATa                      >  1:2949565/5‑91 (MQ=21)
                                                                    cACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATa                      >  2:1488937/1‑91 (MQ=255)
                                                                     aCCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc                     <  1:1048360/91‑1 (MQ=255)
                                                                      ctaTTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACt                    >  1:1021314/3‑91 (MQ=21)
                                                                      ctaTTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACt                    >  1:871295/3‑91 (MQ=21)
                                                                      ccATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACt                    <  1:445324/91‑1 (MQ=255)
                                                                         ttGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTa                 >  1:1227909/1‑91 (MQ=255)
                                                                         ttGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTa                 <  1:311665/91‑1 (MQ=255)
                                                                          tGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaa                >  2:3273750/1‑91 (MQ=37)
                                                                           gAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaaa               >  2:705806/1‑91 (MQ=255)
                                                                           gAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaaa               <  2:2630029/91‑1 (MQ=255)
                                                                           gAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaaa               >  2:1849594/1‑91 (MQ=255)
                                                                            aaaaGAGTATATGACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAg              >  1:3004483/1‑91 (MQ=255)
                                                                            aaaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCCGTTGTCACAATGGTCATAGACACGACATACTTTAAAg              >  2:1052403/1‑91 (MQ=255)
                                                                             aaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAg              <  2:185454/91‑1 (MQ=255)
                                                                              aaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGt            >  2:2360933/1‑91 (MQ=255)
                                                                               agagTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGtt           <  1:2469245/91‑1 (MQ=255)
                                                                                 agTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTc         <  1:1181460/91‑1 (MQ=255)
                                                                                    atataACCTTTTCAATAGTGAGGTATATGATCATAAAAAAAAGCCTGTTGTCACAATCGTCATAGACACGACATACTTTAAAGGTTTCTGa      >  1:3450158/1‑91 (MQ=255)
                                                                                     tataACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaa     >  1:555166/1‑91 (MQ=255)
                                                                                      ataACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaat    <  2:2298354/91‑1 (MQ=255)
                                                                                        aaCCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatat  <  1:2802328/91‑1 (MQ=255)
                                                                                        aaCCTTTTCAATAGTGAGGTATATAATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatat  <  2:2033104/91‑1 (MQ=255)
                                                                                         |                                                                                          
GTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATAT  >  USA300TCH1516_ALE/253902‑254080

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