Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 253,992 C→T 7.2% 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 7.2% 779.9 / 26.8 241intergenic (+205/+157)asbF/USA300HOU_RS011253‑dehydroshikimate dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (125/98);  new base T (8/10);  total (133/108)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

GTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATA  >  USA300TCH1516_ALE/253902‑254081
                                                                                           |                                                                                         
 tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                          <  1:1129276/91‑1 (MQ=255)
 tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                          <  2:313742/91‑1 (MQ=255)
 tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                          <  2:3836560/91‑1 (MQ=255)
   tgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtt                                                                                        <  2:1091178/91‑1 (MQ=255)
   tgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtt                                                                                        <  2:536961/91‑1 (MQ=255)
   tgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtt                                                                                        <  2:3665497/91‑1 (MQ=255)
    gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCttt                                                                                       <  1:2609629/91‑1 (MQ=255)
     tAAAGAAGGTGTACTTTATAGAAGTATATAGCGATGGTCATACCCATTCACAGTAGCAATCCTCACTATTGAAAAGGTTATATACTCtttt                                                                                      <  2:1547909/91‑1 (MQ=25)
      aaaGGAGGTGTACTTTATATAACTATATAGCGATGGTCATACCCATCCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                     >  2:2048330/1‑91 (MQ=255)
      aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                     >  2:1783173/1‑91 (MQ=255)
      aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                     >  1:534285/1‑91 (MQ=255)
      aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                     <  2:3205917/91‑1 (MQ=255)
      aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                     <  1:2293851/91‑1 (MQ=255)
       aagaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCa                                                                                    >  1:2143534/1‑91 (MQ=255)
        agaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCaa                                                                                   <  2:3903558/91‑1 (MQ=25)
        agaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCaa                                                                                   <  1:1508950/91‑1 (MQ=255)
        agaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCaa                                                                                   <  1:3499799/91‑1 (MQ=255)
          aagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATa                                                                                 >  2:724955/1‑91 (MQ=255)
           agGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAg                                                                                <  1:3815048/91‑1 (MQ=255)
           agGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAg                                                                                <  1:1453870/91‑1 (MQ=255)
           agGTGTACTTTAGATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAAtgg                                                                                <  2:1635880/91‑3 (MQ=21)
            ggTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGt                                                                               >  1:692346/1‑91 (MQ=255)
            ggTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTAACAGTAACAATCCTCACCATTGAAACGAGTATATAACCTTTTCAATAGt                                                                               >  2:3006756/1‑91 (MQ=255)
             gtgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTg                                                                              >  1:3938669/1‑91 (MQ=255)
               gtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAg                                                                            <  2:3860964/91‑1 (MQ=255)
                tACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAgg                                                                           >  1:753780/1‑91 (MQ=255)
                  cTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGta                                                                         <  2:2958649/91‑1 (MQ=255)
                  acttATATAAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGt                                                                          >  1:1247476/3‑91 (MQ=255)
                      atataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATg                                                                     >  2:1258507/1‑91 (MQ=255)
                        ataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGat                                                                   >  2:2205965/1‑91 (MQ=255)
                         taAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGata                                                                  <  1:1340952/91‑1 (MQ=255)
                           aaGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataa                                                                 <  2:2769104/91‑1 (MQ=255)
                            aGTCTATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGgtaat                                                                <  1:3122341/91‑1 (MQ=255)
                            aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                                >  2:3067573/1‑91 (MQ=255)
                            aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                                >  1:3845178/1‑91 (MQ=255)
                            aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTAAATGataat                                                                >  2:3696684/1‑91 (MQ=255)
                                  taGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATTGTGAGGTATATGATAATaaaaaa                                                          <  1:400023/91‑1 (MQ=255)
                                  taGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaa                                                          >  1:1272383/1‑91 (MQ=255)
                                    gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                        >  2:1839859/1‑91 (MQ=255)
                                     cGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAg                                                       >  1:3268451/1‑91 (MQ=255)
                                        tGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCt                                                    >  1:1001048/1‑91 (MQ=255)
                                          gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                                  <  1:1803378/91‑1 (MQ=255)
                                          gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                                  <  1:1481900/91‑1 (MQ=255)
                                          gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                                  >  1:3212239/1‑91 (MQ=255)
                                          gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                                  >  1:3058752/1‑91 (MQ=255)
                                           tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgtt                                                 <  1:1982019/91‑1 (MQ=255)
                                            cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttg                                                >  1:996918/1‑91 (MQ=255)
                                            cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttg                                                >  1:1571862/1‑91 (MQ=255)
                                              tACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGCGGTATATGATAATAAAAAAAGCCTGTTGTca                                             >  2:474127/1‑91 (MQ=255)
                                              tACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTc                                              <  2:1170613/91‑1 (MQ=255)
                                                cccATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcac                                            <  1:2862612/91‑1 (MQ=255)
                                                   aTTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTAAATGATAATAAAAAAAAGCCTGTTGTCACAAt                                         <  2:2883067/91‑1 (MQ=255)
                                                    ttCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATg                                        <  1:2455469/91‑1 (MQ=255)
                                                      cacaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGt                                      <  1:3842050/91‑1 (MQ=255)
                                                      cacaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGt                                      <  2:1574105/91‑1 (MQ=255)
                                                      cacaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGt                                      >  2:1030736/1‑91 (MQ=255)
                                                        caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                    >  2:1059595/1‑91 (MQ=255)
                                                        caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                    >  2:225243/1‑91 (MQ=255)
                                                        caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                    <  1:675515/91‑1 (MQ=255)
                                                        caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTCCAATAGCGCGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGCCa                                    >  2:622840/1‑91 (MQ=255)
                                                        caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTGTTCAATAGTGCGGTACATGATAATAAAAAAAAGCCTGTTGTCACAATGGCCa                                    >  2:3819451/1‑91 (MQ=255)
                                                         aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                                   >  2:395607/1‑91 (MQ=255)
                                                         aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                                   >  1:3162596/1‑91 (MQ=255)
                                                          gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                                  >  1:884439/1‑91 (MQ=255)
                                                          gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                                  >  2:3542832/1‑91 (MQ=255)
                                                             aCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGac                               <  1:1659342/91‑1 (MQ=255)
                                                             aCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGac                               <  1:1626344/91‑1 (MQ=255)
                                                             aCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGac                               >  2:430866/1‑91 (MQ=255)
                                                             aCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGac                               >  1:767386/1‑91 (MQ=255)
                                                              tataccTCACTATTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGaca                              >  2:1931054/5‑91 (MQ=17)
                                                                  ccTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGAc                          >  1:1541584/1‑91 (MQ=255)
                                                                  ccTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGAc                          <  1:203218/91‑1 (MQ=255)
                                                                  ccTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGAc                          >  1:2831715/1‑91 (MQ=255)
                                                                   cTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACa                         >  1:2099861/1‑91 (MQ=255)
                                                                    tCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCACAGACACGACAt                        >  1:3575164/1‑91 (MQ=255)
                                                                     cACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATa                       >  1:187114/1‑91 (MQ=255)
                                                                      aCCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc                      >  1:2575786/1‑91 (MQ=255)
                                                                      aCCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc                      >  2:2312320/1‑91 (MQ=255)
                                                                       ccATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACt                     >  2:1657931/1‑91 (MQ=255)
                                                                       ccATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACt                     >  1:3464019/1‑91 (MQ=255)
                                                                            gAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaaa                >  2:1489330/1‑91 (MQ=37)
                                                                            gAAAAGAGTGTATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaaa                <  1:3790852/91‑1 (MQ=255)
                                                                             aaaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAg               <  2:812646/91‑1 (MQ=255)
                                                                             aaaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAg               >  1:1392686/1‑91 (MQ=255)
                                                                              aaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAgg              <  2:2426737/91‑1 (MQ=255)
                                                                              aaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAgg              >  2:2602955/1‑91 (MQ=255)
                                                                              aaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAgg              >  1:1134983/1‑91 (MQ=255)
                                                                                   gTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCt         >  2:2913722/1‑91 (MQ=255)
                                                                                   gTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCt         >  1:928638/1‑91 (MQ=255)
                                                                                     atataACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGa       >  2:2822941/1‑91 (MQ=255)
                                                                                      tataACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaa      >  2:660049/1‑91 (MQ=255)
                                                                                       gtaACCTTTTCAATAGTGAGGTATATGATAATAAAATAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaat     <  2:3910272/90‑1 (MQ=255)
                                                                                       ataACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaat     <  1:2901450/91‑1 (MQ=255)
                                                                                       ataACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaat     >  2:3466272/1‑91 (MQ=255)
                                                                                       ataACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaat     >  2:264506/1‑91 (MQ=255)
                                                                                        taACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaata    <  1:1719968/91‑1 (MQ=255)
                                                                                        taACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaata    <  2:3274291/91‑1 (MQ=255)
                                                                                        taACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaata    <  1:561870/91‑1 (MQ=255)
                                                                                          aCCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatata  >  2:1419835/1‑91 (MQ=255)
                                                                                          aCCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatata  <  2:1003716/91‑1 (MQ=255)
                                                                                           |                                                                                         
GTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATA  >  USA300TCH1516_ALE/253902‑254081

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