Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 917,327 T→G 7.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,3270TG 7.6% 734.7 / 37.2 227intergenic (+86/+244)sufB_2/USA300HOU_RS04545FeS cluster assembly protein SufB/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (113/96);  new base G (5/13);  total (118/109)
Fisher's exact test for biased strand distribution p-value = 4.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 235 total aligned reads displayed.

ATAGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTAA  >  USA300TCH1516_ALE/917238‑917417
                                                                                           |                                                                                          
aTAGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACAtt                                                                                           >  1:1813480/1‑91 (MQ=255)
   gCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACt                                                                                        <  2:1754664/91‑1 (MQ=255)
    cTTTAAGCCGCGTTGTTAAGCCATTATTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAGC‑‑AAACAGCCACAGTAAtg                                                                                       <  2:942181/91‑1 (MQ=255)
    cTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACtg                                                                                       >  2:1072647/1‑91 (MQ=255)
    cTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACtg                                                                                       >  1:1493953/1‑91 (MQ=255)
    cTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACtg                                                                                       >  2:2484610/1‑91 (MQ=255)
      ttAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACtgtg                                                                                     >  1:3383141/1‑91 (MQ=255)
      ttAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTAATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACAt                                                                                            >  1:3072817/1‑84 (MQ=255)
       tAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACAGTAATGTgg                                                                                    <  1:1338316/91‑1 (MQ=255)
       tAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACAGTAATGTgg                                                                                    <  1:204270/91‑1 (MQ=255)
        aaaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGc                                                                                   >  2:1482913/1‑91 (MQ=255)
        aaaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGc                                                                                   >  1:1526021/1‑91 (MQ=255)
        aaaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACAGTAATGTGGc                                                                                   <  1:2513532/91‑1 (MQ=255)
         aaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCt                                                                                  <  1:298846/91‑1 (MQ=255)
         aaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCt                                                                                  <  1:3440632/91‑1 (MQ=255)
         aaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCt                                                                                  >  1:1991727/1‑91 (MQ=255)
         aaCCGCGTGGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCt                                                                                  <  2:3593369/91‑1 (MQ=255)
           ccGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCttt                                                                                >  1:2605748/1‑91 (MQ=255)
            cgcgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACAGTAATGTGGCtttt                                                                               <  2:1330511/91‑1 (MQ=255)
             gcgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCttttt                                                                              <  1:569872/91‑1 (MQ=255)
             gcgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCttttt                                                                              <  1:3815239/91‑1 (MQ=255)
             gcgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCttttt                                                                              >  1:1355664/1‑91 (MQ=255)
             gcgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCttttt                                                                              <  1:1860899/91‑1 (MQ=255)
             gcgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCttttt                                                                              >  2:2825712/1‑91 (MQ=255)
             gcgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCttttt                                                                              >  2:2917063/1‑91 (MQ=255)
                  gttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAAAGCCACAGTAATGTGGCtttttgtt                                                                           <  2:1520270/91‑4 (MQ=25)
                  gttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGtt                                                                         >  2:2725187/1‑91 (MQ=255)
                  gttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGtt                                                                         >  2:1126008/1‑91 (MQ=255)
                  gatAAGCCATTCTTGACTGCCGAAAATGGCTATTGATACCATTTTGCTACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGtt                                                                         >  2:1385073/3‑91 (MQ=255)
                   ttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCATATTACTGTGGCTTTTTTTGttt                                                                        >  2:3120305/1‑91 (MQ=255)
                   ttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGttt                                                                        >  1:741574/1‑91 (MQ=255)
                   ttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGttt                                                                        >  2:2121106/1‑91 (MQ=255)
                   ttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGATTTTTTTGttt                                                                        >  1:3515282/1‑91 (MQ=255)
                     aaGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTta                                                                      >  1:1134059/1‑91 (MQ=255)
                         cATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAAc                                                                  <  1:2423510/91‑1 (MQ=255)
                            tCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTaa                                                               >  2:3679735/1‑91 (MQ=255)
                             cTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTaaa                                                              <  1:3056915/91‑1 (MQ=255)
                              ttGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAAt                                                             >  1:3031916/1‑91 (MQ=255)
                              ttGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAAt                                                             <  1:3754082/91‑1 (MQ=255)
                              ttGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAAt                                                             >  2:1024747/1‑91 (MQ=255)
                                     ccGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTg                                                      <  2:1439562/91‑1 (MQ=255)
                                      cGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTga                                                     >  1:2324081/1‑91 (MQ=255)
                                       gAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgat                                                    <  1:918592/91‑1 (MQ=255)
                                         aaaTGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATGACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgatag                                                  <  1:3363124/91‑1 (MQ=255)
                                          aaTGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgataga                                                 <  1:3252609/91‑1 (MQ=255)
                                          aaTGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgataga                                                 <  2:621788/91‑1 (MQ=255)
                                             ggCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATaa                                              >  2:3642538/1‑91 (MQ=255)
                                             ggCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATaa                                              <  2:3263591/91‑1 (MQ=255)
                                              gCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAg                                             <  1:120545/91‑1 (MQ=255)
                                               cTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGc                                            <  1:1473062/91‑1 (MQ=255)
                                                tATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCt                                           >  2:1869504/1‑91 (MQ=255)
                                                 aTTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCtt                                          >  2:1909218/1‑91 (MQ=255)
                                                 aTTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCtt                                          <  1:1011823/91‑1 (MQ=255)
                                                 aTTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCtt                                          >  1:3832776/1‑91 (MQ=255)
                                                 aTTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCtt                                          <  2:1665057/91‑1 (MQ=255)
                                                 aTTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCtt                                          >  2:3134501/1‑91 (MQ=255)
                                                 aTTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCtt                                          >  2:2608011/1‑91 (MQ=255)
                                                  ttGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCttt                                         >  2:18637/1‑91 (MQ=255)
                                                   tGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTg                                        <  1:3783047/91‑1 (MQ=255)
                                                    gATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGt                                       >  2:2365291/1‑91 (MQ=255)
                                                     aTACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTa                                      <  2:2997839/91‑1 (MQ=255)
                                                      tACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTAc                                     >  1:2283391/1‑91 (MQ=255)
                                                       aCCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACt                                    <  1:74845/91‑1 (MQ=255)
                                                       aCCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACt                                    <  1:824385/91‑1 (MQ=255)
                                                       aCCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACt                                    <  1:1964012/91‑1 (MQ=255)
                                                       aCCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACt                                    <  1:2079731/91‑1 (MQ=255)
                                                         cATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACtta                                  >  1:517442/1‑91 (MQ=255)
                                                          aTTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttat                                 >  2:3227208/1‑91 (MQ=255)
                                                          aTTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttat                                 >  2:2529694/1‑91 (MQ=255)
                                                            tttGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttattt                               >  1:2820427/1‑91 (MQ=255)
                                                             ttGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttattta                              >  2:2281377/1‑91 (MQ=255)
                                                              tGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttatttat                             >  2:1809251/1‑91 (MQ=255)
                                                               gATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTtata                            >  2:685711/1‑91 (MQ=255)
                                                               gATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTtata                            >  2:546454/1‑91 (MQ=255)
                                                                 tACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATc                          <  2:3336506/91‑1 (MQ=255)
                                                                    tGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                       >  1:2553580/1‑91 (MQ=255)
                                                                    tGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                       <  1:2671788/91‑1 (MQ=255)
                                                                    tGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                       <  2:885401/91‑1 (MQ=255)
                                                                        tataAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGa                   <  1:2140592/91‑1 (MQ=255)
                                                                          taAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtt                 <  1:2261134/91‑1 (MQ=255)
                                                                          taAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtt                 <  1:858036/91‑1 (MQ=255)
                                                                          taAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtt                 <  2:2769601/91‑1 (MQ=255)
                                                                             c‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAttttt              <  2:2312594/91‑1 (MQ=255)
                                                                             c‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAttttt              <  2:2051420/91‑1 (MQ=255)
                                                                                aaaaaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAg                        <  2:2029292/80‑1 (MQ=255)
                                                                                aaaaaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGActtttt             >  1:3616843/1‑91 (MQ=255)
                                                                                 aaaaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTg            >  1:1300056/1‑91 (MQ=255)
                                                                                   aaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGAt          >  1:1689565/1‑91 (MQ=255)
                                                                                    aGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGAtt         <  1:2127827/91‑1 (MQ=255)
                                                                                     gCCACAGTAATGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATGTATATAAGTTCGATTTTTTGATTg        >  2:1063723/1‑91 (MQ=255)
                                                                                       cacaTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATCGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTTATTGgt      >  2:480262/1‑91 (MQ=255)
                                                                                       cacaTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgt      >  1:387299/1‑91 (MQ=255)
                                                                                       cacaTTACTGTCGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACCTATTTATATCAGCTCTCCTTTttggttggc      >  2:480431/1‑90 (MQ=255)
                                                                                       cacaTAACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTATGAACTTATTTATATCCGTTCGATTTTTTGATCGgt      >  2:3507095/1‑91 (MQ=255)
                                                                                        acaTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtg     >  1:41025/1‑91 (MQ=255)
                                                                                        acaTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtg     <  1:2061111/91‑1 (MQ=255)
                                                                                        acaTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtg     >  1:2599875/1‑91 (MQ=255)
                                                                                         caTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtgt    >  2:2943110/1‑91 (MQ=255)
                                                                                          aTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTa   <  1:2808957/91‑1 (MQ=255)
                                                                                          aTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTa   <  1:1072647/91‑1 (MQ=255)
                                                                                           |                                                                                          
ATAGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAAC‑‑AAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTAA  >  USA300TCH1516_ALE/917238‑917417

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