Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,600,622:1 +A 8.8% intergenic (‑34/+833) 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,6221.A 8.8% 357.3 / 28.6 91intergenic (‑34/+833)dapF/USA300HOU_RS13395Diaminopimelate epimerase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (56/27);  new base A (2/6);  total (58/33)
Fisher's exact test for biased strand distribution p-value = 2.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

TATCGTCATATTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCT  >  USA300TCH1516_ALE/2600532‑2600712
                                                                                           |                                                                                          
tATCGTCATATTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCctct                                                                                             <  1:3163137/91‑1 (MQ=255)
tATCGTCATATTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCctct                                                                                             >  1:483364/1‑91 (MQ=255)
 aTCGTCATATTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTg                                                                                           >  1:453311/1‑91 (MQ=255)
 aTCGTCATATTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTg                                                                                           >  1:453312/1‑91 (MQ=255)
  tCGTCATATTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCctctag                                                                                           <  1:40562/91‑3 (MQ=255)
  tCGTCATATTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCctctag                                                                                           <  1:827884/91‑3 (MQ=255)
  tCGTCATATTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGt                                                                                          <  1:3269121/91‑1 (MQ=255)
  tCGTCATATTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGt                                                                                          <  1:3266885/91‑1 (MQ=255)
  tCGTCATATTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGt                                                                                          <  1:44213/91‑1 (MQ=255)
   cGTCATATTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTa                                                                                         >  1:3494603/1‑91 (MQ=255)
    gTCATATTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTaa                                                                                        >  1:2313981/1‑91 (MQ=255)
      cATATTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAAtt                                                                                      >  1:2993923/1‑91 (MQ=255)
       atatTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAAttt                                                                                     <  1:1630910/91‑1 (MQ=255)
        tatTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAAtttt                                                                                    >  1:511296/1‑91 (MQ=255)
        tatTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAAtttt                                                                                    >  1:808442/1‑91 (MQ=255)
         atTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTAGTAAttaa                                                                                    <  1:762870/91‑3 (MQ=255)
         atTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTAGTAAttaa                                                                                    <  1:764152/91‑3 (MQ=255)
           tCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTAGTAAttaaaa                                                                                  <  1:490601/91‑5 (MQ=255)
              ccGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTAc                                                                              <  1:1158282/91‑1 (MQ=255)
              ccGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTAc                                                                              <  1:282832/91‑1 (MQ=255)
              ccGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTAc                                                                              <  1:3811047/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:897184/1‑91 (MQ=255)
                               gAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAg                                                             >  1:618467/1‑91 (MQ=255)
                               gAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAg                                                             >  1:618500/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:807904/1‑91 (MQ=255)
                                      cTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGc                                                      >  1:2744839/1‑91 (MQ=255)
                                        ataACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAgttgtttgct                                                    <  1:1195910/91‑5 (MQ=255)
                                        ataACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAgttgtttgct                                                    <  1:3192970/91‑5 (MQ=255)
                                        ataACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGa                                                    >  1:2768277/1‑91 (MQ=255)
                                        ataACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGa                                                    <  1:423895/91‑1 (MQ=255)
                                        ataACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGa                                                    <  1:3713334/91‑1 (MQ=255)
                                         taACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGaa                                                   >  1:979408/1‑91 (MQ=255)
                                         taACCTGCCGATTCATACGACACTCCTTTAGATGTATCTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGaa                                                   >  1:3787064/1‑91 (MQ=255)
                                         taACCTGCCGATTCATACGACACTCCTTTAGATCTATATTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGaa                                                   >  1:3306753/1‑91 (MQ=255)
                                         taACCTGCCGATTCATACGACACTCCTTTAGATCTATATTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGaa                                                   >  1:3306752/1‑91 (MQ=255)
                                           aCCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGc                                                 >  1:3368342/1‑91 (MQ=255)
                                             cTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCaa                                               <  1:568836/91‑1 (MQ=255)
                                              tGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCaaa                                              >  1:3385713/1‑91 (MQ=255)
                                              tGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCaaa                                              >  1:3383429/1‑91 (MQ=255)
                                              tGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCaaa                                              >  1:3383138/1‑91 (MQ=255)
                                               gCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAaac                                             <  1:2544183/91‑1 (MQ=255)
                                               gCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAaac                                             >  1:3918368/1‑91 (MQ=255)
                                               gCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAaac                                             >  1:3918868/1‑91 (MQ=255)
                                                ccGATTCATACGACACTCCTTTAGATGTATTTTTAGATCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAaaca                                            >  1:3482844/1‑91 (MQ=255)
                                                   aTTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACt                                         >  1:466495/1‑91 (MQ=255)
                                                      cATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTa                                      >  1:1897254/1‑91 (MQ=255)
                                                      cATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTa                                      >  1:3676651/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:224037/1‑91 (MQ=255)
                                                           gACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCct                                 >  1:2209836/1‑91 (MQ=255)
                                                           gACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCct                                 >  1:2208242/1‑91 (MQ=255)
                                                                tCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTaaa                            >  1:2876558/1‑91 (MQ=255)
                                                                   tttAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCt                         >  1:2120959/1‑91 (MQ=255)
                                                                   tttAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCt                         <  1:266342/91‑1 (MQ=255)
                                                                   tttAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCt                         >  1:1324738/1‑91 (MQ=255)
                                                                   tttAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCt                         <  1:3693636/91‑1 (MQ=255)
                                                                     tAGATGTATTTTTAGAGCTTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGt                       >  1:2801558/1‑91 (MQ=255)
                                                                       gATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTaa                     >  1:36764/1‑91 (MQ=255)
                                                                       gATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTaa                     >  1:3080679/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)
                                                                                           |                                                                                          
TATCGTCATATTCCCCGAAGGATTATACTTAGAAAATTCTATAACCTGCCGATTCATACGACACTCCTTTAGATGTATTTTTAGAGCCTCTGTAATTTTAATTACTAGAGGCTCTTATGCAGTTGTTGCGAAGCAAACAACTGTATTCCTCTAAACCTGTAATTTTTTAGAAGATGAACCT  >  USA300TCH1516_ALE/2600532‑2600712

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