Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 235,092 C→T 6.5% D490N (GAT→AAT)  gsiA ← Glutathione import ATP‑binding protein GsiA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE235,0920CT 6.5% 196.1 / 4.1 71D490N (GAT→AAT) gsiAGlutathione import ATP‑binding protein GsiA
Reads supporting (aligned to +/- strand):  ref base C (32/34);  new base T (2/3);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

ACAACGCATGACTTTTATTTAGCAATAACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTCTAAAGCGGACACTGCCTCGTCGCAAACAA  >  USA300TCH1516_ALE/234953‑235232
                                                                                                                                           |                                                                                                                                            
acaACGCATGACTTTTATTTAGCAATAACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTCATAATTGCAACGCGATcac                                                                                                                                            <  1:497055/141‑1 (MQ=255)
 caACGCTTGACTTTTATTTAGCAATAACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATcac                                                                                                                                            <  2:287484/141‑1 (MQ=255)
 caACGCATGACTTTTATTTAGCAATAACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATcac                                                                                                                                            <  2:287508/141‑1 (MQ=255)
 caACGCATGACTTTTATTTAGCAATAACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATcac                                                                                                                                            <  2:822871/141‑1 (MQ=255)
  caCGCATGACTTTTATTTAGCAATAACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTCTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATcaca                                                                                                                                           <  2:71815/140‑1 (MQ=255)
  aaCGCATGACTTTTATTTAGCAATAACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATcaca                                                                                                                                           <  2:901908/141‑1 (MQ=255)
  aaCGCATGACTTTTATTTAGCAATAACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATcaca                                                                                                                                           <  2:674199/141‑1 (MQ=255)
  aaCGCATGACTTTTATTTAGCAATAACTGCTACTTCGGGAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATcaca                                                                                                                                           <  1:851373/141‑1 (MQ=255)
     gCATGACTTTTATTTAGCAATAACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGAGTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCACCGCGATCACatat                                                                                                                                        >  2:206214/1‑141 (MQ=255)
     gCATGACTTTTATTTAGCAATAACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATTACatat                                                                                                                                        >  1:774384/1‑141 (MQ=255)
     gCAAGACTTTTATTTAGCAATAACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATTCCatat                                                                                                                                        >  1:654618/1‑141 (MQ=255)
         gACTTTTATTTAGCAATAACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCa                                                                                                                                    <  1:147084/141‑1 (MQ=255)
            tttgATTTAGCAATTACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTg                                                                                                                                 <  1:385220/141‑1 (MQ=255)
                  ttaGCAATAACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGGTTTCGCCATTTTTCATAACTGCAACCCGATCACATATTTCATTGATAacc                                                                                                                           >  1:184466/1‑140 (MQ=255)
                  ttaGCAATAACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGGGATCACATTTTTCATTGATacac                                                                                                                           >  1:184500/1‑138 (MQ=255)
                  ttaGCAATAACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAaca                                                                                                                           >  1:673965/1‑141 (MQ=255)
                       caTAACTGCTACTTCTGAAATAAGTTGACTTGCATAGTCTGACTGCGGATGTTTGATAATATCCTCTGGGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCcat                                                                                                                      <  1:615256/140‑1 (MQ=255)
                          aaCTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCcatgtc                                                                                                                   <  2:18880/141‑1 (MQ=255)
                             tGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCcatgtcatg                                                                                                                >  1:381089/1‑141 (MQ=255)
                               cTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCAtgtg                                                                                                              <  1:415252/141‑1 (MQ=255)
                                tACTTCTGAAAGAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACGGCAACGCGATCACATATTTCATTGATAACACCCATGTCAtgtgt                                                                                                             <  2:702721/141‑1 (MQ=255)
                                  cttcttAAATAATTTGCTTTGCATAGTCTGACTGCTGGTGTTTGATGATATCTTCTGTGTTATTCAGTTCAACGATTTCACCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGtga                                                                                                           <  2:897086/141‑1 (MQ=255)
                                   ttttGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGtgat                                                                                                          <  2:691776/138‑1 (MQ=255)
                                      cgAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGaa                                                                                                       <  1:533102/140‑1 (MQ=255)
                                       gAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCAAATATTTCATTGATAACACCCATGTCATGTGTGATGaat                                                                                                      <  2:327466/141‑1 (MQ=255)
                                          aTAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGgataaa                                                                                                   >  1:483636/1‑141 (MQ=255)
                                           tAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCTTTTTTTTATCGCGATCACATATTTCATTGATAACACCCATGTGATGTGGGATGaataaaa                                                                                                  >  2:568687/1‑140 (MQ=255)
                                           tAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGaataaat                                                                                                  >  1:753155/1‑141 (MQ=255)
                                            aaGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGaataaata                                                                                                 >  1:805553/1‑141 (MQ=255)
                                            aaGTAGCTTTGCATAGTCTGACTGCGGATGTGTGATAATATCTTCGGTGTTATTGAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGAGCACACATTTCATTGATAACACCCATGTCATGTGTGAGGaagaaata                                                                                                 >  2:766481/1‑141 (MQ=255)
                                               ttGCTTTGCATAGCCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGt                                                                                              <  2:482403/141‑1 (MQ=255)
                                                tGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTg                                                                                             >  1:400877/1‑141 (MQ=255)
                                                 gCTTGGCATAGTCTGACTGCGGATGTTTTATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGa                                                                                            <  2:136773/141‑1 (MQ=255)
                                                  ttttGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGAt                                                                                           <  2:74755/140‑1 (MQ=255)
                                                      gCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGGGATGCCg                                                                                       >  2:778730/1‑141 (MQ=255)
                                                               gACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTTTGATGAATAAATAAGTGATGCCGAAGTCTGAc                                                                              >  1:65651/1‑141 (MQ=255)
                                                                    cGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTaa                                                                         >  2:559798/1‑141 (MQ=255)
                                                                          tttGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGttt                                                                   >  2:171101/1‑141 (MQ=255)
                                                                             gATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTtaa                                                                >  2:606348/1‑141 (MQ=255)
                                                                               taataTCTTCTTTGTTATTCAGTTCAACGATTTCGCCAATTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTtaata                                                              <  2:402253/141‑1 (MQ=255)
                                                                               taataTCTTCTGCGTTATTCAGTTCAACGATTTCTCCTTGTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAAGTGTTTtaata                                                              <  2:341591/141‑1 (MQ=255)
                                                                                      ttctGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCgat                                                       >  2:851809/1‑141 (MQ=255)
                                                                                        tttttttATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTCATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAGCTCgatga                                                     <  2:660779/136‑1 (MQ=255)
                                                                                        ttgtgTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCgatga                                                     <  1:769383/140‑1 (MQ=255)
                                                                                        ctGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCgatga                                                     >  2:826673/1‑141 (MQ=255)
                                                                                        ctGTGTTATTCAGTTCAACGATTTCGCCATTTTACATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCgatga                                                     <  1:619254/141‑1 (MQ=255)
                                                                                         tgtgTTATTCAGTTCACCGATTTCTCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCgatgat                                                    <  2:476891/141‑1 (MQ=255)
                                                                                         tgtgTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATTACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCgatgat                                                    <  2:654618/141‑1 (MQ=255)
                                                                                         tgtgTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATTACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCgatgat                                                    <  2:774384/141‑1 (MQ=255)
                                                                                         tgtgTTATTCAGTTCAACGAGTTCGCCATTTTTCATAACTGCAACGCGACCACATATGTCATTGATAACACCCATGTCCTGTGTGATGAATAAATAAGTGCTGCCGAAGCCTAACTGTAATTGTTTTATTAACCCgatata                                                    >  2:815329/1‑138 (MQ=255)
                                                                                          gtgtTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGata                                                   <  2:685157/141‑1 (MQ=255)
                                                                                           tctgcttCACTTCAACGATGTCGCCAGTTTTCATAACTCCAACTCGATCACATATTTCTTTGATATCACGCCTGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTTATTGTTTTAATTACTCGATGatat                                                  <  2:523551/136‑1 (MQ=255)
                                                                                            gtTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATc                                                 >  2:246813/1‑141 (MQ=255)
                                                                                             ttATTCAGTTCAACGCTTTCGCCATTTTTTATAACTGCAACGCGATCACATATTTCATTGATAACCCCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCt                                                <  1:469456/141‑1 (MQ=255)
                                                                                                     ttCAACGCCTTCGCCATTTTTCATAACTGCAACGTGATCACATATTTCCTTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATTCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTttgaat                                        <  1:526782/141‑1 (MQ=255)
                                                                                                      tCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTttgaatt                                       <  1:511459/141‑1 (MQ=255)
                                                                                                          cGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTGTTAATAACTCGATGATATCTTTTTGAATTGaaa                                   <  2:713066/141‑1 (MQ=255)
                                                                                                           gATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAAc                                  >  2:361390/1‑141 (MQ=255)
                                                                                                            aTTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACg                                 >  1:271992/1‑141 (MQ=255)
                                                                                                              ttCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATACCACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTc                               >  2:100601/1‑141 (MQ=255)
                                                                                                              ttCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTc                               >  1:915577/1‑141 (MQ=255)
                                                                                                              tgcGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTc                               <  1:676315/139‑1 (MQ=255)
                                                                                                               tCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTCt                              <  2:63437/141‑1 (MQ=255)
                                                                                                                 gCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTGTaa                            >  1:380743/1‑141 (MQ=255)
                                                                                                                  ccATTTTTCATAACTGCAACGAGCACCCATAGTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAGCGTCTaga                           >  1:218362/1‑139 (MQ=255)
                                                                                                                       tttCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTCTAAAGCGGa                      >  1:721690/1‑141 (MQ=255)
                                                                                                                       tttCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGCGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTCTAAAGCGGa                      >  2:549013/1‑141 (MQ=255)
                                                                                                                       tttCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGGGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTATAATAACTCTACGATATCTTTTTGAATTGAACCGTCTAAACCggc                      >  1:390815/1‑140 (MQ=255)
                                                                                                                          cATAACTGCAACGGGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTACTTGTTTTGCTACCAGAAGGGGAACGTTTTGAATTGAAACGTCTAAAGCAGacac                   >  2:639218/1‑141 (MQ=255)
                                                                                                                           aTAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTCTAAAGCGGACACt                  >  1:183481/1‑141 (MQ=255)
                                                                                                                           aTAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTCTAAAGCGGACACt                  <  2:442436/141‑1 (MQ=255)
                                                                                                                           aTAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTCTAAAGCGGACACt                  <  1:107680/141‑1 (MQ=255)
                                                                                                                               cTGCAACGCGATTACACAGTTCATTGATAACGCCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTCTAAAGCGGACACTGCCt              <  1:741765/141‑1 (MQ=255)
                                                                                                                               cTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTCTAAAGCGGACACTGCCt              <  1:642327/141‑1 (MQ=255)
                                                                                                                                tGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTCTAAAGCGGACACTGCCtc             <  1:143070/141‑1 (MQ=255)
                                                                                                                                  cAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATTACGGTATCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTCTAAAGCGGACACTGCCtcgt           <  1:431417/141‑1 (MQ=255)
                                                                                                                                  cAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTCTAAAACGGACACTGCCtcgt           <  1:676001/141‑1 (MQ=255)
                                                                                                                                   aaCGCGCTCACATATGTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAATCGTCTAAAGGGGACACTGCCtcgtc          <  1:454514/141‑1 (MQ=255)
                                                                                                                                       cgATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTCTAAAGCGGACACTGCCTCGTCGCaa      >  2:614285/1‑141 (MQ=255)
                                                                                                                                           cacaTATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTCTAAAGCGGACACTGCCTCGTCGCAAACaa  <  1:578005/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                            
ACAACGCATGACTTTTATTTAGCAATAACTGCTACTTCTGAAATAAGTTGCTTTGCATAGTCTGACTGCGGATGTTTGATAATATCTTCTGTGTTATTCAGTTCAACGATTTCGCCATTTTTCATAACTGCAACGCGATCACATATTTCATTGATAACACCCATGTCATGTGTGATGAATAAATAAGTGATGCCGAAGTCTAACTGTAATTGTTTTAATAACTCGATGATATCTTTTTGAATTGAAACGTCTAAAGCGGACACTGCCTCGTCGCAAACAA  >  USA300TCH1516_ALE/234953‑235232

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