Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,963,768 C→G 12.9% V365L (GTA→CTA)  hemY ← Protoporphyrinogen oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,963,7680CG12.9% 197.0 / 20.7 70V365L (GTA→CTA) hemYProtoporphyrinogen oxidase
Reads supporting (aligned to +/- strand):  ref base C (29/32);  new base G (3/6);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

ATACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTTT  >  USA300TCH1516_ALE/1963629‑1963906
                                                                                                                                           |                                                                                                                                          
aTACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACa                                                                                                                                           >  1:496469/1‑141 (MQ=255)
aTACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTAATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACa                                                                                                                                           >  1:409752/1‑141 (MQ=255)
 tACTTTTCGGCAAGCGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGACACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACAt                                                                                                                                          <  2:34260/141‑1 (MQ=255)
    ttttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACg                                                                                                                                       <  1:462848/141‑1 (MQ=255)
      ttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCa                                                                                                                                     <  2:129578/141‑1 (MQ=255)
       tCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCAc                                                                                                                                    <  2:287282/141‑1 (MQ=255)
        cGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACg                                                                                                                                   <  2:220297/141‑1 (MQ=255)
         ggCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGTTACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGa                                                                                                                                  <  2:221487/141‑1 (MQ=255)
             aaCGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCa                                                                                                                              <  1:6261/141‑1 (MQ=255)
                  ttGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAacc                                                                                                                         >  2:227647/1‑141 (MQ=255)
                  ttGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATATCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAacc                                                                                                                         >  1:600434/1‑141 (MQ=255)
                     aCAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTATGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTCGTTTACCTACATACGCACGAATCAAAacctta                                                                                                                      >  1:408213/1‑141 (MQ=255)
                         ttGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAaccttacctt                                                                                                                  <  2:36084/141‑1 (MQ=255)
                         ttGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAaccttacctt                                                                                                                  >  2:228592/1‑141 (MQ=255)
                          tGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTc                                                                                                                 <  2:557356/141‑1 (MQ=255)
                               aTTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGgtg                                                                                                            >  1:423350/1‑141 (MQ=255)
                                ttCAGGATCCCCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGgtgt                                                                                                           <  1:511576/141‑1 (MQ=255)
                                    ggATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATCCTAATTCATTATCTGTATGATCGTCAACCACAGTACCACCTGGTTTACCTACATACGCCCGAATCCAACCCTTACCTGCCGctgacgtc                                                                                                       >  2:540858/1‑136 (MQ=255)
                                    ggATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAATATCACCTGGTTTACCTACATACGCACGAGTCAAAACCTTACCTTCCGGGgtatta                                                                                                       >  1:383423/1‑141 (MQ=255)
                                        cACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATg                                                                                                   <  2:496469/141‑1 (MQ=255)
                                                    tgatcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCg                                                                                       <  2:516243/139‑1 (MQ=255)
                                                    tcatcaTTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCg                                                                                       <  1:422161/141‑1 (MQ=255)
                                                    tcatcaTTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCg                                                                                       <  1:119519/141‑1 (MQ=255)
                                                    tcatcaTTTGACTTAAATCGCTACGTACAATCGATACTCATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTGCATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCg                                                                                       >  1:341703/1‑141 (MQ=255)
                                                     catcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACGTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGa                                                                                      >  1:443557/1‑141 (MQ=255)
                                                     catcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGa                                                                                      >  1:113614/1‑141 (MQ=255)
                                                     catcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGa                                                                                      >  1:178634/1‑141 (MQ=255)
                                                      atcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTCTGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGATGTAGTAAATGGCCATTTTTTCGAt                                                                                     <  2:238202/141‑1 (MQ=255)
                                                       tcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGTTg                                                                                    >  1:298057/1‑141 (MQ=255)
                                                        catGTGACTTCAATCTCGACGTACAATCGATGCTGAGTCATTATCTATATGAGCATCAACCACAGTATCACCTGGTTTACCTTCATACGCACAAATCCAAACCTTAGCGTCCGGTGGAGGAGAAGGCCAGTTTTTCgcggg                                                                                   >  2:347121/1‑137 (MQ=255)
                                                         atTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTc                                                                                  >  1:504656/1‑141 (MQ=255)
                                                            tGACTTAAATCTCTACGGACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACAGACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCaa                                                                               <  2:162212/141‑1 (MQ=255)
                                                                ttAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTaa                                                                           >  2:426975/1‑140 (MQ=255)
                                                                     tccctACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCg                                                                      <  1:213333/138‑1 (MQ=255)
                                                                      ctctACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCgg                                                                     <  1:150689/141‑1 (MQ=255)
                                                                       tctACGAACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGt                                                                    <  2:150699/141‑1 (MQ=255)
                                                                         tACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTaa                                                                  <  1:84902/141‑1 (MQ=255)
                                                                               cAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTg                                                            >  1:373274/1‑141 (MQ=255)
                                                                               cAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTg                                                            <  1:174406/141‑1 (MQ=255)
                                                                                 aTCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTa                                                          <  1:227647/141‑1 (MQ=255)
                                                                                    gATAATAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCa                                                       <  1:341253/141‑1 (MQ=255)
                                                                                        cTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCAATCg                                                   >  2:383352/1‑141 (MQ=255)
                                                                                        cTAATTCATTATCTGTACGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCg                                                   <  1:67190/141‑1 (MQ=255)
                                                                                         tAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGt                                                  <  1:333437/141‑1 (MQ=255)
                                                                                          aaTTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCCCTCGtt                                                 >  1:354255/1‑141 (MQ=255)
                                                                                           aTTCATTATCTCTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGTGCCAACACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTtc                                                <  1:598546/141‑1 (MQ=255)
                                                                                                 tATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATAACTCGTTCTCGCaa                                          >  1:55611/1‑141 (MQ=255)
                                                                                                    cTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACGTACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTa                                       >  1:436418/1‑141 (MQ=255)
                                                                                                    cTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATAGCCATTTGTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTa                                       >  2:526858/1‑141 (MQ=255)
                                                                                                       tATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCAGGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGa                                    <  2:232140/141‑1 (MQ=255)
                                                                                                           atcatcAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGcc                                >  2:538930/1‑141 (MQ=255)
                                                                                                           atcatcAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGcc                                >  2:538315/1‑141 (MQ=255)
                                                                                                            tcatcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCa                               <  1:437373/141‑1 (MQ=255)
                                                                                                            tcatcaACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCCCTCGTTCTCGCAATTACGAAGCCa                               >  2:381628/1‑141 (MQ=255)
                                                                                                              ctcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGt                             <  1:86671/140‑1 (MQ=255)
                                                                                                              atcaACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTATGTATCACTCGTTCTCGCACTTACGCAGCCAGt                             >  1:46307/1‑141 (MQ=255)
                                                                                                                ccaCCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTAc                           <  1:208328/139‑1 (MQ=255)
                                                                                                                caACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTAc                           >  1:233986/1‑141 (MQ=255)
                                                                                                                  aCCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACca                         >  2:512778/1‑141 (MQ=255)
                                                                                                                    cacaGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatc                       <  1:531685/141‑1 (MQ=255)
                                                                                                                     acaGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatca                      <  2:381773/141‑1 (MQ=255)
                                                                                                                     acaGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatca                      <  1:299844/141‑1 (MQ=255)
                                                                                                                     aaaGTATCACCTGGTTTACCTACATACTCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatca                      <  2:523363/139‑1 (MQ=255)
                                                                                                                          aTCACCTGGTTCACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGt                 <  2:339141/141‑1 (MQ=255)
                                                                                                                              ccTGGTTTACCTACATTCGCACGAATCATGTCCTTACCGTCCTGTGTAGTAAATGGCCACTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTAttt             <  1:14182/141‑1 (MQ=255)
                                                                                                                              ccTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCACCATGGGTAttt             >  1:65426/1‑141 (MQ=255)
                                                                                                                               cTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTAtttt            >  2:491482/1‑141 (MQ=255)
                                                                                                                                  gTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCaa         <  2:233986/141‑1 (MQ=255)
                                                                                                                                      aCCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTc     <  2:113614/141‑1 (MQ=255)
                                                                                                                                         tACATACGCACGAATCAAAACCTTACCATACTCTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCttt  >  1:145766/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                          
ATACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTTT  >  USA300TCH1516_ALE/1963629‑1963906

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: