Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 19,009 C→T 100% S222F (TCC→TTC)  gdpP → Cyclic‑di‑AMP phosphodiesterase GdpP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE19,0090CT99.4% 530.6 / ‑5.0 153S222F (TCC→TTC) gdpPCyclic‑di‑AMP phosphodiesterase GdpP
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (73/79);  total (74/79)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 153 total aligned reads displayed.

CGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGTTGGTGAAG  >  USA300TCH1516_ALE/18871‑19146
                                                                                                                                          |                                                                                                                                         
cGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCg                                                                                                                                         >  4:360560/1‑141 (MQ=255)
 gACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGa                                                                                                                                        >  2:1159835/1‑141 (MQ=255)
  aCATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGAt                                                                                                                                       <  3:86266/141‑1 (MQ=255)
  aCATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGAt                                                                                                                                       <  2:149499/141‑1 (MQ=255)
    attattTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCa                                                                                                                                     <  1:52774/141‑1 (MQ=255)
         ttttAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCg                                                                                                                                >  3:60739/1‑141 (MQ=255)
          attAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGt                                                                                                                               >  3:385087/2‑141 (MQ=255)
            tAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAg                                                                                                                             <  3:500189/141‑1 (MQ=255)
             aGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGc                                                                                                                            <  1:100961/141‑1 (MQ=255)
             aGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGc                                                                                                                            <  2:771666/141‑1 (MQ=255)
              gATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGcc                                                                                                                           >  1:585239/1‑141 (MQ=255)
              gATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGcc                                                                                                                           <  1:134183/141‑1 (MQ=255)
                  aCTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTAtt                                                                                                                       <  2:1120900/141‑1 (MQ=255)
                     aCGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTaa                                                                                                                    >  2:794833/1‑141 (MQ=255)
                        aTGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATc                                                                                                                 <  2:938927/141‑1 (MQ=255)
                           agaTTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCaaa                                                                                                              <  1:128233/141‑1 (MQ=255)
                            gaTTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTTGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCaaaa                                                                                                             <  3:364421/141‑1 (MQ=255)
                            gaTTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCaaaa                                                                                                             >  2:954925/1‑141 (MQ=255)
                            gaTTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCaaaa                                                                                                             <  1:495125/141‑1 (MQ=255)
                               tACGCAGAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAat                                                                                                          <  1:762244/141‑1 (MQ=255)
                                 cGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAatat                                                                                                        >  1:855905/1‑141 (MQ=255)
                                    aaaaTATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAg                                                                                                     <  2:368978/141‑1 (MQ=255)
                                     aaaTATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGc                                                                                                    >  3:359009/1‑141 (MQ=255)
                                     aaaTATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGc                                                                                                    <  2:178445/141‑1 (MQ=255)
                                     aaaTATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGc                                                                                                    <  1:968475/141‑1 (MQ=255)
                                      aaTATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAAGCAAAAAATATTAGCt                                                                                                   <  1:281293/141‑1 (MQ=255)
                                       atatGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTg                                                                                                  >  2:749227/1‑141 (MQ=255)
                                          tgaatgaTACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACt                                                                                               <  1:300342/141‑1 (MQ=255)
                                           gaatgaTACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCTATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACtt                                                                                              >  4:272611/1‑141 (MQ=255)
                                           gaatgaTACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACtt                                                                                              >  4:325350/1‑141 (MQ=255)
                                             atgaTACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTag                                                                                            >  1:633531/1‑141 (MQ=255)
                                             atgaTACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTag                                                                                            >  2:86155/1‑141 (MQ=255)
                                                 tACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTagaaga                                                                                        >  1:981917/1‑141 (MQ=255)
                                                    gCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATc                                                                                     <  2:849394/141‑1 (MQ=255)
                                                     cAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATACGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCt                                                                                    <  3:412623/141‑1 (MQ=255)
                                                          ttCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAAtt                                                                               <  3:169124/141‑1 (MQ=255)
                                                            cGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTg                                                                             <  1:450549/141‑1 (MQ=255)
                                                            cGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGAAAAAAAAAAAAAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTg                                                                             <  2:377748/141‑1 (MQ=255)
                                                               aaaTCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGata                                                                          >  2:1158221/1‑141 (MQ=255)
                                                                  tCAACTCAATGGTAACGCGTGTCATTAGTCGATGGCCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATGTGATATCt                                                                       >  1:1026518/1‑141 (MQ=255)
                                                                   cAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCtt                                                                      >  1:205099/1‑141 (MQ=255)
                                                                      cTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAg                                                                   >  2:235622/1‑141 (MQ=255)
                                                                      cTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGAATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAg                                                                   >  4:259769/1‑141 (MQ=255)
                                                                       tCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAgg                                                                  >  2:818496/1‑140 (MQ=255)
                                                                       tCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGt                                                                  >  4:438861/1‑141 (MQ=255)
                                                                       tCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGt                                                                  <  2:1036483/141‑1 (MQ=255)
                                                                           tGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAAt                                                              >  2:707281/1‑141 (MQ=255)
                                                                           tGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAAt                                                              >  1:1006725/1‑141 (MQ=255)
                                                                           tGGTAACGCGTGTCATTAGTCGATGGGCAACCTAGTCGAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAAt                                                              >  2:118651/1‑141 (MQ=255)
                                                                            ggTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAAtt                                                             >  4:426388/1‑141 (MQ=255)
                                                                              tAACGCGTGTCATTAGTCGATGGGAAACTGAGTATAATATATTTTTCAAAAGAAACAATATCAAAAAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTaaaaaaa                                                           <  2:209451/141‑5 (MQ=255)
                                                                               aaCGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACg                                                          <  4:58309/141‑1 (MQ=255)
                                                                               aaCGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGAAAATAAAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCaaaaaaa                                                          <  2:558763/141‑6 (MQ=255)
                                                                               aaCGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTAGATCAATTCGTAGACTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTaaaaaaaa                                                          <  2:666521/141‑6 (MQ=255)
                                                                                 cgcgTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTg                                                        >  2:926133/1‑141 (MQ=255)
                                                                                  gcgTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGa                                                       >  2:639416/1‑141 (MQ=255)
                                                                                    gtgtCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGaaa                                                     <  1:794833/141‑1 (MQ=255)
                                                                                    gtgtCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGaaa                                                     <  1:406485/141‑1 (MQ=255)
                                                                                    gtgtCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATAGTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGaaa                                                     >  1:24235/1‑141 (MQ=255)
                                                                                       tCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGaaaaaa                                                  <  4:359009/141‑1 (MQ=255)
                                                                                          ttAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAgtg                                               <  1:1061738/141‑1 (MQ=255)
                                                                                          ttAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAgtg                                               <  2:964073/141‑1 (MQ=255)
                                                                                           tAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAgtgt                                              >  1:1146441/1‑141 (MQ=255)
                                                                                             gTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTg                                            <  1:869022/141‑1 (MQ=255)
                                                                                             gTCGATGGGCAACTGAGTATAATATATTTTTAAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGTTATCTTGAGTAAAAAAAAAAAAAAAAGTGTTg                                            <  4:287723/141‑1 (MQ=255)
                                                                                               cGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGt                                          <  3:434352/141‑1 (MQ=255)
                                                                                               cGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCATATATAAATAAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGt                                          <  4:402635/141‑1 (MQ=255)
                                                                                                  tGGGCAACTGAGTATAATATATTTTTTTTTTTTTTTTGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGCAGAATCTAGATTTGATATCTTGAGTCAATTACGTGAAAATAGTGTTGGTTAt                                       >  2:97067/1‑141 (MQ=255)
                                                                                                  tGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTAt                                       >  3:190620/1‑141 (MQ=255)
                                                                                                  tGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTAt                                       <  1:430317/141‑1 (MQ=255)
                                                                                                   gggCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATc                                      >  1:544214/1‑141 (MQ=255)
                                                                                                    ggCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCg                                     <  3:280389/141‑1 (MQ=255)
                                                                                                    ggCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCg                                     <  4:106622/141‑1 (MQ=255)
                                                                                                       aaCTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGc                                  <  2:315315/141‑1 (MQ=255)
                                                                                                         cggAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGTTTATCGTGccc                                >  2:29538/3‑141 (MQ=255)
                                                                                                         cTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGCAAAAAGCGTTGGTTATCGTGccc                                >  2:41479/1‑141 (MQ=255)
                                                                                                         cTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGccc                                >  1:1095387/1‑141 (MQ=255)
                                                                                                            aGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAAt                             <  1:36122/141‑1 (MQ=255)
                                                                                                             gTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAAtt                            >  4:65570/1‑141 (MQ=255)
                                                                                                             gTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAAtt                            <  1:80052/141‑1 (MQ=255)
                                                                                                              tataATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTa                           >  1:470658/1‑141 (MQ=255)
                                                                                                                taataTATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAAc                         >  3:384982/1‑141 (MQ=255)
                                                                                                                taataTATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAAc                         >  3:310044/1‑141 (MQ=255)
                                                                                                                  atatatTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACAt                       <  3:20548/141‑1 (MQ=255)
                                                                                                                  atatatTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATATAAAAAAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACAt                       <  2:1098659/141‑1 (MQ=255)
                                                                                                                   tatatTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACAtt                      >  1:999802/1‑141 (MQ=255)
                                                                                                                    atatTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTa                     >  1:749086/1‑141 (MQ=255)
                                                                                                                    atatTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTa                     <  1:826699/141‑1 (MQ=255)
                                                                                                                    atatTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATATAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTa                     >  4:483725/1‑141 (MQ=255)
                                                                                                                      atTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTAGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAg                   >  1:701117/1‑141 (MQ=255)
                                                                                                                        ttttCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTa                 <  1:982210/141‑1 (MQ=255)
                                                                                                                           tCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCg              >  3:79221/1‑141 (MQ=255)
                                                                                                                            cAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCgg             >  1:500966/1‑141 (MQ=255)
                                                                                                                             aaaaGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGgt            <  4:388446/141‑1 (MQ=255)
                                                                                                                             aaaaGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGgt            <  3:498317/141‑1 (MQ=255)
                                                                                                                                aGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGtt         >  4:13505/1‑141 (MQ=255)
                                                                                                                                aGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGtt         >  2:142654/1‑141 (MQ=255)
                                                                                                                                 gATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGCTATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGTTg        <  1:86155/141‑1 (MQ=255)
                                                                                                                                 gATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGTTg        <  1:726412/141‑1 (MQ=255)
                                                                                                                                       gTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGTTGGTGAAg  >  3:115839/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                         
CGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGTTGGTGAAG  >  USA300TCH1516_ALE/18871‑19146

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