Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 254,002 A→G 5.1% intergenic (+215/+147) 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_ALE254,0020AG 5.1% 725.6 / 18.7 225intergenic (+215/+147)asbF/USA300HOU_RS011253‑dehydroshikimate dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (126/87);  new base G (9/3);  total (135/90)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 231 total aligned reads displayed.

AAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAA  >  USA300TCH1516_ALE/253912‑254092
                                                                                          |                                                                                           
aagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATa                                                                                             >  1:1274538/1‑91 (MQ=255)
aagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATa                                                                                             >  1:648220/1‑91 (MQ=255)
    tgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAATGGTGa                                                                                         <  2:2455712/91‑1 (MQ=21)
    tgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGa                                                                                         >  1:2583219/1‑91 (MQ=255)
     gtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAg                                                                                        <  2:2363139/91‑1 (MQ=255)
      tACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAATGGTGAgg                                                                                       <  2:583607/91‑1 (MQ=18)
      tACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAgg                                                                                       >  1:3404819/1‑91 (MQ=255)
       aCTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGt                                                                                      >  2:2165677/1‑91 (MQ=255)
        cTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGta                                                                                     >  1:2666582/1‑91 (MQ=255)
        cTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGta                                                                                     >  1:2656230/1‑91 (MQ=255)
           tatataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGCGGTatat                                                                                  >  1:2292353/1‑91 (MQ=255)
              ataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGat                                                                               >  2:2548507/1‑91 (MQ=255)
              ataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATAat                                                                               <  2:1438027/91‑1 (MQ=255)
               taAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGata                                                                              >  2:313482/1‑91 (MQ=255)
                aaGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataa                                                                             <  2:2726656/91‑1 (MQ=255)
                aaGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataa                                                                             <  1:1738770/91‑1 (MQ=255)
                 aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                                            >  2:1408830/1‑91 (MQ=255)
                 aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                                            >  2:1494871/1‑91 (MQ=255)
                  gTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataata                                                                           <  1:2143309/91‑1 (MQ=255)
                  gTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataata                                                                           <  1:2302880/91‑1 (MQ=255)
                   tatataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGAtaataa                                                                          >  2:1607298/1‑91 (MQ=255)
                       taGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaa                                                                      >  1:2937046/1‑91 (MQ=255)
                        aGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaa                                                                     >  1:2982922/1‑91 (MQ=255)
                        aGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaa                                                                     >  2:1737641/1‑91 (MQ=255)
                         gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                                    >  2:428940/1‑91 (MQ=255)
                         gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                                    <  2:616675/91‑1 (MQ=255)
                         gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                                    >  1:2832194/1‑91 (MQ=255)
                         gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                                    >  2:480543/1‑91 (MQ=255)
                         gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                                    <  2:3376288/91‑1 (MQ=255)
                          cGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAg                                                                   <  1:1202849/91‑1 (MQ=255)
                           gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGc                                                                  <  2:1274183/91‑1 (MQ=255)
                           gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGc                                                                  <  1:1016109/91‑1 (MQ=255)
                           gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGc                                                                  <  1:2889058/91‑1 (MQ=255)
                           gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGc                                                                  <  2:1017658/91‑1 (MQ=255)
                           gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGc                                                                  <  2:109461/91‑1 (MQ=255)
                             tGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGCTAATAAAAAAAAGCCt                                                                >  2:3711530/1‑91 (MQ=255)
                                tCATACCCATTCCCAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgtt                                                             >  1:2994273/1‑91 (MQ=255)
                                tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaa                                                                     >  1:1926428/1‑83 (MQ=255)
                                tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaa                                                                     <  2:1926428/83‑1 (MQ=255)
                                 cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttg                                                            >  1:2561875/1‑91 (MQ=255)
                                  aTACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttgt                                                           <  1:3667626/91‑1 (MQ=255)
                                  aTACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttgt                                                           >  2:1853300/1‑91 (MQ=255)
                                    aCCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTca                                                         >  1:3318778/1‑91 (MQ=255)
                                    aCCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTca                                                         <  2:3398171/91‑1 (MQ=255)
                                      ccATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcaca                                                       <  2:3871123/91‑1 (MQ=255)
                                      ccATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcaca                                                       <  2:3009276/91‑1 (MQ=255)
                                      ccATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcaca                                                       <  2:3536536/91‑1 (MQ=255)
                                         ttCACAGTAACAATCCTCACCATTGAAAAGCGTATATAAACTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATg                                                    >  1:1617654/1‑91 (MQ=255)
                                         ttCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATg                                                    >  1:1632726/1‑91 (MQ=255)
                                         ttCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATg                                                    >  2:1789791/1‑91 (MQ=255)
                                          tCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATgg                                                   >  1:909059/1‑91 (MQ=255)
                                            acaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTc                                                 >  2:562287/1‑91 (MQ=255)
                                            acaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTc                                                 >  2:1763986/1‑91 (MQ=255)
                                             caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                                <  1:3452079/91‑1 (MQ=255)
                                              aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                                >  2:3021137/1‑90 (MQ=255)
                                              aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                                               >  2:3207046/1‑91 (MQ=255)
                                              aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                                               >  1:1140625/1‑91 (MQ=255)
                                               gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                                              <  2:1975122/91‑1 (MQ=255)
                                               gTAACAATCCTCACCATTGAAAAGAGTATATAACCCTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCCTa                                              >  2:357335/1‑91 (MQ=255)
                                                 aaCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCCATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGa                                            >  1:1893741/1‑91 (MQ=255)
                                                 aaCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGa                                            >  1:1894393/1‑91 (MQ=255)
                                                   cAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAGCCTGTTGTCACAATGGTCATAGacac                                         <  1:3850542/91‑1 (MQ=255)
                                                     aTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACg                                        >  1:789168/1‑91 (MQ=255)
                                                     aTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACg                                        >  2:2537836/1‑91 (MQ=255)
                                                     aTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACg                                        >  1:3416620/1‑91 (MQ=255)
                                                        cTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACa                                     <  1:1375174/91‑1 (MQ=255)
                                                        cTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTAAAAAAAAAAAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACa                                     <  2:2720530/91‑1 (MQ=255)
                                                         tCACTATTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACAt                                    >  2:2115872/1‑91 (MQ=21)
                                                           aCCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc                                  <  1:706906/91‑1 (MQ=255)
                                                           aCCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc                                  <  1:1048360/91‑1 (MQ=255)
                                                            ctaTTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACt                                 >  1:1021314/3‑91 (MQ=21)
                                                            ccATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACt                                 <  1:445324/91‑1 (MQ=255)
                                                              aTTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACttt                               >  2:3499343/1‑91 (MQ=255)
                                                               ttGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTa                              <  1:311665/91‑1 (MQ=255)
                                                               ttGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTa                              >  1:1227909/1‑91 (MQ=255)
                                                                tGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaa                             >  2:3274564/1‑91 (MQ=25)
                                                                 gAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaaa                            >  2:1848478/1‑91 (MQ=255)
                                                                 gAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaaa                            >  1:3022101/1‑91 (MQ=255)
                                                                  aaaaGAGTATATGACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAg                           >  1:3004483/1‑91 (MQ=255)
                                                                    aaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGt                         >  2:2360566/1‑91 (MQ=255)
                                                                     agagTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGtt                        <  1:2469245/91‑1 (MQ=255)
                                                                      gagTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGttt                       >  2:2554775/1‑91 (MQ=255)
                                                                       agTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcac                                                        >  1:363440/1‑57 (MQ=255)
                                                                       agTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcac                                                        <  2:363440/57‑1 (MQ=255)
                                                                       agTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTc                      <  1:1181460/91‑1 (MQ=255)
                                                                        gTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGttt                       >  1:228683/1‑89 (MQ=255)
                                                                        gTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGttt                       <  2:228683/89‑1 (MQ=255)
                                                                           tataACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaa                 >  1:555166/1‑91 (MQ=255)
                                                                             taACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaata               >  1:2395618/1‑91 (MQ=255)
                                                                              aaCCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatat              <  1:2802328/91‑1 (MQ=255)
                                                                              aaCCTTTTCAATAGTGAGGTATATGAAAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatat              <  2:3785703/91‑1 (MQ=255)
                                                                              aaCCTTTTCAATAGTGAGGTATATAATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatat              <  2:2033104/91‑1 (MQ=255)
                                                                                  ttttCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatataata          <  2:2827063/91‑1 (MQ=255)
                                                                                  ttttCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatataata          <  1:3195664/91‑1 (MQ=255)
                                                                                    ttCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATAtt        >  1:2712388/1‑91 (MQ=255)
                                                                                     tCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATAttt       >  1:2436432/1‑91 (MQ=255)
                                                                                     tCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATAttt       >  1:2691873/1‑91 (MQ=255)
                                                                                        aTAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAAatattatattcag     <  1:3492648/91‑4 (MQ=255)
                                                                                        aTAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAg    >  2:2186251/1‑91 (MQ=255)
                                                                                         tAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGa   <  1:139384/91‑1 (MQ=255)
                                                                                          |                                                                                           
AAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAA  >  USA300TCH1516_ALE/253912‑254092

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: