Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 917,330 C→A 10.8% intergenic (+89/+241) 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,3300CA10.8% 387.1 / 26.3 116intergenic (+89/+241)sufB_2/USA300HOU_RS04545FeS cluster assembly protein SufB/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (54/49);  new base A (6/7);  total (60/56)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

AGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTAAAAA  >  USA300TCH1516_ALE/917240‑917420
                                                                                            |                                                                                          
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:2939589/91‑1 (MQ=255)
   tttAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACtgt                                                                                         >  1:2624031/1‑91 (MQ=255)
    ttAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGTCACATTACtgtg                                                                                        <  1:3198827/91‑1 (MQ=255)
    ttAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACtgtg                                                                                        >  1:1518362/1‑91 (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:3055960/91‑1 (MQ=255)
      aaaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACAGTAATGTGGc                                                                                      <  1:2300848/91‑1 (MQ=255)
       aaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACtgtgacc                                                                                     >  1:1528332/1‑88 (MQ=255)
       aaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCt                                                                                     <  1:976136/91‑1 (MQ=255)
       aaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCt                                                                                     <  1:974480/91‑1 (MQ=255)
       aaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCt                                                                                     <  1:93236/91‑1 (MQ=255)
           gcgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCttttt                                                                                 <  1:529192/91‑1 (MQ=255)
              ttgttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAAAGCCACAGTAATGTGGCtttttg                                                                                <  1:3159256/91‑2 (MQ=37)
              ttgttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTg                                                                              <  1:3711510/91‑1 (MQ=255)
              ttgttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTg                                                                              >  1:2019464/1‑91 (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:1697510/91‑1 (MQ=255)
                     gCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTtata                                                                       <  1:2838527/91‑1 (MQ=255)
                      ccATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATaa                                                                      >  1:3696564/1‑91 (MQ=255)
                       cATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAAc                                                                     <  1:1213016/91‑1 (MQ=255)
                        aTTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACt                                                                    <  1:1132821/91‑1 (MQ=255)
                         ttCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTa                                                                   >  1:1245799/1‑91 (MQ=255)
                         ttCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTa                                                                   >  1:387372/1‑91 (MQ=255)
                         ttCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTa                                                                   >  1:1245001/1‑91 (MQ=255)
                           cTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTaaa                                                                 >  1:3050782/1‑91 (MQ=255)
                           cTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTaaa                                                                 <  1:951708/91‑1 (MQ=255)
                            ttGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAAt                                                                >  1:957402/1‑91 (MQ=255)
                              gACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCg                                                              <  1:1943967/91‑1 (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)
                               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:3957300/91‑1 (MQ=255)
                                             cTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGc                                               <  1:3099234/91‑1 (MQ=255)
                                                 tGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTg                                           <  1:344842/91‑1 (MQ=255)
                                                   aTACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTa                                         <  1:2755940/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:3470463/1‑91 (MQ=255)
                                                      ccATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACtt                                      >  1:2395092/1‑91 (MQ=255)
                                                       cATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACtta                                     >  1:2349707/1‑91 (MQ=255)
                                                          tttGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttattt                                  <  1:2694304/91‑1 (MQ=255)
                                                          tttGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCGTTGTACttattt                                  <  1:3174804/91‑1 (MQ=255)
                                                           ttGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttattta                                 >  1:1314034/1‑91 (MQ=255)
                                                                  tGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                          <  1:3053809/91‑1 (MQ=255)
                                                                  tGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                          <  1:42727/91‑1 (MQ=255)
                                                                  tGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                          >  1:1113696/1‑91 (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:342834/1‑91 (MQ=255)
                                                                         aaC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAttt                   >  1:356895/1‑91 (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:691614/1‑91 (MQ=255)
                                                                            aaaaaaaGCCACAGTAATGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtttt                  >  1:2573210/3‑91 (MQ=38)
                                                                              aaaaaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtttttt                >  1:513211/1‑91 (MQ=255)
                                                                               aaaaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTg               >  1:1267283/1‑91 (MQ=255)
                                                                               aaaaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTg               >  1:564378/1‑91 (MQ=255)
                                                                                aaaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGa              <  1:390465/91‑1 (MQ=255)
                                                                                aaaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGa              <  1:1688587/91‑1 (MQ=255)
                                                                                 aaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGAt             >  1:2274193/1‑91 (MQ=255)
                                                                                 aaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGAt             <  1:194659/91‑1 (MQ=255)
                                                                                  aGCCACAGTAATGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGAtt            >  1:302325/1‑91 (MQ=255)
                                                                                    ccACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTgg          >  1:2294657/1‑91 (MQ=255)
                                                                                    ccACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTgg          <  1:534469/91‑1 (MQ=255)
                                                                                    ccACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTgg          >  1:1818383/1‑91 (MQ=255)
                                                                                    ccACATTACTGTGGCTTTTTTTGTTTTATAACTAAAACGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTgg          >  1:1644254/1‑91 (MQ=255)
                                                                                    ccACAGTAATGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTgg          >  1:3000398/1‑91 (MQ=255)
                                                                                      acaTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtg        <  1:150084/91‑1 (MQ=255)
                                                                                      acaTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtg        >  1:887488/1‑91 (MQ=255)
                                                                                       caTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtgt       <  1:967804/91‑1 (MQ=255)
                                                                                       caTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtgt       <  1:2011908/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)
                                                                                           aCTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTaaaa   <  1:3009771/91‑1 (MQ=255)
                                                                                           aCTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTaaaa   <  1:764165/91‑1 (MQ=255)
                                                                                           aCTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTaaaa   <  1:702739/91‑1 (MQ=255)
                                                                                            cTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTGCGATTTTTTGATTGGTGTaaaaa  <  1:1663598/91‑1 (MQ=255)
                                                                                            |                                                                                          
AGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTAAAAA  >  USA300TCH1516_ALE/917240‑917420

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: