Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 253,983 A→G 7.4% intergenic (+196/+166) 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,9830AG 7.4% 700.1 / 28.4 225intergenic (+196/+166)asbF/USA300HOU_RS011253‑dehydroshikimate dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (122/85);  new base G (9/8);  total (132/93)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 241 total aligned reads displayed.

CAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTC  >  USA300TCH1516_ALE/253893‑254073
                                                                                          |                                                                                           
cAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAaga                                                                                             <  1:1462479/91‑1 (MQ=255)
 aaGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAagag                                                                                            <  1:1004751/91‑1 (MQ=255)
  aGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGt                                                                                           >  1:1188467/1‑91 (MQ=255)
   gTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGta                                                                                          <  2:3643049/91‑1 (MQ=255)
   gTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGta                                                                                          >  1:2824060/1‑91 (MQ=255)
     cATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtata                                                                                        <  1:1044308/91‑1 (MQ=255)
     cATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtata                                                                                        <  1:1323414/91‑1 (MQ=255)
     cATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtata                                                                                        <  2:2246546/91‑1 (MQ=255)
       ttgtttgtAAAGTAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatata                                                                                      >  1:2972186/1‑91 (MQ=255)
         gtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAc                                                                                    <  2:2711482/91‑1 (MQ=255)
         gtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAc                                                                                    <  1:2589343/91‑1 (MQ=255)
         gtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAAACCTCACCATTGAAAAGAGTATATAAc                                                                                    <  1:1305334/91‑1 (MQ=255)
          tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                   <  1:3837455/91‑1 (MQ=255)
           ttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCt                                                                                  >  2:1387571/1‑91 (MQ=255)
            tgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtt                                                                                 <  2:2893625/91‑1 (MQ=255)
             gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCttt                                                                                >  1:397981/1‑91 (MQ=255)
             gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCttt                                                                                >  2:20883/1‑91 (MQ=255)
             gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCttt                                                                                >  1:1703507/1‑91 (MQ=255)
              tAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtttt                                                                               <  1:2833612/91‑1 (MQ=255)
              tAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtttt                                                                               <  1:2253928/91‑1 (MQ=255)
               aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTc                                                                              <  1:1020125/91‑1 (MQ=25)
                   aagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATa                                                                          >  1:648220/1‑91 (MQ=255)
                     ggTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGt                                                                        <  1:2663702/91‑1 (MQ=255)
                      gtgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGGTtatatactc                                                                                   <  2:1731893/79‑4 (MQ=38)
                      gtgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAg                                                                         >  1:2547144/1‑89 (MQ=255)
                       tgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAATGGTGa                                                                      <  2:2455712/91‑1 (MQ=21)
                        gtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAg                                                                     <  2:2363139/91‑1 (MQ=255)
                         tACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAgg                                                                    >  1:3404819/1‑91 (MQ=255)
                          aCTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGt                                                                   >  1:3667798/1‑91 (MQ=255)
                           cTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGta                                                                  >  1:2656230/1‑91 (MQ=255)
                            tttATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGtat                                                                 <  2:2179569/91‑1 (MQ=255)
                              tatataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGCGGTatat                                                               >  1:2292353/1‑91 (MQ=255)
                                 ataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGat                                                            >  1:2566829/1‑91 (MQ=255)
                                 ataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGat                                                            >  2:2548507/1‑91 (MQ=255)
                                  taAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGata                                                           >  2:313482/1‑91 (MQ=255)
                                   aaGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataa                                                          <  1:1738770/91‑1 (MQ=255)
                                    aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                         >  1:1569217/1‑91 (MQ=255)
                                    aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                         >  2:3879045/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)
                                       atataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaa                                                     <  2:2947628/91‑1 (MQ=255)
                                       atataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaa                                                     >  2:200594/1‑91 (MQ=255)
                                          taGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaa                                                  >  1:2937046/1‑91 (MQ=255)
                                           aGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaa                                                 >  1:2983349/1‑91 (MQ=255)
                                            gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                >  2:480543/1‑91 (MQ=255)
                                            gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                >  2:3745182/1‑91 (MQ=255)
                                            gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                >  1:2832194/1‑91 (MQ=255)
                                             cGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAg                                               <  1:1202849/91‑1 (MQ=255)
                                              gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGc                                              <  2:1017658/91‑1 (MQ=255)
                                              gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGc                                              <  1:1016277/91‑1 (MQ=255)
                                              gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGc                                              >  1:2417075/1‑91 (MQ=255)
                                                 ggTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtg                                           >  1:3150703/1‑91 (MQ=255)
                                                  gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                          <  2:301600/91‑1 (MQ=255)
                                                   tCATACCCATTCACAGTAACAATCCTCACCATTGAACAGACTATATACCCTTTTCAATAGTGAGGTACATGATAATAAAAAAAAGCCtgtt                                         >  2:3608081/1‑91 (MQ=255)
                                                   tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaa                                                 >  1:1926428/1‑83 (MQ=255)
                                                   tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaa                                                 <  2:1926428/83‑1 (MQ=255)
                                                    cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGt                                                                                           <  1:1699042/41‑1 (MQ=255)
                                                    cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGt                                                                                           >  2:1699042/1‑41 (MQ=255)
                                                    cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttg                                        >  2:2721941/1‑91 (MQ=255)
                                                       aCCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTca                                     <  2:3398171/91‑1 (MQ=255)
                                                       aCCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTca                                     >  1:1617674/1‑91 (MQ=255)
                                                       aCCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTca                                     >  1:3318778/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)
                                                              cacaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAATGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGt                              >  2:329562/1‑91 (MQ=255)
                                                               acaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTc                             >  2:1763986/1‑91 (MQ=255)
                                                               acaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTc                             >  2:562287/1‑91 (MQ=255)
                                                                caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                            <  1:3452079/91‑1 (MQ=255)
                                                                caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                            >  2:2262420/1‑91 (MQ=255)
                                                                caGTAACAAGCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                            <  2:3522648/91‑1 (MQ=255)
                                                                 aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                            <  1:3021137/90‑1 (MQ=255)
                                                                 aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                           >  1:1140625/1‑91 (MQ=255)
                                                                 aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                           >  2:604055/1‑91 (MQ=255)
                                                                  gTAACAATCCTCACCATTGAAAAGAGTATATAACCCTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCCTa                          >  2:357335/1‑91 (MQ=255)
                                                                    aaCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCCATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGa                        >  1:1893741/1‑91 (MQ=255)
                                                                    aaCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGa                        >  1:1894393/1‑91 (MQ=255)
                                                                      cAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGaca                      <  2:3404819/91‑1 (MQ=255)
                                                                       aaTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGacac                     <  2:1124685/91‑1 (MQ=255)
                                                                        aTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACg                    >  2:2537836/1‑91 (MQ=255)
                                                                        aTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACg                    >  1:3416620/1‑91 (MQ=255)
                                                                           cTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACa                 <  1:1375174/91‑1 (MQ=255)
                                                                             cactaTTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATa               >  1:2949565/5‑91 (MQ=21)
                                                                             cACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATa               >  2:1488937/1‑91 (MQ=255)
                                                                              actaTTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc              >  2:1433467/4‑91 (MQ=21)
                                                                              aCCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc              <  1:706906/91‑1 (MQ=255)
                                                                               ctaTTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACt             >  1:871295/3‑91 (MQ=21)
                                                                               ccATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACt             <  1:445324/91‑1 (MQ=255)
                                                                                  ttGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTa          >  1:1227909/1‑91 (MQ=255)
                                                                                   tGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaa         >  2:3273750/1‑91 (MQ=37)
                                                                                    gAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaaa        >  2:705806/1‑91 (MQ=255)
                                                                                    gAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaaa        >  1:3022101/1‑91 (MQ=255)
                                                                                    gAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaaa        >  2:1849594/1‑91 (MQ=255)
                                                                                      aaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAg       <  2:185454/91‑1 (MQ=255)
                                                                                       aaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACGTACTTTAAAGGt     <  1:2997203/91‑1 (MQ=255)
                                                                                        agagTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGtt    <  1:2469245/91‑1 (MQ=255)
                                                                                         gagTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGttt   <  2:3162124/91‑1 (MQ=255)
                                                                                         gagTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGttt   >  2:2554775/1‑91 (MQ=255)
                                                                                          agTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcac                                    >  1:363440/1‑57 (MQ=255)
                                                                                          |                                                                                           
CAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTC  >  USA300TCH1516_ALE/253893‑254073

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