Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,239,557:1 +C 5.4% intergenic (+24/‑166) 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,5571.C 5.4% 667.0 / 29.0 149intergenic (+24/‑166)USA300HOU_RS06160/recGhypothetical protein/ATP‑dependent DNA helicase RecG
Reads supporting (aligned to +/- strand):  ref base . (96/45);  new base C (5/3);  total (101/48)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 153 total aligned reads displayed.

GATGTAGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTC  >  USA300TCH1516_ALE/1239468‑1239645
                                                                                          |                                                                                        
gATGTAGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAAt                                                                                         >  1:3331534/1‑91 (MQ=255)
gATGTAGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAAt                                                                                         >  1:1326608/1‑91 (MQ=255)
  tGTAGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGa                                                                                       <  1:4055908/91‑1 (MQ=255)
   gTAGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGAt                                                                                      >  1:2899516/1‑91 (MQ=255)
   gTAGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGAt                                                                                      >  1:3638675/1‑91 (MQ=255)
   gTAGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGAt                                                                                      <  1:2046069/91‑1 (MQ=255)
    tAGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATa                                                                                     <  1:1849116/91‑1 (MQ=255)
    tAGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATa                                                                                     <  1:1431337/91‑1 (MQ=255)
     aGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATaa                                                                                    >  1:3737195/1‑91 (MQ=255)
     aGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATaa                                                                                    >  1:764883/1‑91 (MQ=255)
     aGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATaa                                                                                    >  1:2910631/1‑91 (MQ=255)
      gAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATaaa                                                                                   >  1:3475610/1‑91 (MQ=255)
      gAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATaaa                                                                                   >  1:712511/1‑91 (MQ=255)
        aGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatc                                                                                 >  1:707567/1‑91 (MQ=255)
           ggAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatcatc                                                                              >  1:1626649/1‑91 (MQ=255)
             aaGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAg                                                                            >  1:1678458/1‑91 (MQ=255)
              aGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGt                                                                           <  1:3602796/91‑1 (MQ=255)
               gTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGtt                                                                          <  1:1265527/91‑1 (MQ=255)
                ttCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTg                                                                         >  1:2754709/1‑91 (MQ=255)
                 tCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTgg                                                                        >  1:1315925/1‑91 (MQ=255)
                  cATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGt                                                                       >  1:2413178/1‑91 (MQ=255)
                   aTGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTa                                                                      >  1:2365198/1‑91 (MQ=255)
                     gAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGt                                                                    <  1:1454606/91‑1 (MQ=255)
                     gAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGt                                                                    >  1:437858/1‑91 (MQ=255)
                     gAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGt                                                                    >  1:1429516/1‑91 (MQ=255)
                      aaGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGtt                                                                   <  1:218390/91‑1 (MQ=255)
                       aGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGttt                                                                  >  1:268504/1‑91 (MQ=255)
                       aGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGttt                                                                  <  1:3229851/91‑1 (MQ=255)
                       aGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGttt                                                                  >  1:119871/1‑91 (MQ=255)
                        ggTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGtttt                                                                 >  1:1313631/1‑91 (MQ=255)
                           ggACAACCAATTTATCAATATTTCTTTTCAGCAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTa                                                              >  1:3687503/1‑91 (MQ=255)
                            gACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAt                                                             <  1:1996739/91‑1 (MQ=18)
                              cAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAttt                                                           >  1:3140334/1‑91 (MQ=255)
                              cAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAttt                                                           >  1:3140166/1‑91 (MQ=255)
                               aaCCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAtttt                                                          >  1:540937/1‑91 (MQ=255)
                                aCCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTg                                                         >  1:1120606/1‑91 (MQ=255)
                                 ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAttttaa                                                        <  1:631331/91‑3 (MQ=17)
                                 ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAttttaa                                                        <  1:631220/91‑3 (MQ=17)
                                 ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                        <  1:2421881/91‑1 (MQ=255)
                                 ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                        <  1:3866583/91‑1 (MQ=255)
                                 ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                        <  1:689577/91‑1 (MQ=255)
                                 ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                        <  1:689652/91‑1 (MQ=255)
                                 ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                        <  1:1022722/91‑1 (MQ=255)
                                  cAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCt                                                       >  1:2226571/1‑91 (MQ=255)
                                     tttATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTCTGCTAtt                                                    >  1:2082718/1‑91 (MQ=255)
                                       tATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTAtttt                                                  >  1:2235962/1‑91 (MQ=255)
                                        aTCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTa                                                 >  1:647488/1‑91 (MQ=255)
                                         tCAATATTTCTTTTCAGTAGAATAAAAATTTAACATACAAAACTACCAATGATAAATCACCAGTTGGTAGTTTTTTATTTTGCTATTTTAg                                                >  1:3305490/1‑91 (MQ=255)
                                         tCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAg                                                >  1:373030/1‑91 (MQ=255)
                                         tCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAg                                                <  1:2110216/91‑1 (MQ=255)
                                         tCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAg                                                >  1:3049741/1‑91 (MQ=255)
                                         tCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAg                                                >  1:558430/1‑91 (MQ=255)
                                          cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                               >  1:3982496/1‑91 (MQ=255)
                                          cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                               >  1:3143422/1‑91 (MQ=255)
                                          cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                               >  1:847651/1‑91 (MQ=255)
                                          cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                               <  1:1200458/91‑1 (MQ=255)
                                          cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                               >  1:2018392/1‑91 (MQ=255)
                                          cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                               >  1:2018163/1‑91 (MQ=255)
                                           aaTATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTg                                              >  1:1917732/1‑91 (MQ=255)
                                            atatTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGa                                             <  1:3359968/91‑1 (MQ=255)
                                                ttCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATAtt                                         <  1:1248143/91‑1 (MQ=255)
                                                  cTTTTCCGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGc                                       >  1:3231847/1‑91 (MQ=255)
                                                  cTTTTCAGTAGAATAAAAATTTAAAATAAAAACCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGc                                       <  1:3536855/91‑1 (MQ=255)
                                                   ttttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCg                                      >  1:151469/1‑91 (MQ=255)
                                                         gTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTaa                                >  1:3874838/1‑91 (MQ=255)
                                                         gTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTaa                                >  1:2818825/1‑91 (MQ=255)
                                                            gAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAg                             >  1:3946994/1‑91 (MQ=255)
                                                                aaaaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATc                         >  1:970445/1‑91 (MQ=255)
                                                                aaaaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATc                         <  1:1158159/91‑1 (MQ=255)
                                                                  aaaCTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTGTCGt                       <  1:3270860/91‑1 (MQ=255)
                                                                   aaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGtt                      <  1:610184/91‑1 (MQ=255)
                                                                   aaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGtt                      <  1:2737954/91‑1 (MQ=255)
                                                                      ttAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTtctc                   >  1:817582/1‑91 (MQ=255)
                                                                      ttAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTtctc                   >  1:817509/1‑91 (MQ=255)
                                                                      ttAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTtctc                   <  1:737885/91‑1 (MQ=255)
                                                                       taaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCg                  <  1:2686399/91‑1 (MQ=255)
                                                                        aaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGa                 >  1:920734/1‑91 (MQ=255)
                                                                        aaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGa                 >  1:533608/1‑91 (MQ=255)
                                                                        aaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGa                 >  1:533506/1‑91 (MQ=255)
                                                                         aaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAg                >  1:473140/1‑91 (MQ=255)
                                                                         aaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAg                >  1:473042/1‑91 (MQ=255)
                                                                          aataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGt               >  1:1146502/1‑91 (MQ=255)
                                                                          aataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGt               >  1:3472271/1‑91 (MQ=255)
                                                                           ataaaaAACTACCAACTGATGATCTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGtt              >  1:3288385/1‑91 (MQ=18)
                                                                           ataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGtt              >  1:208192/1‑91 (MQ=255)
                                                                             aaaaaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGc            >  1:1056004/1‑91 (MQ=18)
                                                                             aaaaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGc            <  1:3090731/91‑1 (MQ=255)
                                                                              aaaaaCTACCACTGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCt           <  1:3764534/91‑1 (MQ=255)
                                                                              aaaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCt           >  1:655976/1‑91 (MQ=255)
                                                                              aaaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCt           >  1:1716161/1‑91 (MQ=255)
                                                                              aaaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCt           <  1:2513344/91‑1 (MQ=255)
                                                                              aaaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCt           >  1:3721250/1‑91 (MQ=255)
                                                                               aaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTa          >  1:428530/1‑91 (MQ=255)
                                                                                aaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa         >  1:2004606/1‑91 (MQ=18)
                                                                                aaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa         <  1:2046591/91‑1 (MQ=255)
                                                                                 aaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAAc        <  1:1121031/91‑1 (MQ=255)
                                                                                    tACCAATGATAAAGCATCAGTTGGTAGTTTTTTATTTTGCTATGTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAAt     <  1:4150861/91‑1 (MQ=255)
                                                                                     aCCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATg    >  1:1375042/1‑91 (MQ=255)
                                                                                      ccAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGt   <  1:2494706/91‑1 (MQ=255)
                                                                                       cAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTc  <  1:2786600/91‑1 (MQ=255)
                                                                                          |                                                                                        
GATGTAGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTC  >  USA300TCH1516_ALE/1239468‑1239645

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