Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,239,562 A→G 5.7% intergenic (+29/‑161) 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,5620AG 5.7% 463.3 / 12.9 139intergenic (+29/‑161)USA300HOU_RS06160/recGhypothetical protein/ATP‑dependent DNA helicase RecG
Reads supporting (aligned to +/- strand):  ref base A (88/43);  new base G (5/3);  total (93/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 150 total aligned reads displayed.

TAGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTC  >  USA300TCH1516_ALE/1239472‑1239651
                                                                                           |                                                                                         
tAGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATa                                                                                           <  1:1849116/91‑1 (MQ=255)
tAGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATa                                                                                           <  1:1431337/91‑1 (MQ=255)
 aGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATaa                                                                                          >  1:3737195/1‑91 (MQ=255)
 aGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATaa                                                                                          >  1:2910631/1‑91 (MQ=255)
 aGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATaa                                                                                          >  1:1447468/1‑91 (MQ=255)
  gAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATaaa                                                                                         >  1:792997/1‑91 (MQ=255)
  gAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATaaa                                                                                         >  1:3475610/1‑91 (MQ=255)
  gAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTATTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATaaa                                                                                         >  1:3296932/1‑91 (MQ=255)
    aGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatc                                                                                       >  1:707567/1‑91 (MQ=255)
         aaGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAg                                                                                  >  1:1678458/1‑91 (MQ=255)
          aGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGt                                                                                 <  1:3602796/91‑1 (MQ=255)
            ttCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTg                                                                               >  1:2754709/1‑91 (MQ=255)
            ttCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTg                                                                               >  1:1719081/1‑91 (MQ=255)
              cATGAAGGTGGACAACCAATTTATCAATATTTCTTTTGAGTAGAATAAAAATTTAAAATAAAAAAGTACCAATGATAAATCATCAGTTGGt                                                                             >  1:3864014/1‑91 (MQ=255)
              cATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGt                                                                             >  1:2413178/1‑91 (MQ=255)
               aTGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTa                                                                            >  1:2365198/1‑91 (MQ=255)
                tGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAg                                                                           >  1:704084/1‑91 (MQ=255)
                 gAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGt                                                                          <  1:728717/91‑1 (MQ=255)
                 gAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGt                                                                          >  1:671248/1‑91 (MQ=255)
                 gAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGt                                                                          >  1:437858/1‑91 (MQ=255)
                 gAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGt                                                                          >  1:1429424/1‑91 (MQ=255)
                   aGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGttt                                                                        <  1:3229851/91‑1 (MQ=255)
                   aGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGttt                                                                        >  1:119871/1‑91 (MQ=255)
                   aGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGttt                                                                        >  1:268504/1‑91 (MQ=255)
                    ggTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGtttt                                                                       >  1:1313631/1‑91 (MQ=255)
                       ggACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTa                                                                    >  1:3687339/1‑91 (MQ=255)
                        gACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAt                                                                   <  1:1996739/91‑1 (MQ=18)
                          cAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAttt                                                                 >  1:202198/1‑91 (MQ=255)
                          cAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAttt                                                                 >  1:2206013/1‑91 (MQ=255)
                          cAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAttt                                                                 >  1:3140166/1‑91 (MQ=255)
                          cAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAttt                                                                 >  1:3140334/1‑91 (MQ=255)
                           aaCCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAtttt                                                                >  1:540937/1‑91 (MQ=255)
                            aCCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTg                                                               >  1:1120606/1‑91 (MQ=255)
                             ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAttttaa                                                              <  1:631220/91‑3 (MQ=17)
                             ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                              <  1:3866583/91‑1 (MQ=255)
                             ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                              <  1:2421881/91‑1 (MQ=255)
                             ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                              <  1:1022722/91‑1 (MQ=255)
                             ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                              <  1:689577/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:2110216/91‑1 (MQ=255)
                                     tCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAg                                                      >  1:616493/1‑91 (MQ=255)
                                     tCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAg                                                      >  1:3049741/1‑91 (MQ=255)
                                      cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                     >  1:3143422/1‑91 (MQ=255)
                                      cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                     >  1:847651/1‑91 (MQ=255)
                                      cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                     >  1:847415/1‑91 (MQ=255)
                                      cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                     <  1:1200458/91‑1 (MQ=255)
                                      cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                     >  1:2017972/1‑91 (MQ=255)
                                      cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                     >  1:2018163/1‑91 (MQ=255)
                                       aaTATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTg                                                    >  1:1661706/1‑91 (MQ=255)
                                        atatTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGa                                                   <  1:3359968/91‑1 (MQ=255)
                                            ttCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATAtt                                               <  1:1248143/91‑1 (MQ=255)
                                              cTTTTCAGTAGAATAAAAATTTAAAATAAAAACCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGc                                             <  1:3536855/91‑1 (MQ=255)
                                               ttttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCg                                            >  1:151469/1‑91 (MQ=255)
                                                 ttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCggg                                          >  1:4033071/1‑91 (MQ=255)
                                                  tCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGt                                         <  1:420693/91‑1 (MQ=255)
                                                   cAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGtt                                        <  1:1467679/91‑1 (MQ=255)
                                                     gTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTaa                                      >  1:2818825/1‑91 (MQ=255)
                                                        gAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAg                                   >  1:3946994/1‑91 (MQ=255)
                                                         gaTAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGt                                  >  1:3249236/2‑91 (MQ=255)
                                                         aaTAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGt                                  <  1:1290008/91‑1 (MQ=255)
                                                           tAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTAt                                >  1:1925046/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:2737954/91‑1 (MQ=255)
                                                                  ttAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTtctc                         >  1:817509/1‑91 (MQ=255)
                                                                  ttAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTtctc                         <  1:737885/91‑1 (MQ=255)
                                                                   taaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCg                        <  1:2686399/91‑1 (MQ=255)
                                                                   taaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCg                        >  1:472664/1‑91 (MQ=255)
                                                                    aaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGa                       >  1:920858/1‑91 (MQ=255)
                                                                    aaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGa                       >  1:533639/1‑91 (MQ=255)
                                                                    aaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGa                       >  1:533608/1‑91 (MQ=255)
                                                                     aaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAg                      >  1:3734291/1‑91 (MQ=255)
                                                                     aaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAg                      >  1:473140/1‑91 (MQ=255)
                                                                     aaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAg                      >  1:3925124/1‑91 (MQ=255)
                                                                      aataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGt                     >  1:1146502/1‑91 (MQ=255)
                                                                      aataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGt                     >  1:3472271/1‑91 (MQ=255)
                                                                       ataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGtt                    >  1:208192/1‑91 (MQ=255)
                                                                       ataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGtt                    >  1:572843/1‑91 (MQ=255)
                                                                         aaaaaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGc                  >  1:1056004/1‑91 (MQ=18)
                                                                          aaaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCt                 <  1:2513344/91‑1 (MQ=255)
                                                                          aaaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCt                 >  1:3721250/1‑91 (MQ=255)
                                                                          aaaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCt                 >  1:3721274/1‑91 (MQ=255)
                                                                            aaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa               >  1:494789/1‑91 (MQ=18)
                                                                            aaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa               <  1:2046591/91‑1 (MQ=255)
                                                                             aaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAAc              <  1:1121031/91‑1 (MQ=255)
                                                                               cTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACaa            >  1:639258/1‑91 (MQ=255)
                                                                                tACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAAt           <  1:3107292/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)
                                                                                        gATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGt    <  1:2167517/91‑1 (MQ=255)
                                                                                        gATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGt    >  1:417438/1‑91 (MQ=255)
                                                                                          tAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTc  >  1:1245402/1‑91 (MQ=255)
                                                                                          tAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTc  <  1:3202942/91‑1 (MQ=255)
                                                                                           |                                                                                         
TAGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTC  >  USA300TCH1516_ALE/1239472‑1239651

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: