Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 419,965 A→T 7.1% intergenic (+183/+83) rpsR → / ← USA300HOU_RS01950 30S ribosomal protein S18/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE419,9650AT 7.1% 695.4 / 29.2 211intergenic (+183/+83)rpsR/USA300HOU_RS0195030S ribosomal protein S18/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (110/85);  new base T (14/2);  total (124/87)
Fisher's exact test for biased strand distribution p-value = 1.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 233 total aligned reads displayed.

GTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGTTTATACTT  >  USA300TCH1516_ALE/419876‑420055
                                                                                            |                                                                                          
gTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTtaaat                                                                                           <  1:3749567/91‑1 (MQ=255)
   ttGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTtaaataaa                                                                                        >  1:2904354/1‑91 (MQ=255)
     gtgtTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAGAAt                                                                                      >  1:3061164/1‑91 (MQ=255)
     gtgtTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAAt                                                                                      >  1:790825/1‑91 (MQ=255)
      tgtTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATg                                                                                     >  1:181044/1‑91 (MQ=255)
      tgtTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATg                                                                                     >  1:2405310/1‑91 (MQ=255)
       gtTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGc                                                                                    <  1:332599/91‑1 (MQ=255)
        ttCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCa                                                                                   <  1:1796057/91‑1 (MQ=255)
          cGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATa                                                                                 >  1:292013/1‑91 (MQ=255)
             aaTTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGct                                                                              <  2:3278753/91‑1 (MQ=255)
             aaTTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGct                                                                              >  2:1027476/1‑91 (MQ=255)
              aTTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGctc                                                                             >  2:3725215/1‑91 (MQ=255)
              aTTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGctc                                                                             <  1:2050880/91‑1 (MQ=255)
               tttGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGctct                                                                            <  1:3633276/91‑1 (MQ=255)
               tttGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGctct                                                                            <  2:2965783/91‑1 (MQ=255)
                ttGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCtctt                                                                           <  2:3855291/91‑1 (MQ=255)
                 tGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAT‑‑‑CTACGCATTTTACATAAAATGCATAGCtcttc                                                                          <  2:165628/91‑1 (MQ=255)
                 tGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCtcttc                                                                          <  1:2181529/91‑1 (MQ=255)
                 tGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCtcttc                                                                          >  2:816457/1‑91 (MQ=255)
                 tGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCtcttc                                                                          <  2:1756169/91‑1 (MQ=255)
                  gCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCtcttct                                                                         >  1:3816375/1‑91 (MQ=255)
                     caAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCtttt                                                                      >  2:827392/1‑91 (MQ=255)
                     aaaaGATATGTATATATAGAAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCtttt                                                                      <  2:3633960/90‑1 (MQ=255)
                       aaGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAG‑‑‑CTACGCAGTTTAAATAAAATGCATAGCTCTTCTTTTTc                                                                    <  2:1426013/91‑1 (MQ=255)
                       aaGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAAGAGCTATGCATTTTATTTAAAATGCGTAGTTCtttttt                                                                       <  1:2486327/91‑1 (MQ=9)
                       aaGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTc                                                                    >  1:3850476/1‑91 (MQ=255)
                        aGATATGTTTATATTGCAAAAATAATATGAATTGAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCt                                                                   <  2:2464044/91‑1 (MQ=255)
                           tatGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGc                                                                >  1:872520/1‑91 (MQ=255)
                             tGTTTATATTGCAAACCTAATATGAATTTAGATGCATCAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCATTTTCTTGCac                                                              >  2:2579798/1‑90 (MQ=255)
                             tGTTTATATTGCAAAAATAATATGAATTTAGATGCATCAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCAt                                                              >  1:3344687/1‑91 (MQ=255)
                               tttATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATAc                                                            >  2:1202217/1‑91 (MQ=255)
                               tttATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATAc                                                            <  1:3128305/91‑1 (MQ=255)
                               tttATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATAc                                                            >  1:414949/1‑91 (MQ=255)
                                ttATATTGCAAAAATAATATGGATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATGTTAATTAAAATGCATAGCTCTTCTTTTTCTTGCATACg                                                           <  1:3304664/91‑1 (MQ=255)
                                ttATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACg                                                           >  1:2459104/1‑91 (MQ=255)
                                 tataTTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCACTTTAAATAAAATGCCTAGCTCTTCTTTTTCCTGCATACGa                                                          >  2:772883/1‑91 (MQ=255)
                                      tGCAAAAATAATATGAATTTAGATGCATAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTaaa                                                    <  1:528888/91‑1 (MQ=255)
                                        cAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTGCGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTaaaa                                                   <  2:1118904/91‑1 (MQ=255)
                                        cAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTaaaa                                                   >  1:1996871/1‑91 (MQ=255)
                                         aaaaaTAATATGAATTTAGATGCATAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATa                                                 >  1:2898571/1‑91 (MQ=255)
                                              taataTGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTc                                             >  2:2364687/1‑91 (MQ=255)
                                              taataTGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTc                                             >  1:265782/1‑91 (MQ=255)
                                                atatGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAAATcgc                                           >  2:2403891/1‑91 (MQ=255)
                                                 tatGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTcgcg                                          >  2:2224091/1‑91 (MQ=255)
                                                 tatGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTcgcg                                          <  2:1471278/91‑1 (MQ=255)
                                                 tatGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTcgcg                                          >  1:142364/1‑91 (MQ=255)
                                                 tatGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTcgcg                                          >  1:826849/1‑91 (MQ=255)
                                                  atGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCga                                         >  1:2239243/1‑91 (MQ=255)
                                                  atGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCga                                         >  2:3357128/1‑91 (MQ=255)
                                                  atGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCga                                         >  1:497085/1‑91 (MQ=255)
                                                    gAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCgaga                                       <  2:1943623/91‑1 (MQ=255)
                                                    gAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCgaga                                       <  2:678018/91‑1 (MQ=255)
                                                    gAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCgaga                                       >  1:2752946/1‑91 (MQ=255)
                                                      aTTTAGATGCATAACAAAAGAA‑‑‑CGTCGCATTTTAAATAAAATGGATAGCTGTTCCTTTTCTTGCATACGAATTAAAATAACTCTCGAGAcc                                     <  2:3680914/91‑1 (MQ=255)
                                                         tAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCtat                                  <  2:3235963/91‑1 (MQ=255)
                                                         tAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCtat                                  <  1:3176213/91‑1 (MQ=255)
                                                          aGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCtata                                 <  2:3103706/91‑1 (MQ=255)
                                                           ggggCAGAAAAGAAGAA‑‑‑CTTCGCGTTGTAAATAAAATGGATAGCTCTTCTGTTTCTGGCATACGAATGAAAATCACTCGCGAGACCTATaa                                <  2:120598/88‑1 (MQ=255)
                                                           gATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATaa                                >  1:1228836/1‑91 (MQ=255)
                                                           gATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATaa                                <  2:1811272/91‑1 (MQ=255)
                                                           gATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATaa                                >  1:3116421/1‑91 (MQ=255)
                                                           gATGCATAAAAAAAGAA‑‑‑CTACGCATGTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATaa                                <  2:3602708/91‑1 (MQ=255)
                                                            aTGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAg                               <  2:567830/91‑1 (MQ=255)
                                                            aTGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAg                               >  2:2106655/1‑91 (MQ=255)
                                                               cATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAAGTAAAATAACTCGCGAGACCTATAAGtct                            <  2:3171558/91‑1 (MQ=255)
                                                                aTAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGtctc                           <  2:1981586/91‑1 (MQ=255)
                                                                  aaaaaaaGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCtt                         <  2:3402191/91‑1 (MQ=255)
                                                                   aaaaaaGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCttt                        <  2:3384504/91‑1 (MQ=255)
                                                                     aaaaGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTcc                      >  1:1886987/1‑91 (MQ=255)
                                                                     aaaaGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTcc                      >  1:2322127/1‑91 (MQ=255)
                                                                      gaaGAG‑‑‑CTATGCATTTTATTTAAAATGCGTAGTTCTTTTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCt                     >  2:782807/2‑91 (MQ=18)
                                                                      gaaGAG‑‑‑CTATGCATTTTATTTAAAATGCGTAATTCTTTTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCt                     >  1:3506946/2‑91 (MQ=18)
                                                                      aaaGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCt                     <  2:3459642/91‑1 (MQ=255)
                                                                       aaGAG‑‑‑CTATGCATTTTATTTAAAATGCGTAGTTCTTTTTTTTCTTGCCTACGAATTAAAATAACTCGCGCGCCCTATACGTCTCTTTCCTc                    >  2:426631/1‑91 (MQ=12)
                                                                       aaGAG‑‑‑CTATGCATTTTATTTAAAATGCGTAGTTCTTTTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTc                    >  1:1706919/1‑91 (MQ=18)
                                                                       aaGAG‑‑‑CTATGCATTTTATTTAAAATGCGTAGTTCTTTTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTc                    >  1:2974805/1‑91 (MQ=18)
                                                                         gAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCAc                  >  2:1319289/1‑91 (MQ=255)
                                                                         gAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCAc                  >  2:2847610/1‑91 (MQ=255)
                                                                          aa‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACt                 >  1:1115756/1‑91 (MQ=255)
                                                                           a‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCTTACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTa                >  2:2625755/1‑91 (MQ=255)
                                                                           a‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATACGTCTCTTTCATCACTa                >  2:2893014/1‑91 (MQ=255)
                                                                           a‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTa                <  1:1853765/91‑1 (MQ=255)
                                                                              gcTATGCATTTTATTTAAAATGCGTAGTTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTa                >  1:666227/2‑91 (MQ=25)
                                                                               cTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAg               >  2:2179828/1‑91 (MQ=255)
                                                                               cTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAg               <  1:1587486/91‑1 (MQ=255)
                                                                                 aCGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGAt             <  2:2705931/91‑1 (MQ=255)
                                                                                    cATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGt          <  2:1498391/91‑1 (MQ=255)
                                                                                    cATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGt          <  1:935120/91‑1 (MQ=255)
                                                                                     aTTTTATTTAAAATGCGTAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGtt         >  1:1331601/1‑91 (MQ=255)
                                                                                     aTTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGtt         <  2:1583339/91‑1 (MQ=255)
                                                                                        ttAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGTTtat      <  2:1326318/91‑1 (MQ=255)
                                                                                        ttAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGTTtat      <  1:2023971/91‑1 (MQ=255)
                                                                                         taaataaaATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGTTtata     >  1:1924548/1‑91 (MQ=255)
                                                                                          aaataaaATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGTTTATAc    >  1:1050478/1‑91 (MQ=255)
                                                                                          aaataaaATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGTTTATAc    >  1:1667347/1‑91 (MQ=255)
                                                                                           aataaaATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGTTTATACt   <  1:2931679/91‑1 (MQ=255)
                                                                                           aataaaATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGTTTATACt   <  1:3558688/91‑1 (MQ=255)
                                                                                           aataaaATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGTTTATACt   <  1:1568565/91‑1 (MQ=255)
                                                                                           aataaaATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGTTTATACt   >  2:1379489/1‑91 (MQ=255)
                                                                                            ataaaATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGTTTATACtt  >  2:3054908/1‑91 (MQ=255)
                                                                                            |                                                                                          
GTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAA‑‑‑CTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATAGTTTATACTT  >  USA300TCH1516_ALE/419876‑420055

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