Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,600,630 T→A 5.9% intergenic (‑42/+825) dapF ← / ← USA300HOU_RS13395 Diaminopimelate epimerase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,600,6300TA 5.9% 294.2 / 8.0 84intergenic (‑42/+825)dapF/USA300HOU_RS13395Diaminopimelate epimerase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (54/25);  new base A (2/3);  total (56/28)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TATTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTC  >  USA300TCH1516_ALE/2600540‑2600720
                                                                                           |                                                                                          
tatTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAAtttt                                                                                            >  1:511296/1‑91 (MQ=255)
tatTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAAtttt                                                                                            >  1:808442/1‑91 (MQ=255)
      ccGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTAc                                                                                      <  1:1158282/91‑1 (MQ=255)
      ccGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTAc                                                                                      <  1:3811047/91‑1 (MQ=255)
      ccGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTAc                                                                                      <  1:282832/91‑1 (MQ=255)
         aaGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTAGTAATTAAAATTAcag                                                                                    <  1:3014363/91‑3 (MQ=38)
         aaGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTag                                                                                   <  1:1414598/91‑1 (MQ=255)
           ggATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTagag                                                                                 <  1:3741044/91‑1 (MQ=255)
            gATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTAGTAATTAAAATTACAGAgg                                                                                <  1:155293/91‑1 (MQ=37)
              ttATACTTAGATAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGct                                                                              <  1:3763611/91‑1 (MQ=255)
              ttATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGct                                                                              >  1:4020659/1‑91 (MQ=255)
              ttATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGct                                                                              >  1:4021208/1‑91 (MQ=255)
               tataCTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGctc                                                                             >  1:661668/1‑91 (MQ=255)
               tataCTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGctc                                                                             <  1:1373982/91‑1 (MQ=255)
                ataCTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTAGTAATTAAAATTACAGAGGctct                                                                            <  1:3616487/91‑1 (MQ=25)
                ataCTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGctct                                                                            >  1:984495/1‑91 (MQ=255)
                ataCTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGctct                                                                            <  1:2292376/91‑1 (MQ=255)
                ataCTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGACTCTGTAATTTTAATTACTAGAGGctct                                                                            <  1:987817/91‑1 (MQ=255)
                 taCTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCtt                                                                           >  1:1333478/1‑91 (MQ=255)
                  aCTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTa                                                                          >  1:1312470/1‑91 (MQ=255)
                     tAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGc                                                                       >  1:180170/1‑91 (MQ=255)
                     tAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGc                                                                       >  1:2293538/1‑91 (MQ=255)
                       gAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAg                                                                     >  1:3119359/1‑91 (MQ=255)
                       gAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAg                                                                     >  1:618467/1‑91 (MQ=255)
                       gAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAg                                                                     >  1:618500/1‑91 (MQ=255)
                       gAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAg                                                                     >  1:897184/1‑91 (MQ=255)
                         aaaTTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAgtt                                                                   >  1:71663/1‑91 (MQ=255)
                          aaTTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAgttg                                                                  >  1:39126/1‑91 (MQ=255)
                           aTTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAgttgt                                                                 >  1:3513419/1‑91 (MQ=255)
                             tCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGttgttg                                                               >  1:758334/1‑91 (MQ=255)
                              cTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGc                                                              >  1:2744839/1‑91 (MQ=255)
                              cTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGc                                                              >  1:807904/1‑91 (MQ=255)
                                ataACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAgttgtttgct                                                            <  1:3192970/91‑5 (MQ=255)
                                ataACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAgttgtttgct                                                            <  1:1195910/91‑5 (MQ=255)
                                ataACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGa                                                            <  1:3713334/91‑1 (MQ=255)
                                ataACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGa                                                            >  1:2768277/1‑91 (MQ=255)
                                ataACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGa                                                            <  1:423895/91‑1 (MQ=255)
                                 taACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGaa                                                           >  1:979408/1‑91 (MQ=255)
                                 taACCTGCCGATTCATACGACACTCCTTTAGATGTATCTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGaa                                                           >  1:3787064/1‑91 (MQ=255)
                                 taACCTGCCGATTCATACGACACTCCTTTAGATCTATATTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGaa                                                           >  1:3306752/1‑91 (MQ=255)
                                 taACCTGCCGATTCATACGACACTCCTTTAGATCTATATTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGaa                                                           >  1:3306753/1‑91 (MQ=255)
                                   aCCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGc                                                         >  1:3368342/1‑91 (MQ=255)
                                     cTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCaa                                                       <  1:568836/91‑1 (MQ=255)
                                      tGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCaaa                                                      >  1:3383138/1‑91 (MQ=255)
                                      tGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCaaa                                                      >  1:3383429/1‑91 (MQ=255)
                                      tGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCaaa                                                      >  1:3385713/1‑91 (MQ=255)
                                       gCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAaac                                                     <  1:2544183/91‑1 (MQ=255)
                                       gCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAaac                                                     >  1:3918868/1‑91 (MQ=255)
                                       gCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAaac                                                     >  1:3918368/1‑91 (MQ=255)
                                        ccGATTCATACGACACTCCTTTAGATGTATTTTTAGATCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAaaca                                                    >  1:3482844/1‑91 (MQ=255)
                                           aTTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACt                                                 >  1:466495/1‑91 (MQ=255)
                                              cATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTa                                              >  1:3676651/1‑91 (MQ=255)
                                              cATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTa                                              >  1:1897254/1‑91 (MQ=255)
                                               aTACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTAt                                             >  1:3095303/1‑91 (MQ=255)
                                               aTACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTAt                                             >  1:3586881/1‑91 (MQ=255)
                                                  cGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTcc                                          <  1:751400/91‑1 (MQ=255)
                                                   gACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCct                                         >  1:2953793/1‑91 (MQ=255)
                                                   gACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCct                                         >  1:2209836/1‑91 (MQ=255)
                                                   gACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCct                                         >  1:224037/1‑91 (MQ=255)
                                                   gACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCct                                         >  1:2208242/1‑91 (MQ=255)
                                                        tCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTaaa                                    >  1:2876558/1‑91 (MQ=255)
                                                           tttAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCt                                 <  1:266342/91‑1 (MQ=255)
                                                           tttAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCt                                 >  1:2120959/1‑91 (MQ=255)
                                                           tttAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCt                                 <  1:3693636/91‑1 (MQ=255)
                                                           tttAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCt                                 >  1:1324738/1‑91 (MQ=255)
                                                             tAGATGTATTTTTAGAGCTTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGt                               >  1:2801558/1‑91 (MQ=255)
                                                               gATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTaa                             >  1:3080679/1‑91 (MQ=255)
                                                               gATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTaa                             >  1:36764/1‑91 (MQ=255)
                                                                    aTTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAAttttt                        >  1:527057/1‑91 (MQ=255)
                                                                     tttttAGAGCCTCTGTAATTTTAATTACTAGATGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTTAACCTGTAAtttttt                       <  1:3533030/91‑1 (MQ=255)
                                                                       tttAGAGTCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTag                     <  1:2072238/91‑1 (MQ=255)
                                                                       tttAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTag                     <  1:834359/91‑1 (MQ=255)
                                                                       tttAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTag                     <  1:957372/91‑1 (MQ=255)
                                                                       tttAGAGCCACTGTAAGTTTAATTACTAGAGGCTCTTATGCAGTAGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATGTTTTag                     <  1:3620876/91‑1 (MQ=21)
                                                                           gagCCTGTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTagaaga                 >  1:3232209/1‑91 (MQ=255)
                                                                           gagCCTGTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTagaaga                 >  1:3232357/1‑91 (MQ=255)
                                                                           gagCCTATGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTagaaga                 >  1:3232208/1‑91 (MQ=255)
                                                                            agCCTCTAGTAATTAAAATTACAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGAt                >  1:871228/1‑91 (MQ=17)
                                                                            agCCTCTAGTAATTAAAATTACAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGAt                >  1:871243/1‑91 (MQ=17)
                                                                            agCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGAt                >  1:3757715/1‑91 (MQ=255)
                                                                              ccTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGa              >  1:2023772/1‑91 (MQ=255)
                                                                                tctGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAAc            >  1:3446565/1‑91 (MQ=255)
                                                                                  tGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCt          >  1:493858/1‑91 (MQ=255)
                                                                                  tGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCt          <  1:2007075/91‑1 (MQ=255)
                                                                                     tAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTaa        <  1:858219/91‑1 (MQ=255)
                                                                                       aTTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCACACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATa      >  1:3600778/1‑91 (MQ=255)
                                                                                        ttttAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCtagtag     <  1:799287/91‑1 (MQ=255)
                                                                                        ttttAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAg     <  1:1519126/91‑1 (MQ=255)
                                                                                         tttAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGc    <  1:3151163/91‑1 (MQ=255)
                                                                                          ttAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGct   <  1:2766482/91‑1 (MQ=255)
                                                                                          ttAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGct   >  1:3183388/1‑91 (MQ=255)
                                                                                           tAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGctc  <  1:1479732/91‑1 (MQ=255)
                                                                                           tAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGctc  <  1:1478356/91‑1 (MQ=255)
                                                                                           |                                                                                          
TATTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCTAATAGCTC  >  USA300TCH1516_ALE/2600540‑2600720

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