Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 168,343 G→A 7.1% intergenic (+310/+38) yfkN_1 → / ← USA300TCH1516_00148 Trifunctional nucleotide phosphoesterase protein YfkN/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE168,3430GA 7.1% 453.4 / 16.3 140intergenic (+310/+38)yfkN_1/USA300TCH1516_00148Trifunctional nucleotide phosphoesterase protein YfkN/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (68/62);  new base A (2/8);  total (70/70)
Fisher's exact test for biased strand distribution p-value = 9.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 145 total aligned reads displayed.

ATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTT  >  USA300TCH1516_ALE/168253‑168433
                                                                                          |                                                                                          
aTTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAg                                                                                            >  1:3078209/1‑91 (MQ=255)
 tttAAAACTCTTTCAGGATTATGTTAATGAAAATAGAGAATATGATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGc                                                                                           <  1:4152827/91‑1 (MQ=255)
 tttAAAACTCTTTCAGGATTATGTTAATGAAAATAGAGAATATGATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGc                                                                                           <  1:4151882/91‑1 (MQ=255)
 tttAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGc                                                                                           <  1:3093220/91‑1 (MQ=255)
 tttAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGc                                                                                           <  1:3093017/91‑1 (MQ=255)
  ttAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCt                                                                                          <  1:1894480/91‑1 (MQ=255)
     aaaCTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGaa                                                                                       <  1:1670466/91‑1 (MQ=255)
       aCTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGaaaa                                                                                     <  1:266881/91‑1 (MQ=255)
       aCTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGaaaa                                                                                     <  1:4115430/91‑1 (MQ=255)
         tctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGa                                                                                   <  1:3131601/91‑1 (MQ=255)
         tctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGa                                                                                   >  1:2019794/1‑91 (MQ=255)
          ctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGat                                                                                  >  1:1583481/1‑91 (MQ=255)
          ctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGat                                                                                  >  1:1608124/1‑91 (MQ=255)
          ctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGat                                                                                  >  1:2903947/1‑91 (MQ=255)
          ctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGat                                                                                  <  1:3928115/91‑1 (MQ=255)
           tttAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGatc                                                                                 >  1:2131861/1‑91 (MQ=255)
           tttAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGatc                                                                                 >  1:51713/1‑91 (MQ=255)
               aGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGAtcatca                                                                             >  1:2284914/1‑91 (MQ=255)
               aGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGAtcatca                                                                             >  1:247129/1‑91 (MQ=255)
               aGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGAtcatca                                                                             >  1:2823368/1‑91 (MQ=255)
                 gATTATGTTAATGAAAATAGAGAATATAAGATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGatgatcttt                                                                           <  1:3326014/91‑4 (MQ=38)
                  aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAttt                                                                          >  1:2615235/1‑91 (MQ=255)
                  aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAttt                                                                          >  1:2589040/1‑91 (MQ=255)
                  aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                          <  1:432457/91‑1 (MQ=38)
                  aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                          <  1:2108585/91‑1 (MQ=38)
                  aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                          <  1:45724/91‑1 (MQ=38)
                    tATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTCa                                                                        <  1:1330293/91‑1 (MQ=37)
                      tGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGt                                                                      >  1:3169606/1‑91 (MQ=255)
                      tGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGt                                                                      >  1:3167026/1‑91 (MQ=255)
                       gTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGtt                                                                     >  1:1923203/1‑91 (MQ=255)
                        ttAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTc                                                                    <  1:1757696/91‑1 (MQ=255)
                         tAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCg                                                                   <  1:4023370/91‑1 (MQ=255)
                          aaTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGa                                                                  >  1:1512637/1‑91 (MQ=255)
                          aaTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGa                                                                  <  1:2685490/91‑1 (MQ=255)
                          aaTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGa                                                                  <  1:1485845/91‑1 (MQ=255)
                          aaTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTCAGCTCGa                                                                  <  1:666004/91‑1 (MQ=25)
                           aTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAg                                                                 >  1:2932130/1‑91 (MQ=255)
                            tGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAgg                                                                >  1:3903650/1‑91 (MQ=255)
                               aaaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                             >  1:2719580/1‑91 (MQ=255)
                               aaaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                             >  1:3538921/1‑91 (MQ=255)
                               aaaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                             <  1:1391000/91‑1 (MQ=255)
                                aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttt                                                            <  1:874976/91‑1 (MQ=255)
                                aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttt                                                            <  1:2166680/91‑1 (MQ=255)
                                aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttt                                                            <  1:155164/91‑1 (MQ=255)
                                   agagAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTtttta                                                         >  1:2404449/1‑91 (MQ=255)
                                    gagaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttat                                                        >  1:1496318/1‑91 (MQ=255)
                                    gagaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttat                                                        >  1:4151875/1‑91 (MQ=255)
                                     agaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatt                                                       >  1:467926/1‑91 (MQ=255)
                                      gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAACAGATCATCATTTCAGTTCGAGGTTttttattt                                                      >  1:3532671/1‑91 (MQ=255)
                                      gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                      >  1:2612246/1‑91 (MQ=255)
                                      gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                      >  1:3615799/1‑91 (MQ=255)
                                      gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                      >  1:2526917/1‑91 (MQ=255)
                                      gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                      <  1:2460874/91‑1 (MQ=255)
                                       aatataatatAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatttt                                                     >  1:3005504/1‑91 (MQ=255)
                                        atataatataACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatttta                                                    >  1:2504852/1‑91 (MQ=255)
                                             atataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctc                                               >  1:2985013/1‑91 (MQ=255)
                                             atataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctc                                               <  1:678389/91‑1 (MQ=255)
                                             atataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctc                                               <  1:85276/91‑1 (MQ=255)
                                               ataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCtt                                             <  1:2825468/91‑1 (MQ=255)
                                               ataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCtt                                             <  1:2824575/91‑1 (MQ=255)
                                                taaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCttt                                            <  1:2095825/91‑1 (MQ=255)
                                                  aCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGt                                          <  1:1428408/91‑1 (MQ=255)
                                                  aCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGt                                          <  1:1474307/91‑1 (MQ=255)
                                                     aaTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATg                                       <  1:1506042/91‑1 (MQ=255)
                                                      aTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATgg                                      >  1:624871/1‑91 (MQ=255)
                                                      aTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATgg                                      >  1:96507/1‑91 (MQ=255)
                                                       tCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATggg                                     >  1:1529554/1‑91 (MQ=255)
                                                        cTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGt                                    <  1:4127482/91‑1 (MQ=255)
                                                        cTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGt                                    <  1:4127439/91‑1 (MQ=255)
                                                         tatataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTa                                   <  1:3400376/91‑1 (MQ=255)
                                                          atataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAt                                  >  1:2257374/1‑91 (MQ=255)
                                                           tataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAtt                                 >  1:826304/1‑91 (MQ=255)
                                                           tataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAtt                                 >  1:222210/1‑91 (MQ=255)
                                                            ataaataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTc                                >  1:2275094/1‑91 (MQ=255)
                                                             caaataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCt                               <  1:4097435/90‑1 (MQ=255)
                                                                 taaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTa                           >  1:2979845/1‑91 (MQ=255)
                                                                  aaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTAc                          <  1:1050916/91‑1 (MQ=255)
                                                                   aGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACt                         <  1:1021796/91‑1 (MQ=255)
                                                                   aGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACt                         >  1:2879403/1‑91 (MQ=255)
                                                                       aCCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaa                     >  1:3609822/1‑91 (MQ=255)
                                                                        ccAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaaa                    <  1:1847281/91‑1 (MQ=255)
                                                                            tatTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTg                >  1:3787203/1‑91 (MQ=255)
                                                                            tatTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTg                >  1:3787098/1‑91 (MQ=255)
                                                                               ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGttt             >  1:1985483/1‑91 (MQ=255)
                                                                                tAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtttt            >  1:3715657/1‑91 (MQ=255)
                                                                                tAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtttt            <  1:621451/91‑1 (MQ=255)
                                                                                 aaaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTc           <  1:3986316/91‑1 (MQ=255)
                                                                                 aaaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTc           <  1:2872601/91‑1 (MQ=255)
                                                                                 aaaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTc           <  1:2197789/91‑1 (MQ=255)
                                                                                  aaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTcc          <  1:2398677/91‑1 (MQ=255)
                                                                                  aaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTcc          >  1:2764186/1‑91 (MQ=255)
                                                                                   aCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCt         <  1:1600144/91‑1 (MQ=255)
                                                                                    ccTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTg        <  1:4098634/91‑1 (MQ=255)
                                                                                     cTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGt       <  1:1674233/91‑1 (MQ=255)
                                                                                     cTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGt       <  1:1009095/91‑1 (MQ=255)
                                                                                     cTCGAACTGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGt       >  1:1763664/1‑91 (MQ=25)
                                                                                     cTCGAACTGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGt       >  1:1587355/1‑91 (MQ=25)
                                                                                          gCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGtttttt  >  1:4047049/1‑91 (MQ=255)
                                                                                          gCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGtttttt  >  1:4046648/1‑91 (MQ=255)
                                                                                          gCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGtttaaa  <  1:1228253/91‑4 (MQ=255)
                                                                                          |                                                                                          
ATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTT  >  USA300TCH1516_ALE/168253‑168433

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: