Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,825,587 C→T 5.2% intergenic (‑402/+446) cap8A_2 ← / ← icaR Capsular polysaccharide type 8 biosynthesis protein cap8A/Biofilm operon icaADBC HTH‑type negative transcriptional regulator IcaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,825,5870CT 5.2% 838.2 / 21.6 251intergenic (‑402/+446)cap8A_2/icaRCapsular polysaccharide type 8 biosynthesis protein cap8A/Biofilm operon icaADBC HTH‑type negative transcriptional regulator IcaR
Reads supporting (aligned to +/- strand):  ref base C (121/117);  new base T (3/10);  total (124/127)
Fisher's exact test for biased strand distribution p-value = 8.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

Only 100 of 267 total aligned reads displayed.

CTTTACAACACCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATA  >  USA300TCH1516_ALE/2825497‑2825677
                                                                                          |                                                                                          
ctTTACAACACCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCacac                                                                                            >  2:3625656/1‑91 (MQ=255)
ctTTACAACACCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCacac                                                                                            <  1:737568/91‑1 (MQ=255)
  ttACAACACCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCt                                                                                          <  1:1104049/91‑1 (MQ=255)
  ttACAACACCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCt                                                                                          <  1:629876/91‑1 (MQ=255)
   tACAACACCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCtt                                                                                         <  1:2234169/91‑1 (MQ=255)
    acaacaCCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCttt                                                                                        <  1:397806/91‑1 (MQ=255)
     caacaCCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCtttt                                                                                       <  1:1989110/91‑1 (MQ=255)
       acacCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTaa                                                                                     >  1:379110/1‑91 (MQ=255)
        cacCGAGCGTGCCTGTACTCGGTAATTCAATACCTTGCGTAACGCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTaaa                                                                                    <  1:3962339/91‑1 (MQ=255)
        cacCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTaaa                                                                                    <  2:1641282/91‑1 (MQ=255)
        cacCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTaaa                                                                                    <  1:2374552/91‑1 (MQ=255)
        cacCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTaaa                                                                                    <  2:1044449/91‑1 (MQ=255)
         acCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACATCTTTTtaaa                                                                                   <  1:3565424/91‑1 (MQ=255)
         acCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTaaaa                                                                                   >  2:3116146/1‑91 (MQ=255)
          ccGAGCGTGCCCGTACTCGGTAATTCAGTACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACATCTTTTtaaaa                                                                                  <  1:1494965/91‑1 (MQ=255)
          ccGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTaaaaa                                                                                  >  2:1169004/1‑91 (MQ=255)
            gAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGa                                                                                >  2:3254929/1‑91 (MQ=255)
             aGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGAt                                                                               <  2:245570/91‑1 (MQ=255)
              gCGTGCCCGTACTCGGTAATTCAATACCTGGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGAtg                                                                              <  2:2604185/91‑1 (MQ=255)
               cGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACATCTTTTTAAAAGgtgt                                                                             <  2:1142566/91‑1 (MQ=255)
               cGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGAtgt                                                                             >  2:435920/1‑91 (MQ=255)
                gTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGAtgtg                                                                            <  2:961675/91‑1 (MQ=255)
                 tGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTgg                                                                           >  1:3718057/1‑91 (MQ=255)
                 tGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTgg                                                                           >  2:2893892/1‑91 (MQ=255)
                   cccGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGt                                                                         >  1:164655/1‑91 (MQ=255)
                   cccGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTCTAAAAAGATGTGGGt                                                                         >  2:784788/1‑91 (MQ=255)
                    ccGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTa                                                                        >  2:329552/1‑91 (MQ=255)
                       tACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAAtt                                                                     <  2:2568642/91‑1 (MQ=255)
                         cTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTta                                                                   >  2:872393/1‑91 (MQ=255)
                         cTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTta                                                                   >  1:3082301/1‑91 (MQ=255)
                          tCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTtat                                                                  >  1:2528959/1‑91 (MQ=255)
                          tCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTtat                                                                  <  1:1600685/91‑1 (MQ=255)
                           cGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTtata                                                                 <  1:1139592/91‑1 (MQ=255)
                                aTTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAAtt                                                            <  2:2172265/91‑1 (MQ=255)
                                aTTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAAtt                                                            <  1:3158421/91‑1 (MQ=255)
                                  tCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAAtttt                                                          <  1:2112894/91‑1 (MQ=255)
                                   cAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAAttttt                                                         >  1:63642/1‑91 (MQ=255)
                                   cAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAAttttt                                                         >  1:3519248/1‑91 (MQ=255)
                                    aaTACCTTGCGCAACCCGTCACAGTGAGTTGGGTTAATGATAATAAAGCGCACACCTTTGAAAAAGATGTGGGTAATTTATATAATTttta                                                        <  2:3011178/91‑1 (MQ=255)
                                     aTACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttat                                                       >  1:457432/1‑91 (MQ=255)
                                     aTACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttat                                                       >  1:1804491/1‑91 (MQ=255)
                                     aTACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttat                                                       <  1:943769/91‑1 (MQ=255)
                                     aTACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttat                                                       <  1:3177580/91‑1 (MQ=255)
                                     aTACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttat                                                       >  2:1796817/1‑91 (MQ=255)
                                      tACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttatt                                                      >  1:258299/1‑91 (MQ=255)
                                       aCCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttattt                                                     <  2:919237/91‑1 (MQ=255)
                                        ccTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttattta                                                    <  1:1469144/91‑1 (MQ=255)
                                        ccTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttattta                                                    <  2:2722004/91‑1 (MQ=255)
                                            gCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACAtt                                                >  2:2643421/1‑91 (MQ=255)
                                                aaCCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTTGGTAATTTATATAATTTTTATTTACATTTTTa                                            >  1:3353900/1‑91 (MQ=255)
                                                  cccGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTCTTTACATTTTTAAc                                          >  1:1191242/1‑91 (MQ=255)
                                                   ccGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACt                                         >  2:1833057/1‑91 (MQ=255)
                                                    cGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTACTTTATATAATTTTTATTTACATTTTTAACtt                                        >  2:2230159/1‑91 (MQ=255)
                                                     gTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACGTTTTTAACTta                                       >  1:3587092/1‑91 (MQ=255)
                                                      tCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTtat                                      <  1:174102/91‑1 (MQ=255)
                                                      tCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTtat                                      >  1:1056965/1‑91 (MQ=255)
                                                      tCACTGGGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTtat                                      <  1:1435552/91‑1 (MQ=255)
                                                       cACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTtata                                     >  1:2767654/1‑91 (MQ=255)
                                                        aCTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATaa                                    <  1:3329155/91‑1 (MQ=255)
                                                        aCTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATaa                                    >  2:3765688/1‑91 (MQ=255)
                                                        aCTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATaa                                    <  1:2434551/91‑1 (MQ=255)
                                                         cTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATaaa                                   <  1:2049550/91‑1 (MQ=255)
                                                          tgtgAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATaaaa                                  >  1:2182909/1‑91 (MQ=255)
                                                           gtgAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATaaaca                                 >  2:834535/1‑89 (MQ=255)
                                                           gtgAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATaaaaa                                 >  1:2278236/1‑91 (MQ=255)
                                                           gtgAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTAAATATTTAACTTATaaaaa                                 >  2:837568/1‑91 (MQ=255)
                                                             gAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATaaaaaaa                               <  1:1954726/91‑1 (MQ=255)
                                                               gTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAgc                             <  2:3826131/91‑1 (MQ=255)
                                                               gTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAgc                             >  1:3130902/1‑91 (MQ=255)
                                                               gTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAgc                             >  2:1011476/1‑91 (MQ=255)
                                                                ttGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAgcg                            <  1:998079/91‑1 (MQ=255)
                                                                ttGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTAGGTAATTTATATAATTTTTATTTACATTTTAAAAAAAAAAAAAAAgcg                            <  2:3914222/91‑1 (MQ=255)
                                                                 tGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAgcgc                           >  2:666959/1‑91 (MQ=255)
                                                                  gggTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGcc                          <  1:1276760/91‑1 (MQ=255)
                                                                   ggTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCt                         >  1:2567700/1‑91 (MQ=255)
                                                                    gTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTa                        <  2:3388953/91‑1 (MQ=255)
                                                                     ttAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTAt                       >  1:160534/1‑91 (MQ=255)
                                                                       aaTGATAATAAAGCCCACACCTTTTAAAACGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGAGCCTATGt                     >  2:2166257/1‑91 (MQ=255)
                                                                       aaTGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATACTTTTTCTTTACATTTTTAACTTATAAAAAAAGCGCCTAtac                     >  2:3241161/1‑89 (MQ=255)
                                                                       aaTGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGt                     >  1:3840244/1‑91 (MQ=255)
                                                                        aTGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTTTTTACATTTTTAACTTATAAAAAAAGCGCCTATGTc                    >  2:1209918/1‑91 (MQ=255)
                                                                         tGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCa                   <  1:2752431/91‑1 (MQ=255)
                                                                            taataaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAACAAAGCGCCTATGTCATGa                >  2:2117670/1‑91 (MQ=255)
                                                                             aataaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGAt               >  1:820154/1‑91 (MQ=255)
                                                                             aataaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGAt               <  2:2064632/91‑1 (MQ=255)
                                                                             aataaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGAt               <  1:3931361/91‑1 (MQ=255)
                                                                              ataaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGAtt              <  2:1358252/91‑1 (MQ=255)
                                                                               taaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGAttt             >  1:2937661/1‑91 (MQ=255)
                                                                               taaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGAttt             <  1:2892420/91‑1 (MQ=255)
                                                                                 aaGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTAc           >  2:3523854/1‑91 (MQ=255)
                                                                                 aaGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTAc           >  2:1310138/1‑91 (MQ=255)
                                                                                   gCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCa         >  2:2308834/1‑91 (MQ=255)
                                                                                    cccACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCAt        <  1:1716107/91‑1 (MQ=255)
                                                                                    cccACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCAt        <  1:2643749/91‑1 (MQ=255)
                                                                                     ccACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATc       >  1:2773968/1‑91 (MQ=255)
                                                                                       acacCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATcac     >  1:2101466/1‑91 (MQ=255)
                                                                                       acacCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATcac     <  2:2864701/91‑1 (MQ=255)
                                                                                        cacCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATcaca    >  1:1250596/1‑91 (MQ=255)
                                                                                        cacCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATcaca    <  1:203117/91‑1 (MQ=255)
                                                                                          ccTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATa  >  1:811257/1‑91 (MQ=255)
                                                                                          |                                                                                          
CTTTACAACACCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATA  >  USA300TCH1516_ALE/2825497‑2825677

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