Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,608,816 A→G 6.3% intergenic (+97/+163) USA300TCH1516_02485 → / ← USA300HOU_RS13450 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,608,8160AG 6.3% 668.7 / 26.2 221intergenic (+97/+163)USA300TCH1516_02485/USA300HOU_RS13450hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (110/97);  new base G (6/8);  total (116/105)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.90e-01

Only 100 of 240 total aligned reads displayed.

TTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCAC  >  USA300TCH1516_ALE/2608726‑2608906
                                                                                          |                                                                                          
ttCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCcaca                                                                                            >  2:2394143/1‑91 (MQ=255)
ttCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCcaca                                                                                            <  2:3006057/91‑1 (MQ=255)
    tttATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCa                                                                                        >  1:2271095/1‑91 (MQ=255)
     ttATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCaa                                                                                       <  1:1923819/91‑1 (MQ=255)
        tatataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAatat                                                                                    >  2:3623500/1‑91 (MQ=255)
        tatataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAatat                                                                                    <  2:3551966/91‑1 (MQ=255)
         atataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAAtata                                                                                   >  1:973733/1‑91 (MQ=255)
            taGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGa                                                                                >  1:544511/1‑91 (MQ=255)
            taGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTCTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGa                                                                                >  1:3872168/1‑91 (MQ=255)
             ggCAGACATGATAGAATGTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGAc                                                                               <  2:1383877/90‑1 (MQ=255)
             aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTATATTGGat                                                                               <  1:1704224/91‑2 (MQ=255)
             aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTATATTGGat                                                                               <  1:3684213/91‑2 (MQ=255)
              gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                              <  2:3180675/91‑1 (MQ=255)
              gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                              <  1:1653847/91‑1 (MQ=255)
              gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                              >  1:60708/1‑91 (MQ=255)
              gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                              <  1:2818040/91‑1 (MQ=255)
              gCAGACATGAAAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                              <  2:1581968/91‑1 (MQ=255)
               cAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGt                                                                             <  1:3155626/91‑1 (MQ=255)
                 gACATGATAGAATTTTATATGTAAAGCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTgg                                                                           <  2:892480/91‑1 (MQ=255)
                   cATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCt                                                                         >  2:2035585/1‑91 (MQ=255)
                      gATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCtttt                                                                      >  2:1281453/1‑91 (MQ=255)
                        tAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAg                                                                    <  2:658361/91‑1 (MQ=255)
                         aGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGt                                                                   >  1:707584/1‑91 (MQ=255)
                          gAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTa                                                                  >  1:1446832/1‑91 (MQ=255)
                              tttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGt                                                              >  2:2593238/1‑91 (MQ=255)
                               ttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGtt                                                             >  2:2669076/1‑91 (MQ=255)
                               ttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTCAGTATCGtt                                                             >  1:3709593/1‑91 (MQ=255)
                                tataTGTAAATCTTGTAGGTAATCGTTTTAAAGATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTc                                                            <  1:3308218/91‑1 (MQ=255)
                                  tatGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCa                                                          <  2:1079691/91‑1 (MQ=255)
                                   atGTAAATCTTGTAGGTATTCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAg                                                         >  2:2016601/1‑91 (MQ=255)
                                   atGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAg                                                         <  2:3810529/91‑1 (MQ=255)
                                   atGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAg                                                         >  1:640083/1‑91 (MQ=255)
                                       aaaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGAttt                                                     <  2:1123421/91‑1 (MQ=255)
                                       aaaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGAttt                                                     <  2:105057/91‑1 (MQ=255)
                                       aaaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGAttt                                                     <  2:2192361/91‑1 (MQ=255)
                                        aaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTtg                                                    <  1:821873/91‑1 (MQ=255)
                                           cTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgg                                                 >  2:1397770/1‑91 (MQ=255)
                                           cTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgg                                                 <  2:3248531/91‑1 (MQ=255)
                                             tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgggg                                               >  1:1035918/1‑91 (MQ=255)
                                             tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgggg                                               <  2:1655337/91‑1 (MQ=255)
                                             tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgggg                                               >  2:456911/1‑91 (MQ=255)
                                              gTAGGTAATCGTTTTAAAAATAATCTAAGTATGACTATAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGACTTGTGGGGa                                              >  2:2752759/1‑91 (MQ=255)
                                              gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGa                                              >  2:1826386/1‑91 (MQ=255)
                                              gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGa                                              >  1:2028591/1‑91 (MQ=255)
                                               tAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCTAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAg                                             >  2:3487731/1‑91 (MQ=255)
                                               tAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAg                                             <  2:2325856/91‑1 (MQ=255)
                                               tAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAg                                             >  2:1688930/1‑91 (MQ=255)
                                               tAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAg                                             >  2:1502744/1‑91 (MQ=255)
                                                 ggTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAAGATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGtt                                           <  2:500439/91‑1 (MQ=255)
                                                  gTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTa                                          <  2:3804551/91‑1 (MQ=255)
                                                    aaTCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAgg                                        <  1:1373913/91‑1 (MQ=255)
                                                     aTCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAggg                                       <  1:1471333/91‑1 (MQ=255)
                                                     aTCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAggg                                       <  1:1070588/91‑1 (MQ=255)
                                                      tCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGt                                      <  2:1766724/91‑1 (MQ=255)
                                                      tCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGt                                      >  1:2932831/1‑91 (MQ=255)
                                                       cGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGtt                                     >  2:3293163/1‑91 (MQ=255)
                                                         ttttAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTa                                   >  2:2274835/1‑91 (MQ=255)
                                                          ttttAAAATAATATAAGCATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTATGGTTTAc                                  >  2:3436056/1‑91 (MQ=255)
                                                           ttAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTTCAGATTTGTGGGGAGTTAGGGTTTACt                                 >  2:2914573/1‑91 (MQ=255)
                                                           ttAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACt                                 >  1:2085195/1‑91 (MQ=255)
                                                            tAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACtt                                >  2:3720055/1‑91 (MQ=255)
                                                              aaaaTAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaa                              >  2:3871432/1‑91 (MQ=255)
                                                              aaaaTAATATAAGCATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaa                              >  2:2365606/1‑91 (MQ=255)
                                                               aaaTAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaat                             <  2:3713542/91‑1 (MQ=255)
                                                                  taataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaattaa                          >  2:141643/1‑91 (MQ=255)
                                                                  taataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaattaa                          >  2:702824/1‑91 (MQ=255)
                                                                   aataTAAGTATGACTAAAGCCACAGCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTaaa                         <  1:3420417/91‑1 (MQ=255)
                                                                    atatAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAg                        >  1:1794600/1‑91 (MQ=255)
                                                                     tataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAgg                       >  1:1014982/1‑91 (MQ=255)
                                                                      ataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAggg                      >  2:2803351/1‑91 (MQ=255)
                                                                      ataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAggg                      >  1:2575014/1‑91 (MQ=255)
                                                                       taAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAgggg                     >  2:2030946/1‑91 (MQ=255)
                                                                       taAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAgggg                     >  1:515476/1‑91 (MQ=255)
                                                                        aaGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGt                    >  1:3846568/1‑91 (MQ=255)
                                                                        aaGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGt                    >  2:3018621/1‑91 (MQ=255)
                                                                         aGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGACTTGTGGGGAGTTAGGGTTTACTTAATTAAAggggcc                   >  1:442970/1‑89 (MQ=255)
                                                                          gTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCt                  <  2:251171/91‑1 (MQ=255)
                                                                          gTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCt                  >  1:1828776/1‑91 (MQ=255)
                                                                           tATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTg                 <  1:3204270/91‑1 (MQ=255)
                                                                            aTGACTAAGGCCGCATAGTATGTGGTACGTGGCTTTTAGTATCGTTCCAGGTTTGTGGGGAGTTAGGGTTTACTTAATCAAATGGGTCTGc                <  2:1279037/91‑1 (MQ=255)
                                                                             tGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGcc               <  1:2231050/91‑1 (MQ=255)
                                                                             tGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGcc               >  1:1752646/1‑91 (MQ=255)
                                                                              gACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGccc              <  2:1133380/91‑1 (MQ=255)
                                                                                cTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaatta                           <  2:1838468/76‑1 (MQ=255)
                                                                                 tAAAGCCACATCCAATATAGGACGTGGCTTTTGGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCa           <  2:2454122/91‑1 (MQ=255)
                                                                                 tAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGTCCTCa           >  2:3580776/1‑91 (MQ=255)
                                                                                 tAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCa           <  1:78952/91‑1 (MQ=255)
                                                                                    aGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAAtt        >  1:1885381/1‑91 (MQ=255)
                                                                                     gCCACGTCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg       >  1:3042426/1‑91 (MQ=38)
                                                                                     gCCACGTCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATCTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg       >  2:2833859/1‑91 (MQ=38)
                                                                                     gCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg       >  1:35404/1‑91 (MQ=255)
                                                                                     gCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg       <  2:794201/91‑1 (MQ=255)
                                                                                      ccACGTCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGc      >  2:242421/1‑91 (MQ=38)
                                                                                      ccACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGc      >  2:733577/1‑91 (MQ=255)
                                                                                        acaTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTc    >  2:2121872/1‑91 (MQ=255)
                                                                                        acaTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTc    <  1:2707087/91‑1 (MQ=255)
                                                                                        acaTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTc    >  1:1612388/1‑91 (MQ=255)
                                                                                         caTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCa   <  1:1360441/91‑1 (MQ=255)
                                                                                         caTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGCGGTCTGCCCTCAATTGCTCa   >  1:2183716/1‑91 (MQ=255)
                                                                                          aTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCAc  <  1:2632043/91‑1 (MQ=255)
                                                                                          |                                                                                          
TTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCAC  >  USA300TCH1516_ALE/2608726‑2608906

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: