Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,899,097 C→G 6.3% intergenic (‑88/+405) ribD ← / ← USA300HOU_RS09405 Riboflavin biosynthesis protein RibD/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,899,0970CG 6.3% 501.7 / 19.1 143intergenic (‑88/+405)ribD/USA300HOU_RS09405Riboflavin biosynthesis protein RibD/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (67/67);  new base G (4/5);  total (71/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

Only 100 of 143 total aligned reads displayed.

TAAATTGACTCAAATGATCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTAACTGTCGGCTCTA  >  USA300TCH1516_ALE/1899009‑1899186
                                                                                        |                                                                                         
tAAATTGACTCAAATGATCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTa                                                                                         >  1:2283221/1‑91 (MQ=255)
tAAATTGACTCAAATGATCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTa                                                                                         <  1:3001585/91‑1 (MQ=255)
tAAATTGACTCAAATGATCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTa                                                                                         <  1:1894047/91‑1 (MQ=255)
     tGACTCAAATGATCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGt                                                                                    <  1:203115/91‑1 (MQ=255)
     tGACTCAAATGATCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGt                                                                                    <  1:3465542/91‑1 (MQ=255)
     tGACTCAAATGATCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGt                                                                                    <  1:876968/91‑1 (MQ=255)
      gACTCAAACGATCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTa                                                                                   >  1:1543935/1‑91 (MQ=255)
       aCTCAAATGATCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTac                                                                                  <  1:208032/91‑1 (MQ=255)
        cTCAAATGATCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTaca                                                                                 <  1:3319153/91‑1 (MQ=255)
         tCAAATGATCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTacac                                                                                <  1:3145705/91‑1 (MQ=255)
           aaaTGATCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACAt                                                                              <  1:1852611/91‑1 (MQ=255)
              tGATCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACAtttt                                                                           <  1:755339/91‑1 (MQ=255)
              tGATCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACAtttt                                                                           >  1:2711085/1‑91 (MQ=255)
               gtgagattctagagccgacagttaaagtctggatgggagaaagaatgttaattatcgacaaagataatgtAGCGTATTTGTACACAttttt                                                                          >  1:631193‑M1/71‑91 (MQ=255)
               gtgagattctagagccgacagttaaagtctggatgggagaaagaatgttaattatcgacaaagataatgtAGCGTATTTGTACACAttttt                                                                          >  1:667816‑M1/71‑91 (MQ=255)
               gtgagatgctagagccgacagttaaagtctggatgggagaaagaatgttaattatcgacaaagataatgtAGCGTATTTGTACACAttttt                                                                          >  1:2086705‑M1/71‑91 (MQ=255)
               gATCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACAttttt                                                                          <  1:1088154/91‑1 (MQ=255)
                 tCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTAc                                                                        >  1:3911770/1‑91 (MQ=255)
                 tCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTAc                                                                        >  1:1635945/1‑91 (MQ=255)
                  cACCTCTATAAACAAAAATCCTCGCGTCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACa                                                                       <  1:1837383/91‑1 (MQ=255)
                  cACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAATGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACa                                                                       >  1:2440537/1‑91 (MQ=255)
                  cACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACa                                                                       >  1:2434133/1‑91 (MQ=255)
                  cACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACa                                                                       >  1:560944/1‑91 (MQ=255)
                  cACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACa                                                                       >  1:3664989/1‑91 (MQ=255)
                  cACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACa                                                                       >  1:3884077/1‑91 (MQ=255)
                  cACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACa                                                                       >  1:2433398/1‑91 (MQ=255)
                  cACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACa                                                                       >  1:1444082/1‑91 (MQ=255)
                    ccTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACaaa                                                                     >  1:291606/1‑91 (MQ=255)
                    ccTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACaaa                                                                     >  1:2805817/1‑91 (MQ=255)
                    ccTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACaaa                                                                     >  1:752048/1‑91 (MQ=255)
                     ttctagagccgacagttaaagtctggatgggagaaagaatgttaattatcgacaaagataatgtAGCGTATTTGTACACATTTTTACAAAt                                                                    >  1:3810315‑M1/65‑91 (MQ=255)
                     ttctagagccgacagttaaagtctggatgggagaaagaatgttaattatcgacaaagataatgtAGCCTATTTGTACACATTTTTACAAAt                                                                    >  1:3164098‑M1/65‑91 (MQ=255)
                         ataAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTAAAAATGTGTACAAATaggc                                                                <  1:2771713/91‑4 (MQ=37)
                         ataAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTAAAAATGTGTACAAATaggc                                                                <  1:2771700/91‑4 (MQ=37)
                           aaaCAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTAAAAATGTGTACAAATAGGCtt                                                              <  1:1124614/91‑2 (MQ=37)
                           aaaCAAAAATCCGCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTAAAAATGTGTACAAATAGGCtt                                                              <  1:1123815/91‑2 (MQ=25)
                            aaCAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTAc                                                             >  1:2047488/1‑91 (MQ=255)
                            aaCAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTAc                                                             <  1:242344/91‑1 (MQ=255)
                            aaCAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTAc                                                             >  1:3719491/1‑91 (MQ=255)
                             aCAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACa                                                            >  1:2262308/1‑91 (MQ=255)
                             aCAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACa                                                            >  1:3262002/1‑91 (MQ=255)
                              cAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACAt                                                           <  1:604343/91‑1 (MQ=255)
                              cAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACAt                                                           >  1:2390242/1‑91 (MQ=255)
                              cAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACAt                                                           <  1:228376/91‑1 (MQ=255)
                                aaaaTCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTa                                                         >  1:3105645/1‑91 (MQ=255)
                                aaaaTCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTa                                                         >  1:3105606/1‑91 (MQ=255)
                                  aaTCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATc                                                       <  1:2874832/91‑1 (MQ=255)
                                     ccTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCttt                                                    <  1:636830/91‑1 (MQ=255)
                                     ccTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCCATTTGTACACATTTTTACAAATACGCTACATTATCttt                                                    <  1:3826394/91‑1 (MQ=255)
                                      cTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTg                                                   >  1:1094435/1‑91 (MQ=255)
                                      cTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTg                                                   <  1:4111455/91‑1 (MQ=255)
                                      cTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTg                                                   >  1:948221/1‑91 (MQ=255)
                                       tCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGt                                                  >  1:1954332/1‑91 (MQ=255)
                                       tCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGt                                                  >  1:3133208/1‑91 (MQ=255)
                                        cacaTCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTc                                                 >  1:763980/1‑91 (MQ=255)
                                           atcatGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGAt                                              <  1:2688764/91‑1 (MQ=255)
                                           atcatGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGAt                                              <  1:2688810/91‑1 (MQ=255)
                                            tcatGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATa                                             <  1:536119/91‑1 (MQ=255)
                                                gAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTa                                         >  1:756311/1‑91 (MQ=255)
                                                  aTTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAAc                                       <  1:540488/91‑1 (MQ=255)
                                                  aTTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAAc                                       <  1:2585639/91‑1 (MQ=255)
                                                   ttAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACa                                      <  1:2640737/91‑1 (MQ=255)
                                                       gATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACAttct                                  <  1:3647190/91‑1 (MQ=255)
                                                       gATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACAttct                                  <  1:165020/91‑1 (MQ=255)
                                                        aTGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACAttctt                                 <  1:2026229/91‑1 (MQ=255)
                                                        aTGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACAttctt                                 <  1:1534255/91‑1 (MQ=255)
                                                             aGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTcc                            <  1:746027/91‑1 (MQ=255)
                                                              ggAGAAAAATTTATCGTTAAATAAGCCTATTTGCACACATTTTTACAAATACGCTACATGATCTTTGTCGATAATTAACATTCTTTCTccc                           <  1:3682434/91‑1 (MQ=255)
                                                              ggACAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTccc                           <  1:3264860/91‑1 (MQ=255)
                                                               gagaAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCa                          <  1:861903/91‑1 (MQ=255)
                                                               gagaAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCa                          >  1:3637263/1‑91 (MQ=255)
                                                               gagaAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTCTCGATAATTAACATTCTTTCTCCCa                          >  1:3849004/1‑91 (MQ=255)
                                                               gagaAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTCTCGATAATTAACATTCTTTCTCCCa                          >  1:3850014/1‑91 (MQ=255)
                                                                 gaAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATc                        <  1:1027088/91‑1 (MQ=255)
                                                                 gaAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATc                        <  1:95218/91‑1 (MQ=255)
                                                                 gaAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATc                        <  1:1027767/91‑1 (MQ=255)
                                                                 gaAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATc                        >  1:3827207/1‑91 (MQ=255)
                                                                  aaaaaTTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATcc                       <  1:12334/91‑1 (MQ=255)
                                                                       tttATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACt                  <  1:595091/91‑1 (MQ=255)
                                                                        ttATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACtt                 >  1:1182264/1‑91 (MQ=255)
                                                                         tATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACttt                >  1:2986696/1‑91 (MQ=255)
                                                                           tCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTaa              <  1:771182/91‑1 (MQ=255)
                                                                           tCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTaa              >  1:1963445/1‑91 (MQ=255)
                                                                           tCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTaa              >  1:1959271/1‑91 (MQ=255)
                                                                           tCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTaa              >  1:1799916/1‑91 (MQ=255)
                                                                           tCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTaa              >  1:3927457/1‑91 (MQ=255)
                                                                             gTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTAACt            >  1:579828/1‑91 (MQ=255)
                                                                               tAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTAACTGt          <  1:3110131/91‑1 (MQ=255)
                                                                                aaaTAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTAACTGTc         <  1:1364433/91‑1 (MQ=255)
                                                                                 atgtAGCGTATTTGTAAAAATGTGTACAAATAGGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTAACTGTCg        <  1:2264370‑M1/87‑1 (MQ=255)
                                                                                 aaTAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTAACTGTCg        <  1:1174239/91‑1 (MQ=255)
                                                                                 aaTAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTAACTGTCg        <  1:1174304/91‑1 (MQ=255)
                                                                                  tgtAGCGTATTTGTAAAAATGTGTACAAATAGGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTAACTGTCgg       <  1:1591477‑M1/88‑1 (MQ=255)
                                                                                  aTAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTAACTGTCgg       >  1:1012846/1‑91 (MQ=255)
                                                                                  aTAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTAACTGTCgg       >  1:4027066/1‑91 (MQ=255)
                                                                                    tAGCGTATTTGTAAAAATGTGTACAAATAGGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTAACTGTCGGct     <  1:3973540‑M1/90‑1 (MQ=255)
                                                                                      GCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTAACTGTCGGctct   <  1:1389977‑M1/91‑1 (MQ=255)
                                                                                      GCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTAACTGTCGGctct   >  1:607696‑M1/1‑91 (MQ=255)
                                                                                      GCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTAACTGTCGGctct   >  1:3249736‑M1/1‑91 (MQ=255)
                                                                                       CCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTAACTGTCGGCTCTa  <  1:4106426‑M1/91‑1 (MQ=255)
                                                                                        |                                                                                         
TAAATTGACTCAAATGATCACCTCTATAAACAAAAATCCTCACATCATGAATTAAGATGCAAGGAGAAAAATTTATCGTTAAATAAGCCTATTTGTACACATTTTTACAAATACGCTACATTATCTTTGTCGATAATTAACATTCTTTCTCCCATCCAGACTTTAACTGTCGGCTCTA  >  USA300TCH1516_ALE/1899009‑1899186

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