Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 917,327 T→G 11.6% intergenic (+86/+244) sufB_2 → / ← USA300HOU_RS04545 FeS cluster assembly protein SufB/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE917,3270TG11.6% 364.3 / 31.8 112intergenic (+86/+244)sufB_2/USA300HOU_RS04545FeS cluster assembly protein SufB/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (52/47);  new base G (6/7);  total (58/54)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 115 total aligned reads displayed.

ATAGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTAA  >  USA300TCH1516_ALE/917238‑917417
                                                                                           |                                                                                          
  aGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTAc                                                                                         <  1:2661314/91‑1 (MQ=255)
  aGCTTGAAACCGCGTTGTTAAGCCAGTCTTGACTTCCGAGAATGGCTATTGATGCCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTAc                                                                                         <  1:3864131/91‑1 (MQ=255)
   gCTTTAAACGGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACt                                                                                        <  1:3965786/91‑1 (MQ=255)
     tttAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACtgt                                                                                      >  1:2624031/1‑91 (MQ=255)
     tttAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACtgt                                                                                      <  1:2939589/91‑1 (MQ=255)
      ttAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGTCACATTACtgtg                                                                                     <  1:3198827/91‑1 (MQ=255)
      ttAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACtgtg                                                                                     >  1:1518362/1‑91 (MQ=255)
      ttAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACAGTAAtgtg                                                                                     <  1:3140135/91‑1 (MQ=255)
      ttAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACAGTAAtgtg                                                                                     <  1:3892007/91‑1 (MQ=255)
        aaaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGc                                                                                   >  1:293734/1‑91 (MQ=255)
        aaaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGc                                                                                   >  1:1930653/1‑91 (MQ=255)
        aaaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACAGTAATGTGGc                                                                                   <  1:2300848/91‑1 (MQ=255)
        aaaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACAGTAATGTGGc                                                                                   <  1:3055960/91‑1 (MQ=255)
         aaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCt                                                                                  <  1:976136/91‑1 (MQ=255)
         aaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCt                                                                                  <  1:974480/91‑1 (MQ=255)
             gcgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCttttt                                                                              <  1:529192/91‑1 (MQ=255)
                ttgttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAAAGCCACAGTAATGTGGCtttttg                                                                             <  1:3159256/91‑2 (MQ=37)
                ttgttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTg                                                                           >  1:2878837/1‑91 (MQ=255)
                ttgttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTg                                                                           <  1:3711510/91‑1 (MQ=255)
                ttgttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACAGTAATGTGGCtttttgtt                                                                           <  1:4075292/91‑4 (MQ=38)
                 tgttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGt                                                                          <  1:2135786/91‑1 (MQ=255)
                  gttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGtt                                                                         >  1:1781442/1‑91 (MQ=255)
                  gttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGtt                                                                         >  1:2049917/1‑91 (MQ=255)
                   ttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATAAAACAAAAAAAGCCACAGTAATGTGGCtttttgtta                                                                          <  1:877259/91‑5 (MQ=18)
                    tAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGtttt                                                                       <  1:1123408/91‑1 (MQ=255)
                     aaGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTta                                                                      >  1:2699162/1‑91 (MQ=255)
                     aaGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTta                                                                      >  1:2517928/1‑91 (MQ=255)
                       gCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTtata                                                                    <  1:2838527/91‑1 (MQ=255)
                       gCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTtata                                                                    <  1:1697510/91‑1 (MQ=255)
                         cATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAAc                                                                  <  1:1427289/91‑1 (MQ=255)
                          aTTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACt                                                                 <  1:1132821/91‑1 (MQ=255)
                           ttCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTa                                                                >  1:1245799/1‑91 (MQ=255)
                           ttCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTa                                                                >  1:1669528/1‑91 (MQ=255)
                           ttCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTa                                                                >  1:387372/1‑91 (MQ=255)
                           ttCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTa                                                                >  1:1245001/1‑91 (MQ=255)
                             cTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTaaa                                                              <  1:951708/91‑1 (MQ=255)
                             cTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTaaa                                                              >  1:3050782/1‑91 (MQ=255)
                              ttGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAAt                                                             >  1:957402/1‑91 (MQ=255)
                                gACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCg                                                           >  1:1154070/1‑91 (MQ=255)
                                gACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCg                                                           <  1:2656597/91‑1 (MQ=255)
                                gACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCg                                                           <  1:2656598/91‑1 (MQ=255)
                                gACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCg                                                           <  1:1943967/91‑1 (MQ=255)
                                 aCTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCgg                                                          >  1:2790454/1‑91 (MQ=255)
                                  cTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGGTTTATTACCAAAACggg                                                         >  1:4053079/1‑90 (MQ=255)
                                       gAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgat                                                    >  1:1035495/1‑91 (MQ=255)
                                       gAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgat                                                    <  1:1700123/91‑1 (MQ=255)
                                        aaaaTGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgata                                                   <  1:2933791/91‑1 (MQ=255)
                                            tGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgatagata                                               <  1:3956750/91‑1 (MQ=255)
                                            tGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgatagata                                               <  1:3957300/91‑1 (MQ=255)
                                               cTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGc                                            <  1:3099234/91‑1 (MQ=255)
                                                   tGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTg                                        <  1:344842/91‑1 (MQ=255)
                                                     aTACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATGACTGTGGCCTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTa                                      >  1:3303949/1‑91 (MQ=255)
                                                       aCCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACt                                    >  1:513383/1‑91 (MQ=255)
                                                       aCCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACt                                    >  1:513462/1‑91 (MQ=255)
                                                        ccATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACtt                                   >  1:540996/1‑91 (MQ=255)
                                                        ccATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACtt                                   >  1:2395092/1‑91 (MQ=255)
                                                        ccATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACtt                                   >  1:3470463/1‑91 (MQ=255)
                                                         cATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACtta                                  >  1:2349707/1‑91 (MQ=255)
                                                         cATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACtta                                  >  1:3784731/1‑91 (MQ=255)
                                                            tttGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttattt                               <  1:2694304/91‑1 (MQ=255)
                                                             ttGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttattta                              >  1:1314034/1‑91 (MQ=255)
                                                                    tGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                       <  1:3053809/91‑1 (MQ=255)
                                                                    tGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                       >  1:1113696/1‑91 (MQ=255)
                                                                    tGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                       <  1:42727/91‑1 (MQ=255)
                                                                    tGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                       <  1:2104105/91‑1 (MQ=255)
                                                                     gAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGtt                      >  1:1743600/1‑91 (MQ=255)
                                                                      aaTATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTc                     <  1:3794421/91‑1 (MQ=255)
                                                                         ataAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAt                  >  1:4026730/1‑91 (MQ=255)
                                                                         ataAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAt                  <  1:2219755/91‑1 (MQ=255)
                                                                          taAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtt                 >  1:3468982/1‑91 (MQ=255)
                                                                          taAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtt                 >  1:342834/1‑91 (MQ=255)
                                                                           aaC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAttt                >  1:356895/1‑91 (MQ=255)
                                                                            aC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtttt               <  1:2665246/91‑1 (MQ=255)
                                                                            aC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtttt               <  1:828018/91‑1 (MQ=255)
                                                                            aC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtttt               <  1:1903357/91‑1 (MQ=255)
                                                                            aC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtttt               <  1:2466070/91‑1 (MQ=255)
                                                                            aC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtttt               <  1:2465940/91‑1 (MQ=255)
                                                                            aC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTAGTACTTATTTATATCAGTTCGCtttt               >  1:637762/1‑91 (MQ=255)
                                                                             c‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAttttt              >  1:1177153/1‑91 (MQ=255)
                                                                             c‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAttttt              >  1:1177186/1‑91 (MQ=255)
                                                                             c‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAttttt              >  1:691614/1‑91 (MQ=255)
                                                                                aaaaaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtttttt             >  1:513211/1‑91 (MQ=255)
                                                                                 aaaaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTg            >  1:1267283/1‑91 (MQ=255)
                                                                                  aaaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGa           <  1:1688587/91‑1 (MQ=255)
                                                                                   aaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGAt          >  1:2274193/1‑91 (MQ=255)
                                                                                   aaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGAt          <  1:194659/91‑1 (MQ=255)
                                                                                   aaGCCACAGTAATGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGAt          >  1:2044286/1‑91 (MQ=255)
                                                                                      ccACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTgg       <  1:534469/91‑1 (MQ=255)
                                                                                      ccACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTgg       >  1:2293358/1‑91 (MQ=255)
                                                                                      ccACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTgg       >  1:1818383/1‑91 (MQ=255)
                                                                                      ccACATTACTGTGGCTTTTTTTGTTTTATAACTAAAACGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTgg       >  1:1644254/1‑91 (MQ=255)
                                                                                      ccACAGTAATGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTgg       >  1:3000398/1‑91 (MQ=255)
                                                                                       cacaGTAATGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgt      >  1:937441/1‑91 (MQ=255)
                                                                                        acaTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtg     <  1:150084/91‑1 (MQ=255)
                                                                                        acaTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtg     >  1:887488/1‑91 (MQ=255)
                                                                                         caTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtgt    <  1:2011908/91‑1 (MQ=255)
                                                                                         caTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtgt    <  1:967804/91‑1 (MQ=255)
                                                                                         caTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtgt    >  1:3715589/1‑91 (MQ=255)
                                                                                           ttACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTaa  >  1:903745/1‑91 (MQ=255)
                                                                                           ttACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTaa  >  1:903865/1‑91 (MQ=255)
                                                                                           |                                                                                          
ATAGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTAA  >  USA300TCH1516_ALE/917238‑917417

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: