Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 253,984 G→T 7.6% intergenic (+197/+165) 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,9840GT 7.6% 769.9 / 37.7 225intergenic (+197/+165)asbF/USA300HOU_RS011253‑dehydroshikimate dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (122/85);  new base T (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 = 9.74e-01

Only 100 of 242 total aligned reads displayed.

AAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCT  >  USA300TCH1516_ALE/253894‑254074
                                                                                          |                                                                                           
aaGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAagag                                                                                             <  1:1004751/91‑1 (MQ=255)
aaGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAagag                                                                                             >  2:3160492/1‑91 (MQ=255)
 aGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGt                                                                                            >  1:1188467/1‑91 (MQ=255)
  gTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGta                                                                                           <  2:3643049/91‑1 (MQ=255)
   tCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtat                                                                                          <  1:162087/91‑1 (MQ=255)
    cATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtata                                                                                         <  1:1044308/91‑1 (MQ=255)
      ttgtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatata                                                                                       >  1:253264/1‑91 (MQ=255)
        gtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAc                                                                                     <  2:2711482/91‑1 (MQ=255)
        gtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAc                                                                                     <  1:2589343/91‑1 (MQ=255)
         tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATAACCATTCACAGTAACAATCCTCACCATTTAAAAGAGTATATAAcc                                                                                    <  1:3837431/91‑1 (MQ=255)
          ttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCt                                                                                   <  1:296192/91‑1 (MQ=255)
          ttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCt                                                                                   >  2:1387571/1‑91 (MQ=255)
           tgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtt                                                                                  <  2:2893625/91‑1 (MQ=255)
            gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCttt                                                                                 >  1:1703415/1‑91 (MQ=255)
            gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCttt                                                                                 >  1:169661/1‑91 (MQ=255)
            gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCttt                                                                                 >  1:397981/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)
              aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                               >  1:13176/1‑91 (MQ=255)
              aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                               >  1:13012/1‑91 (MQ=255)
                agaatGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCaa                                                                             <  1:1332870/91‑1 (MQ=25)
                 gaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAAt                                                                            <  2:1525163/91‑1 (MQ=25)
                 gaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAAt                                                                            >  1:3009276/1‑91 (MQ=255)
                  aagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATa                                                                           >  1:648149/1‑91 (MQ=255)
                    ggTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGt                                                                         <  1:2663702/91‑1 (MQ=255)
                     gtgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGGTtatatactc                                                                                    <  2:1731893/79‑4 (MQ=38)
                     gtgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTg                                                                        >  1:2827063/1‑91 (MQ=255)
                      tgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAATGGTGa                                                                       <  2:2455712/91‑1 (MQ=21)
                       gtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAg                                                                      <  2:2363139/91‑1 (MQ=255)
                         aCTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGt                                                                    >  1:2669447/1‑91 (MQ=255)
                         aCTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGt                                                                    >  1:3667798/1‑91 (MQ=255)
                          cTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGta                                                                   >  1:2655875/1‑91 (MQ=255)
                           tttATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGtat                                                                  <  2:2179569/91‑1 (MQ=255)
                              atataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATg                                                               >  1:2287357/1‑91 (MQ=255)
                                ataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGat                                                             >  2:2548507/1‑91 (MQ=255)
                                ataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGat                                                             >  1:3357106/1‑91 (MQ=255)
                                  aaGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataa                                                           <  1:1738770/91‑1 (MQ=255)
                                  aaGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataa                                                           >  2:1762456/1‑91 (MQ=255)
                                   aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                          >  2:3879045/1‑91 (MQ=255)
                                   aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                          >  2:1408830/1‑91 (MQ=255)
                                   aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                          >  1:1569217/1‑91 (MQ=255)
                                   aGTATATAGCGATGGTCATACCCATACACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                          <  1:3606606/91‑1 (MQ=255)
                                    gTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataata                                                        <  1:2302880/91‑1 (MQ=255)
                                    gTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataata                                                        <  1:2143309/91‑1 (MQ=255)
                                      atataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaa                                                      >  2:200594/1‑91 (MQ=255)
                                      atataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaa                                                      <  2:2947628/91‑1 (MQ=255)
                                          aGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaa                                                  >  1:2982922/1‑91 (MQ=255)
                                           gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                 <  2:3376288/91‑1 (MQ=255)
                                           gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                 >  2:480543/1‑91 (MQ=255)
                                           gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                 >  2:3745182/1‑91 (MQ=255)
                                            cGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAg                                                <  1:1202849/91‑1 (MQ=255)
                                             gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGc                                               <  2:1017658/91‑1 (MQ=255)
                                             gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGc                                               <  2:109461/91‑1 (MQ=255)
                                             gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGc                                               <  1:1016277/91‑1 (MQ=255)
                                                  tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaa                                                  >  1:1926428/1‑83 (MQ=255)
                                                  tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaa                                                  <  2:1926428/83‑1 (MQ=255)
                                                   cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGt                                                                                            >  2:1699042/1‑41 (MQ=255)
                                                   cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttg                                         >  1:2561875/1‑91 (MQ=255)
                                                   cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttg                                         >  2:2721941/1‑91 (MQ=255)
                                                    aTACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttgt                                        <  2:3591365/91‑1 (MQ=255)
                                                    aTACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTAAGGTATATGATAATAAAAAAAAGCCtgttgt                                        <  2:2499087/91‑1 (MQ=255)
                                                      aCCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTca                                      >  1:3318778/1‑91 (MQ=255)
                                                       cccATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcac                                     >  1:2386753/1‑91 (MQ=255)
                                                        ccATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcaca                                    <  2:3009276/91‑1 (MQ=255)
                                                          aTTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACGTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAAt                                  <  2:3402761/91‑1 (MQ=255)
                                                           ttCACAGTAACAATCCTCACCATTGAAAAGCGTATATAAACTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATg                                 >  1:1617654/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)
                                                               caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                             <  1:1772591/91‑1 (MQ=255)
                                                               caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                             <  1:3452079/91‑1 (MQ=255)
                                                                aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                            >  2:3207046/1‑91 (MQ=255)
                                                                aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                            >  1:1140625/1‑91 (MQ=255)
                                                                aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                            >  2:604055/1‑91 (MQ=255)
                                                                   aaCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGa                         >  1:1894371/1‑91 (MQ=255)
                                                                       aTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACg                     >  1:789168/1‑91 (MQ=255)
                                                                       aTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACg                     >  1:3416523/1‑91 (MQ=255)
                                                                        tCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGa                    >  2:2261533/1‑91 (MQ=255)
                                                                          cTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACa                  >  1:2906168/1‑91 (MQ=255)
                                                                           tCACTATTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACAt                 >  2:2115872/1‑91 (MQ=21)
                                                                            cactaTTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATa                >  1:2949565/5‑91 (MQ=21)
                                                                            cACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATa                >  2:1488937/1‑91 (MQ=255)
                                                                             aCCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc               <  1:706906/91‑1 (MQ=255)
                                                                              ctaTTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACt              >  2:3545310/3‑91 (MQ=21)
                                                                              ccATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACt              <  1:445324/91‑1 (MQ=255)
                                                                                aTTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACttt            >  2:3499343/1‑91 (MQ=255)
                                                                                 ttGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTa           >  1:1227909/1‑91 (MQ=255)
                                                                                 ttGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTa           <  1:311665/91‑1 (MQ=255)
                                                                                  tGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaa          >  2:3274564/1‑91 (MQ=25)
                                                                                   gAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaaa         >  1:3022101/1‑91 (MQ=255)
                                                                                   gAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaaa         >  2:1848478/1‑91 (MQ=255)
                                                                                   gAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaaa         >  2:705806/1‑91 (MQ=255)
                                                                                    aaaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAg        >  1:2466619/1‑91 (MQ=255)
                                                                                     aaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAg        <  2:185454/91‑1 (MQ=255)
                                                                                      aaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACGTACTTTAAAGGt      <  1:2997203/91‑1 (MQ=255)
                                                                                        gagTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGttt    >  2:2554775/1‑91 (MQ=255)
                                                                                          gTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGttt    <  2:228683/89‑1 (MQ=255)
                                                                                          gTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCt  <  1:138052/91‑1 (MQ=255)
                                                                                          gTATATAACCTTTTCAATAGTGAGGTATATGAAAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCt  <  2:3054561/91‑1 (MQ=255)
                                                                                          |                                                                                           
AAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCT  >  USA300TCH1516_ALE/253894‑254074

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: