Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 253,992 C→T 6.6% intergenic (+205/+157) 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,9920CT 6.6% 728.6 / 23.2 225intergenic (+205/+157)asbF/USA300HOU_RS011253‑dehydroshikimate dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (126/82);  new base T (8/7);  total (135/90)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 235 total aligned reads displayed.

GTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAA  >  USA300TCH1516_ALE/253902‑254082
                                                                                          |                                                                                          
gtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAc                                                                                            <  2:2711482/91‑1 (MQ=255)
gtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAc                                                                                            <  2:1399556/91‑1 (MQ=255)
 tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATAACCATTCACAGTAACAATCCTCACCATTTAAAAGAGTATATAAcc                                                                                           <  1:3837431/91‑1 (MQ=255)
  ttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCt                                                                                          <  1:296192/91‑1 (MQ=255)
   tgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtt                                                                                         <  2:2893625/91‑1 (MQ=255)
    gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCttt                                                                                        >  1:169661/1‑91 (MQ=255)
    gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCttt                                                                                        >  1:1438027/1‑91 (MQ=255)
    gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCttt                                                                                        >  1:397981/1‑91 (MQ=255)
     tAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtttt                                                                                       <  1:2253928/91‑1 (MQ=255)
     tAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtttt                                                                                       <  1:2833612/91‑1 (MQ=255)
      aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTc                                                                                      <  1:1020057/91‑1 (MQ=37)
      aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                      >  1:13012/1‑91 (MQ=255)
      aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                      >  1:2630029/1‑91 (MQ=255)
        agaagGTGTACTTTATATAAGTATATAGCGATGGTCGTACCCATTCACACTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCaa                                                                                    <  1:3355507/91‑1 (MQ=255)
          aagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATa                                                                                  >  1:2298354/1‑91 (MQ=255)
           agGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAg                                                                                 <  1:1329929/91‑1 (MQ=255)
              tgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGa                                                                              >  1:2583219/1‑91 (MQ=255)
               gtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAATGGTGAg                                                                             <  2:1798118/91‑1 (MQ=21)
               gtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAg                                                                             <  2:2363139/91‑1 (MQ=255)
                 aCTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGt                                                                           >  1:3667798/1‑91 (MQ=255)
                  cTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGta                                                                          >  1:2666582/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)
                          aaGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataa                                                                  <  1:1738770/91‑1 (MQ=255)
                          aaGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataa                                                                  <  2:2726656/91‑1 (MQ=255)
                           aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                                 >  2:1494871/1‑91 (MQ=255)
                            gTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataata                                                                <  1:2143309/91‑1 (MQ=255)
                            gTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataata                                                                >  1:3318406/1‑91 (MQ=255)
                             tatataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGAtaataa                                                               >  2:1607298/1‑91 (MQ=255)
                              atataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaa                                                              <  2:2947628/91‑1 (MQ=255)
                                  aGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaa                                                          >  1:2983349/1‑91 (MQ=255)
                                   gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                         >  2:480543/1‑91 (MQ=255)
                                   gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                         >  1:2832194/1‑91 (MQ=255)
                                   gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                         >  2:428940/1‑91 (MQ=255)
                                    cGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAg                                                        <  1:1202553/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)
                                         gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATCATAATAAAAAAAAGCCtgt                                                   <  2:3387836/91‑1 (MQ=255)
                                          tCATACCCATTCACAGTAACAATCCTCACCATTGAACAGACTATATACCCTTTTCAATAGTGAGGTACATGATAATAAAAAAAAGCCtgtt                                                  >  2:3608081/1‑91 (MQ=255)
                                          tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaa                                                          >  1:1926428/1‑83 (MQ=255)
                                          tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAACCAAAGCCtgtt                                                  >  2:2462055/1‑91 (MQ=255)
                                           cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttg                                                 <  2:2566829/91‑1 (MQ=255)
                                           cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttg                                                 >  1:2561875/1‑91 (MQ=255)
                                            aTACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttgt                                                >  2:1853300/1‑91 (MQ=255)
                                               cccATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcac                                             >  1:2386753/1‑91 (MQ=255)
                                                   ttCACAGTAACAATCCTCACCATTGAAAAGCGTATATAAACTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATg                                         >  1:1617654/1‑91 (MQ=255)
                                                   ttCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATg                                         >  1:1632726/1‑91 (MQ=255)
                                                     cacaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAATGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGt                                       >  2:329562/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)
                                                       caGTAACAAGCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                     <  2:3522648/91‑1 (MQ=255)
                                                        aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                     <  1:3021137/90‑1 (MQ=255)
                                                        aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                                    >  1:1140625/1‑91 (MQ=255)
                                                        aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                                    >  1:744395/1‑91 (MQ=255)
                                                        aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                                    >  2:604055/1‑91 (MQ=255)
                                                        aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                                    >  2:1714549/1‑91 (MQ=255)
                                                         gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                                   <  2:1975122/91‑1 (MQ=255)
                                                           aaCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGa                                 >  1:1894371/1‑91 (MQ=255)
                                                             cAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGaca                               <  2:3404819/91‑1 (MQ=255)
                                                              aaTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGacac                              <  2:1124685/91‑1 (MQ=255)
                                                               aTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACg                             >  1:3416523/1‑91 (MQ=255)
                                                               aTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACg                             >  2:2537836/1‑91 (MQ=255)
                                                                tCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGa                            >  2:2061822/1‑91 (MQ=255)
                                                                 ccTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGAc                           <  1:2658840/91‑1 (MQ=255)
                                                                  cTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACa                          >  1:2906168/1‑91 (MQ=255)
                                                                  cTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACa                          <  1:1375174/91‑1 (MQ=255)
                                                                   tCACTATTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACAt                         >  2:2115872/1‑91 (MQ=21)
                                                                    cactaTTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATa                        >  1:2949565/5‑91 (MQ=21)
                                                                     actaTTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc                       >  2:1433467/4‑91 (MQ=21)
                                                                     aCCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc                       <  1:706409/91‑1 (MQ=255)
                                                                      ctaTTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACt                      >  1:1021314/3‑91 (MQ=21)
                                                                      ctaTTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACt                      >  2:3545310/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)
                                                                            aaaaGAGTATATGACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAg                >  1:3004483/1‑91 (MQ=255)
                                                                            aaaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAg                >  1:935574/1‑91 (MQ=255)
                                                                            aaaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCCGTTGTCACAATGGTCATAGACACGACATACTTTAAAg                >  2:1052403/1‑91 (MQ=255)
                                                                              aaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGt              >  2:2360566/1‑91 (MQ=255)
                                                                               agagTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGtt             >  2:188240/1‑91 (MQ=255)
                                                                               agagTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGtt             <  1:2469245/91‑1 (MQ=255)
                                                                                gagTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGttt            <  2:3162124/91‑1 (MQ=255)
                                                                                 agTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcac                                             >  1:363440/1‑57 (MQ=255)
                                                                                 agTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTc           <  1:1180680/91‑1 (MQ=255)
                                                                                  gTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGttt            <  2:228683/89‑1 (MQ=255)
                                                                                  gTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCt          >  1:1271712/1‑91 (MQ=255)
                                                                                  gTATATAACCTTTTCAATAGTGAGGTATATGAAAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCt          <  2:3054561/91‑1 (MQ=255)
                                                                                    atataACCTTTTCAATAGTGAGGTATATGATCATAAAAAAAAGCCTGTTGTCACAATCGTCATAGACACGACATACTTTAAAGGTTTCTGa        >  1:3450158/1‑91 (MQ=255)
                                                                                     tataACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaa       >  1:555166/1‑91 (MQ=255)
                                                                                        aaCCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatat    <  2:1774186/91‑1 (MQ=255)
                                                                                        aaCCTTTTCAATAGTGAGGTATATGAAAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatat    <  2:3785703/91‑1 (MQ=255)
                                                                                          ccTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatataa  >  1:2719075/1‑91 (MQ=255)
                                                                                          ccTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatataa  >  1:1840307/1‑91 (MQ=255)
                                                                                          |                                                                                          
GTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAA  >  USA300TCH1516_ALE/253902‑254082

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: