Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 643,620 G→A 8.0% intergenic (‑268/+281) yjjP_1 ← / ← USA300HOU_RS03070 Inner membrane protein YjjP/Putative heme‑dependent peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE643,6200GA 8.0% 208.3 / 6.1 72intergenic (‑268/+281)yjjP_1/USA300HOU_RS03070Inner membrane protein YjjP/Putative heme‑dependent peroxidase
Reads supporting (aligned to +/- strand):  ref base G (29/37);  new base A (2/4);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGGAATTTTTAGTAGTTAAAATTCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACATAT  >  USA300TCH1516_ALE/643483‑643757
                                                                                                                                         |                                                                                                                                         
atcGGAATTTTTAGTAGTTAAAATTCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTa                                                                                                                                         <  1:118663/140‑1 (MQ=255)
 tcGGAATTATTAGTAGTTAAAATTCACTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAg                                                                                                                                        <  2:676745/140‑1 (MQ=255)
    gAATTTTTAGTAGTTAAAATTCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTTGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGa                                                                                                                                     >  1:357982/1‑140 (MQ=255)
      aTTTTTAGTAGTTAAAATTCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACt                                                                                                                                   <  1:583881/140‑1 (MQ=255)
          ttAGTAGTTAAAATTCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAGTATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGtat                                                                                                                               >  2:214270/1‑140 (MQ=255)
          ttAGTAGTTAAAATTCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTCTAGAGTAGCGACTGtat                                                                                                                               >  2:211266/1‑140 (MQ=255)
          ttAGTAGTTAAAATTCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAAAGTAGCGCCTGtat                                                                                                                               >  2:732149/1‑140 (MQ=255)
              tagtTAAAATTCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAAtt                                                                                                                           <  1:398671/140‑1 (MQ=255)
                   aaaaTTCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTAt                                                                                                                      <  2:471412/140‑1 (MQ=255)
                     aaTTCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGTGTATCGACTGTATAATTTCTATTg                                                                                                                    <  1:420166/140‑1 (MQ=255)
                     aaTTCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTg                                                                                                                    >  2:317507/1‑140 (MQ=255)
                       tgctctTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAg                                                                                                                  <  2:731878/138‑1 (MQ=255)
                        tctcTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAgg                                                                                                                 <  2:524975/140‑1 (MQ=255)
                        tctcTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAgg                                                                                                                 <  1:669494/140‑1 (MQ=255)
                         tcctggAAGTATAAAAAGGAAAAATTATAGAAAAAAACTCGACAAAAATATGTATCAAATATTTCAAATAGGTTAAAATAATATTAGAAAAAAAAAAAAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGt                                                                                                                <  2:699438/135‑1 (MQ=255)
                         ctctTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGt                                                                                                                <  2:370547/140‑1 (MQ=255)
                         ctctTGAAGAATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGt                                                                                                                <  2:8725/140‑1 (MQ=255)
                          tctTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGtt                                                                                                               >  1:645827/1‑140 (MQ=255)
                            ttGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTaa                                                                                                             >  1:758496/1‑140 (MQ=255)
                            ttGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTaa                                                                                                             <  2:649663/140‑1 (MQ=255)
                                aGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGtt                                                                                                         <  1:467614/140‑1 (MQ=255)
                                aGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAAGCTATATTGAGGACAACGtt                                                                                                         >  2:495968/1‑140 (MQ=255)
                                  tataAAAAGAAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTta                                                                                                       <  2:404998/140‑1 (MQ=255)
                                    taAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTtata                                                                                                     >  1:695451/1‑140 (MQ=255)
                                       aaaGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGt                                                                                                  <  1:635782/140‑1 (MQ=255)
                                       aaaGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTACTATTGAGGTTAACGTTTATATGt                                                                                                  >  2:3705/1‑140 (MQ=255)
                                       aaaGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAAGATTTCAAAGAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGt                                                                                                  <  1:80244/140‑1 (MQ=255)
                                             aaaaTTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATAAGTAGTGAt                                                                                            >  2:203752/1‑140 (MQ=255)
                                              aaaTTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACAATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTATAGTGAGGTTAACGTATATATGTAGTGAta                                                                                           >  2:125186/1‑140 (MQ=255)
                                                attattGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTATATATGAAGTGAtagt                                                                                         >  2:372474/1‑140 (MQ=255)
                                                    ttGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGAtagtatgt                                                                                     >  1:239434/1‑137 (MQ=255)
                                                      gAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTaa                                                                                   <  1:22120/140‑1 (MQ=255)
                                                            aaTTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTtct                                                                             <  1:373597/140‑1 (MQ=255)
                                                              ttcaataaaaaTATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGGATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATCTAGTGATAGTAGTTAAAGTTCTcc                                                                           <  2:297722/134‑1 (MQ=255)
                                                                  aCAAAAATATGATTCAAATATTTTAAATAGGTTAAAATAATAGTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAg                                                                       >  2:402167/1‑140 (MQ=255)
                                                                   cAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAgg                                                                      >  2:595583/1‑140 (MQ=255)
                                                                   cAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAATAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAgg                                                                      >  2:16356/1‑140 (MQ=255)
                                                                        atatGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGac                                                                 <  1:534281/140‑1 (MQ=255)
                                                                         tatGTTACAAATATTTCAAATAGGTTAAAATAAAATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGact                                                                <  2:568528/140‑1 (MQ=255)
                                                                             tttcaaatatttcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTc                                                            >  1:717194/1‑140 (MQ=255)
                                                                               tcaaatatttcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCgg                                                          >  1:684119/1‑140 (MQ=255)
                                                                                caaatatttcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAATAGAGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCggg                                                         >  2:713448/1‑140 (MQ=255)
                                                                                 acatatttcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGt                                                        <  2:640726/138‑1 (MQ=255)
                                                                                 aagtatttcaaataGGTTAAAAGAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGt                                                        <  1:195807/137‑1 (MQ=255)
                                                                                 aaatatttcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAATAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTTATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGt                                                        <  2:164810/140‑1 (MQ=255)
                                                                                  aatatttcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTa                                                       <  2:175565/140‑1 (MQ=255)
                                                                                   atatttcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGTAAGACTACTCGGGTac                                                      >  2:558315/1‑140 (MQ=255)
                                                                                        tcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCACGGAAGACTACTCGGGTACACttt                                                 >  1:101349/1‑140 (MQ=255)
                                                                                        tcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACttt                                                 >  1:39057/1‑140 (MQ=255)
                                                                                          aaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGc                                               <  1:179984/140‑1 (MQ=255)
                                                                                           aataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAATAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCt                                              <  2:581557/140‑1 (MQ=255)
                                                                                               ggTTAAAATAATATTAGTTAAGCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCGATTGAGGTTAACATTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGa                                          >  2:378586/1‑140 (MQ=255)
                                                                                                 ttAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGc                                        >  1:373564/1‑140 (MQ=255)
                                                                                                  tAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCa                                       >  2:719656/1‑140 (MQ=255)
                                                                                                     aataatATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAg                                    >  2:390682/1‑140 (MQ=255)
                                                                                                            ttAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACtt                             <  1:595583/140‑1 (MQ=255)
                                                                                                            ttAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACtt                             <  1:247438/140‑1 (MQ=255)
                                                                                                            ttAGTTAACCATTACAAAAATTGTATAGAATAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACtt                             <  1:494251/140‑1 (MQ=255)
                                                                                                             tAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACttt                            >  2:227589/1‑140 (MQ=255)
                                                                                                             tAGTTAACCACTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACttt                            <  1:94723/140‑1 (MQ=255)
                                                                                                                  aaCCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTTATAGTTGTTTTTGTTCTCGCAATGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTAt                       <  1:211266/140‑1 (MQ=255)
                                                                                                                  aaCCATTACAAAAATTGTATAGAATAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTAt                       <  1:16356/140‑1 (MQ=255)
                                                                                                                  aaCCATTACAAAAATGGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTAt                       <  1:214270/140‑1 (MQ=255)
                                                                                                                        tACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGataa                 >  1:318399/1‑140 (MQ=255)
                                                                                                                         aCAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGataat                <  2:443096/140‑1 (MQ=255)
                                                                                                                         aCAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGataat                >  1:9211/1‑140 (MQ=255)
                                                                                                                         aCAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGataat                <  1:203752/140‑1 (MQ=255)
                                                                                                                           aaaaaTTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATAc              <  2:10056/140‑1 (MQ=255)
                                                                                                                            aagatggTAGAGAGGAGCGACGGTATAAGTTCTATTGAGGTGAACGTTGATATGTAGGGATAGTAGTTAAAGTTCGCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACa             <  1:307794/134‑1 (MQ=255)
                                                                                                                             aaaTTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACAt            <  1:298790/140‑1 (MQ=255)
                                                                                                                                  gTATAGAGTAGCGACTGTATACCTACAATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGca       >  2:546337/1‑140 (MQ=255)
                                                                                                                                    ataGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGcaca     <  1:475031/140‑1 (MQ=255)
                                                                                                                                       gagTAGAGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACatat  <  2:641273/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                         
ATCGGAATTTTTAGTAGTTAAAATTCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACATAT  >  USA300TCH1516_ALE/643483‑643757

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