Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,239,566 C→T 7.8% intergenic (+33/‑157) USA300HOU_RS06160 → / → recG hypothetical protein/ATP‑dependent DNA helicase RecG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,239,5660CT 7.8% 733.5 / 31.1 227intergenic (+33/‑157)USA300HOU_RS06160/recGhypothetical protein/ATP‑dependent DNA helicase RecG
Reads supporting (aligned to +/- strand):  ref base C (110/99);  new base T (9/9);  total (119/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

AGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTCAACTT  >  USA300TCH1516_ALE/1239476‑1239656
                                                                                           |                                                                                          
aGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatc                                                                                            >  2:3042142/1‑91 (MQ=255)
aGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatc                                                                                            >  1:711325/1‑91 (MQ=255)
aGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatc                                                                                            >  2:14538/1‑91 (MQ=255)
aGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatc                                                                                            >  2:1616792/1‑91 (MQ=255)
     aaGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAg                                                                                       >  2:1906368/1‑91 (MQ=255)
      aGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGt                                                                                      >  1:3171533/1‑91 (MQ=255)
       gTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGtt                                                                                     <  2:1702709/91‑1 (MQ=255)
       gTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGtt                                                                                     >  2:2126667/1‑91 (MQ=255)
         tCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTgg                                                                                   <  1:1450705/91‑1 (MQ=255)
          cATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGt                                                                                  >  2:3041361/1‑91 (MQ=255)
           aTGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTa                                                                                 >  1:3807260/1‑91 (MQ=255)
           aTGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTa                                                                                 >  2:536903/1‑91 (MQ=255)
           aTGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTa                                                                                 >  1:293052/1‑91 (MQ=255)
           aTGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTa                                                                                 >  1:2522200/1‑91 (MQ=255)
             gAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGt                                                                               <  2:2494804/91‑1 (MQ=255)
             gAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGt                                                                               <  1:1268199/91‑1 (MQ=255)
              aaGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGtt                                                                              <  1:2565319/91‑1 (MQ=21)
              aaGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGtt                                                                              >  1:996159/1‑91 (MQ=255)
                ggTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGtttt                                                                            <  2:1880289/91‑1 (MQ=255)
                ggTGGACAACCAATTTATCAATATTTCTTTTCAGCAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGtttt                                                                            >  2:3396286/1‑91 (MQ=255)
                 ttGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGttttt                                                                           <  2:3853740/90‑1 (MQ=18)
                 gTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGttttt                                                                           <  2:321157/91‑1 (MQ=255)
                  tGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGtttttt                                                                          >  1:62901/1‑91 (MQ=255)
                    gACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAt                                                                        <  2:2634454/91‑1 (MQ=18)
                     aCAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAtt                                                                       <  1:1072607/91‑1 (MQ=17)
                       aaCCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAtttt                                                                     <  1:1594614/91‑1 (MQ=255)
                        aCCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAtttta                                                                    <  2:2739780/91‑2 (MQ=17)
                        aCCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTg                                                                    <  1:962872/91‑1 (MQ=255)
                         ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                                   <  1:968107/91‑1 (MQ=255)
                         ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                                   >  1:3334634/1‑91 (MQ=255)
                         ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                                   <  2:3342831/91‑1 (MQ=255)
                          cAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCt                                                                  <  2:1791674/91‑1 (MQ=255)
                           aaTTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTa                                                                 >  2:1735672/1‑91 (MQ=255)
                             tttATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTAtt                                                               <  2:3276818/91‑1 (MQ=255)
                                aTCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTa                                                            >  1:3462370/1‑91 (MQ=255)
                                  cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTTCTATTTTAGt                                                          >  2:2595991/1‑91 (MQ=255)
                                  cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                          >  1:219229/1‑91 (MQ=255)
                                  cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                          <  2:826824/91‑1 (MQ=255)
                                   aaTATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTg                                                         <  1:2934297/91‑1 (MQ=255)
                                   aaTATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTg                                                         <  1:573887/91‑1 (MQ=255)
                                   aaTATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTg                                                         >  1:2775223/1‑91 (MQ=255)
                                     tatTTCTTTTCAGTAGAATAAAGATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTATTTTATTTTGCTATTTTAGTGat                                                       <  2:3512055/91‑1 (MQ=255)
                                         tCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTg                                                   >  2:3492851/1‑91 (MQ=255)
                                         tCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTg                                                   >  1:3471070/1‑91 (MQ=255)
                                         tCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTg                                                   >  1:185153/1‑91 (MQ=255)
                                           ttttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCg                                                 >  1:544064/1‑91 (MQ=255)
                                            tttCAGTAGAATAAAAGTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGctg                                                <  2:3712769/91‑3 (MQ=255)
                                            tttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCgg                                                >  2:477638/1‑91 (MQ=255)
                                             ttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCggg                                               >  1:3696421/1‑91 (MQ=255)
                                             ttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCggg                                               <  2:2585099/91‑1 (MQ=255)
                                              tCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGt                                              <  1:3266353/91‑1 (MQ=255)
                                                 gTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTaa                                           >  2:2008079/1‑91 (MQ=255)
                                                   aGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTaaaa                                         >  1:1960384/1‑91 (MQ=255)
                                                   aGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTaaaa                                         <  2:2047336/91‑1 (MQ=255)
                                                   aGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTaaaa                                         <  2:2314266/91‑1 (MQ=255)
                                                    gAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAg                                        >  1:2265623/1‑91 (MQ=255)
                                                       tAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTAt                                     <  1:3380951/91‑1 (MQ=255)
                                                        aaaaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATc                                    <  2:119237/91‑1 (MQ=255)
                                                        aaaaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATc                                    <  1:2685235/91‑1 (MQ=255)
                                                        aaaaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATc                                    <  1:2474201/91‑1 (MQ=255)
                                                         aaaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCg                                   >  2:364143/1‑91 (MQ=255)
                                                         aaaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCg                                   >  1:341898/1‑91 (MQ=255)
                                                          aaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGt                                  >  1:3143263/1‑91 (MQ=255)
                                                          aaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGt                                  <  1:3736419/91‑1 (MQ=255)
                                                           aaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGtt                                 >  2:926744/1‑91 (MQ=255)
                                                           aaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGtt                                 <  2:2009361/91‑1 (MQ=255)
                                                           aaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCCGGCTAAAAGTATCGtt                                 >  2:2976825/1‑91 (MQ=255)
                                                             tttAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTtct                               >  2:3107932/1‑91 (MQ=255)
                                                                 aaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAg                           >  2:2364634/1‑91 (MQ=255)
                                                                  aataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGt                          >  2:1875402/1‑91 (MQ=255)
                                                                   ataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGtt                         <  1:1722802/91‑1 (MQ=255)
                                                                      aaaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCt                      >  1:457647/1‑91 (MQ=255)
                                                                       aaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTa                     >  1:2034142/1‑91 (MQ=255)
                                                                        aaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa                    >  1:3657240/1‑91 (MQ=18)
                                                                        aaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa                    >  1:3142961/1‑91 (MQ=18)
                                                                        aaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa                    >  1:1645083/1‑91 (MQ=21)
                                                                        aaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa                    >  2:2303288/1‑91 (MQ=18)
                                                                        aaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa                    <  2:86586/91‑1 (MQ=255)
                                                                        aaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa                    >  2:276254/1‑91 (MQ=255)
                                                                         aaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAAc                   >  1:419838/1‑91 (MQ=255)
                                                                         aaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAAc                   >  1:1713417/1‑91 (MQ=255)
                                                                         aaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAAc                   >  1:761564/1‑91 (MQ=255)
                                                                          aCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACa                  >  2:115424/1‑91 (MQ=255)
                                                                          aCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACa                  <  1:2769600/91‑1 (MQ=255)
                                                                          aCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACa                  <  1:3214498/91‑1 (MQ=255)
                                                                           cTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACaa                 >  2:2847506/1‑91 (MQ=255)
                                                                           cTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACaa                 >  2:164844/1‑91 (MQ=255)
                                                                           cTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACaa                 >  1:3253683/1‑91 (MQ=255)
                                                                            tgccAATGATAAATCATCAGTTGGGAGTTTTTTATTTTGCTATTTTAGTGTTATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAAt                <  2:3599512/89‑1 (MQ=255)
                                                                             aCCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGATCAAAGTATCGTTCTCGAGTTGCTAACAATg               >  1:3107108/1‑91 (MQ=255)
                                                                                   tGATAAATCGTCAGTTGGTAGTTTTTTATTTTCCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATg          >  2:3805866/1‑91 (MQ=255)
                                                                                   tGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATg          <  1:1882178/91‑1 (MQ=255)
                                                                                     aTAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGtt        <  2:2672053/91‑1 (MQ=255)
                                                                                     aTAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGtt        <  2:1380683/91‑1 (MQ=255)
                                                                                      tAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTc       <  2:130295/91‑1 (MQ=255)
                                                                                        aaTCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTCaa     <  1:1188027/91‑1 (MQ=255)
                                                                                        aaTCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTCaa     >  2:2405017/1‑91 (MQ=255)
                                                                                         atcatcAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTCAAc    >  1:434649/1‑91 (MQ=255)
                                                                                          tcatcaGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTCAACt   >  2:441807/1‑91 (MQ=255)
                                                                                           catcaGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTCAACtt  >  1:658615/1‑91 (MQ=255)
                                                                                           |                                                                                          
AGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTCAACTT  >  USA300TCH1516_ALE/1239476‑1239656

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: