Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 491,927 T→A 9.9% L118L (CTT→CTA sle1_1 → N‑acetylmuramoyl‑L‑alanine amidase sle1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE491,9270TA 9.9% 243.6 / 2.5 98L118L (CTT→CTAsle1_1N‑acetylmuramoyl‑L‑alanine amidase sle1
Reads supporting (aligned to +/- strand):  ref base T (37/49);  new base A (10/2);  total (47/51)
Fisher's exact test for biased strand distribution p-value = 1.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCAAGTAATTCTACGAGTAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCATACTATACAGTACAAGCAGG  >  USA300TCH1516_ALE/491789‑492067
                                                                                                                                          |                                                                                                                                             
cTCAAGTAATTCTACGAGTAATAGTAGCCGTCCATCAACGAACCCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTaa                                                                                                                                             <  2:181275/141‑1 (MQ=255)
    aGTAATTCTACGAGTAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGt                                                                                                                                         >  1:1105776/1‑141 (MQ=255)
      tAATTCTACGAGTAATAGTAGCAGTCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGCCTAAATGGctt                                                                                                                                       >  2:41722/1‑141 (MQ=255)
       aaTTCTACGTGTAATAGTCGCCGTCCCTCAACGATCTCAGCTGTCGGATCATAGTACACAGTACAATCAGGCGGCGCATTATGATTAATCGCATCAAAATATGTTACAACTTACCAAAACATTATGCGACTTAATGGTTTa                                                                                                                                      <  2:921424/141‑1 (MQ=255)
         ttCTACGAGTAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTaaa                                                                                                                                    <  1:823921/141‑1 (MQ=255)
           cTACGAGTAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCTACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAaata                                                                                                                                  <  2:431313/141‑1 (MQ=255)
           cTACGAGTAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGCTACAACTTACCAAAACATTATGCGACTTAATGGTTTAaata                                                                                                                                  <  1:863917/141‑1 (MQ=255)
              cGAGTAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAAAATTATGCGACTTAATGGTTTAAATAAtt                                                                                                                               >  2:787692/1‑141 (MQ=255)
                   aaTAGTAGCCTGCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCATGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCttt                                                                                                                          <  2:333831/141‑1 (MQ=255)
                   aaTAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAACATATGGTACAACTTCCCAAAACATTATGCGACTTAATGGTTTAAATAATTTCttt                                                                                                                          <  1:699698/141‑1 (MQ=255)
                      agtagCCGTCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTGTCATTAATCGCATCAAAATATGGTACAACTTATCAAAACATTATGAGACTAAATGGCTTAAATAATTTCtttatt                                                                                                                       >  2:268026/1‑141 (MQ=11)
                      agtagCCGTCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCtttatt                                                                                                                       >  1:298010/1‑141 (MQ=255)
                       gtagCCGTCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGCCGCAGAATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTGAATAATTTctttcttt                                                                                                                      >  1:979549/1‑141 (MQ=18)
                         agCCGTCCATCAACGAACTCAGTTTGCGGATAATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTttatttat                                                                                                                    <  1:903801/141‑1 (MQ=255)
                         agCCGTCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTttatttat                                                                                                                    >  2:175218/1‑141 (MQ=255)
                           ccGTCCATCCACGAACTTAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATcc                                                                                                                  <  1:393007/141‑1 (MQ=255)
                            ctgccaccaACGAACTCCGGTGGCGGATCATACTACACAGTACAAGCCGGCGACTCATTTTCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCa                                                                                                                 <  2:138209/138‑1 (MQ=18)
                            cGTCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCa                                                                                                                 <  1:274434/141‑1 (MQ=255)
                            cGTCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTAAAACTTACCAGACCATTAGACTACTTAATGGTTTAAATAATTTCTTTATTTATCCa                                                                                                                 <  1:1093174/141‑1 (MQ=18)
                            cGTCCATCAACGAACTCAGGTGGCAGATCATACTACACAGTACAAGCAGGAGACTCATTATCATTACTCGCATCAAAATATGGTACAACTTACCAAAACATTATGAGACTTAATGGTTTAAATAATTTCTTTATTTATCCa                                                                                                                 >  1:808205/1‑141 (MQ=255)
                             gTCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAg                                                                                                                <  1:619212/141‑1 (MQ=255)
                               ccATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTAAATGGTTTAAATAATTTCTTTATTTATCCAGGt                                                                                                              >  2:165255/1‑141 (MQ=255)
                                cATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAGTAGCATCAAGATATGGGACAACGTACCAAACCATTATGCGACTTAATGGTTTAAATAATTTCTTTATATATCCAGGTc                                                                                                             >  2:650065/1‑141 (MQ=17)
                                 aTCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCa                                                                                                           <  1:955073/141‑1 (MQ=255)
                                  ccAACGATCTCAGGTCGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTAAATGGTTAAAATAATTTCTTTATTTATCCAGGTCaa                                                                                                          <  2:210889/140‑1 (MQ=255)
                                  acAACGGACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCaa                                                                                                          <  2:203332/140‑1 (MQ=255)
                                   cAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGGCTCCTTGTCATTAATTGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCaaa                                                                                                         >  2:148932/1‑141 (MQ=255)
                                   cAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTGTCATTAATTGCATCAAAATATGGTACAACTTACCAAAAAATTATGCGACTAAATGGTTTAAATAATTTATTTATGTATCCAGGTCaaa                                                                                                         >  2:1038246/1‑141 (MQ=14)
                                   cAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTGTCATTAATTGCATCAAAATATGGTACAACTTACCAAAAAATTATGCGACTAAATGGCTTAAATAATTTCTTTATATATCCAGGTCaaa                                                                                                         >  2:548285/1‑141 (MQ=11)
                                   cAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTGTCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGATTAAATGGTTTAAATAATTTCTTTATTTATCCGGGTCaaa                                                                                                         >  2:521699/1‑141 (MQ=255)
                                   cAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAAGTTCTTTATTTATTCGGGTCaaa                                                                                                         >  1:123097/1‑141 (MQ=255)
                                          cgcAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTaaa                                                                                                  <  1:16614/139‑1 (MQ=255)
                                          cTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTaaa                                                                                                  <  2:1093775/141‑1 (MQ=255)
                                           tCAGGTGGCGGATCATACTACACAGTACAAGCAGGCTACTCATTATCATTAATCGAATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTaaaa                                                                                                 >  1:858029/1‑141 (MQ=255)
                                           tCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTaaaa                                                                                                 >  1:76210/1‑141 (MQ=255)
                                           tCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAAGTaaaa                                                                                                 >  1:411587/1‑141 (MQ=255)
                                            cAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAg                                                                                                <  2:325582/141‑1 (MQ=255)
                                            cAGGTGGCCGATCATACTACACAGTACAAGCAGGCGACTATGCTTTATTAAGCCCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATACAGGTCAACAATTAAAAg                                                                                                <  2:1059281/141‑1 (MQ=21)
                                              ggTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGCCTTAATGGTTTAAATAATTTATTTATTTATCCAGGTCAAAAATTAAAAGTa                                                                                              >  1:582756/1‑141 (MQ=255)
                                              ggTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGATTAAATGGTTTAAATAATTTCTTTATTTATCCGGGTCAAAAATTAAAAGTa                                                                                              >  2:644761/1‑141 (MQ=255)
                                              ggTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCGAGGTCAAAAATTAAAAGTa                                                                                              >  1:137308/1‑141 (MQ=255)
                                              ggTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTa                                                                                              >  1:299358/1‑141 (MQ=255)
                                              ggTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTAAATGGTTTAAATAATTTCTTTATATATCCAGGTCAAAAATTAAAAGTa                                                                                              >  2:529740/1‑141 (MQ=255)
                                               gTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTGTAATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGCCTAAATGGCTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTaa                                                                                             >  2:811586/1‑140 (MQ=14)
                                               gTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTAt                                                                                             >  1:159317/1‑141 (MQ=255)
                                               gTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGATCAAAAATTGAAATTaa                                                                                             >  2:42763/1‑140 (MQ=255)
                                                   cGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTGTGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAgg                                                                                         <  1:414857/141‑1 (MQ=255)
                                                   cGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAgg                                                                                         <  1:1102119/141‑1 (MQ=255)
                                                    ggATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTATCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATAAGGt                                                                                        >  2:636751/1‑141 (MQ=255)
                                                     gATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTa                                                                                       >  2:547902/1‑141 (MQ=255)
                                                      aTCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTAc                                                                                      <  2:985696/141‑1 (MQ=255)
                                                         aTACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGc                                                                                   <  1:986299/141‑1 (MQ=255)
                                                         aTACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGc                                                                                   <  1:34800/141‑1 (MQ=255)
                                                                  aGTACCAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAg                                                                          <  1:932274/141‑1 (MQ=255)
                                                                   gTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGt                                                                         <  1:1016299/141‑1 (MQ=255)
                                                                    tACACGCAGGCGACTCATTATCACTAATCGCATCAAAATCTGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTa                                                                        <  1:230127/141‑1 (MQ=255)
                                                                     aCAAGCAGGCTACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTaa                                                                       <  1:277360/141‑1 (MQ=255)
                                                                     aCAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTaa                                                                       <  1:506756/141‑1 (MQ=255)
                                                                     aCAAGCAGGCGACTCATTATCAGTAATGGCATCAAAGTATGGTACAACTTAGCAAAACACTATGCGACTTAATGGTTTAAACAATTTCTTGCTTGATCCAGGTCAAGAATTAAAAGTATCAGGTACTGCTAGCTCAAGTac                                                                       <  2:336606/141‑2 (MQ=255)
                                                                       aaGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACg                                                                     <  2:383988/141‑1 (MQ=255)
                                                                        aGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGCCTAAATGGTTTAAATAATTTCTTTATTTATCCAGGTAACAAATTAAAAGTATCAGGTACTGCTAGCTGAAGTAACGc                                                                    >  2:105013/1‑141 (MQ=21)
                                                                        aGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGc                                                                    >  1:464719/1‑141 (MQ=255)
                                                                            ggCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCg                                                                >  1:150421/1‑141 (MQ=255)
                                                                              cGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGa                                                               <  1:236003/141‑1 (MQ=255)
                                                                               gACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGc                                                             <  1:562946/141‑1 (MQ=255)
                                                                               gACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGc                                                             >  2:888427/1‑141 (MQ=255)
                                                                                ccTCATTATCATTAATCGCATCAAAAGATTGCACAACCTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTACTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAATAACGCTGCGAGCa                                                            <  2:751213/140‑1 (MQ=255)
                                                                                aCTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCa                                                            <  1:837995/141‑1 (MQ=255)
                                                                                    attatcattaATCGCATCACAATATGGTACCACTTACCCACACATTCTGCCAATTAATCGTGTAAGTAATTTCTTTATCTATCCAGGGGAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAAtag                                                        <  2:1062375/141‑1 (MQ=255)
                                                                                    attatcattaATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAAtag                                                        <  2:1108717/141‑1 (MQ=255)
                                                                                          attaATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAACTTAAAAGTATCAGGTACGGCTAGCGCAAGTAACGCTGCGAGCAATAGTAGCCg                                                  <  1:1093859/141‑1 (MQ=255)
                                                                                          attaATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCg                                                  <  1:563581/141‑1 (MQ=255)
                                                                                          attaATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCg                                                  <  1:175218/141‑1 (MQ=255)
                                                                                               tCGCATCAAAATATGGTACAACTTACCAAGACATCATGCGACTTATTGGTTTGAATCGGCGCTTTCTTTCTTCCGCTCAAAAATTGAAAGTATCAGGTACTGCGAGCCCATGCCGCGCTGCGAGCAATAGTAGCCGTCCAt                                             <  2:979549/141‑1 (MQ=255)
                                                                                                cGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCACGTAACGCTGCGAGCAATAGTAGccttgcttc                                            >  2:803987/1‑134 (MQ=255)
                                                                                                cGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCACGTCAAAAATTAAAAGTATCAGGTACTACTAGCTCAACTAACGATGCGAGCAATAGTACACGTCGAta                                            >  2:53832/1‑140 (MQ=25)
                                                                                                 gCATCAAAATATGGTAGAAGTTTCGAAAACATTATGCGACTTAATGGTTTAAATACTTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTACATCa                                           >  1:265546/1‑141 (MQ=255)
                                                                                                    tCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTGATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATCGTAGCCGTCCATCTACg                                        >  2:740970/1‑141 (MQ=255)
                                                                                                     cAAAAGATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGa                                       >  2:927323/1‑141 (MQ=255)
                                                                                                         atatGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTc                                   >  1:97676/1‑141 (MQ=255)
                                                                                                            tgatACAACTGAGCATAACATTACTAAACTTAAGGGTTTAAATAATTTCTTTATTTATCCAAGTTAAACAATACAAGTATCAGATACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAgg                                <  2:274712/138‑1 (MQ=25)
                                                                                                            tGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAgg                                <  1:1018116/141‑1 (MQ=255)
                                                                                                            tGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAATAACGCTGCGAGCAATAGTAGCCGTCCATAAACGAACTCAgg                                >  2:145838/1‑141 (MQ=255)
                                                                                                            tGGTACAACTGACCAAAACATTATGCCACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCCGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAgg                                <  1:233433/141‑1 (MQ=255)
                                                                                                               ttcatcctACCAAAACCACATGCAACTTCATCGTGTAAGTAATTTCTTCCTTTATCCAGGTTAAACCGTAATAGTACTAGATGCTGCTAGCTCAAGTAATGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAGGTgg                             <  2:252455/134‑1 (MQ=12)
                                                                                                                 cgaCTTGCCAAAACATTATGCGACTTAATGGATTAAATAATTTCTCTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAGGTGGCg                           <  2:569440/139‑1 (MQ=255)
                                                                                                                    cTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGAt                        <  1:927323/141‑1 (MQ=255)
                                                                                                                     ttACCAAAACATTATGCGATGTAATGGTTTAAATAATTTCTTTATTTATCCAGGCCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATc                       >  1:716980/1‑141 (MQ=255)
                                                                                                                     ttACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATc                       >  1:628364/1‑141 (MQ=255)
                                                                                                                       aCCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAGGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCAt                     <  1:599241/141‑1 (MQ=255)
                                                                                                                       aCCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCAt                     <  1:796457/141‑1 (MQ=255)
                                                                                                                        ccAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAAATCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCATa                    >  2:721989/1‑141 (MQ=255)
                                                                                                                             ccATTATGCAATTAAATGGTTTAAACAATTTCTCAATTTATCCAGGTCTAAAAATAAAAGTACCAGAAACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCATACtata               <  2:1014980/140‑1 (MQ=25)
                                                                                                                             aCATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCATACtata               >  1:252319/1‑141 (MQ=255)
                                                                                                                              cATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCTAGCAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCATAcgaaac              >  2:924680/1‑136 (MQ=255)
                                                                                                                                ttATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATAATACTATACaa            >  2:78328/1‑140 (MQ=255)
                                                                                                                                ttATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCGCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCATACTATACAg            >  1:1037203/1‑141 (MQ=255)
                                                                                                                                    gCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAGGTGGGGGATCATACTATACAGTACa        >  2:846748/1‑141 (MQ=255)
                                                                                                                                          tAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCAGCAACGAACTCAGGTGGCGGATCATACTATACAGTACAAGCAgg  <  1:650065/141‑1 (MQ=255)
                                                                                                                                          tAATGGTTTAAATAATTTCTTCATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCATACTATACAGTACAAGCAgg  <  2:1081648/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                             
CTCAAGTAATTCTACGAGTAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCATACTACACAGTACAAGCAGGCGACTCATTATCATTAATCGCATCAAAATATGGTACAACTTACCAAAACATTATGCGACTTAATGGTTTAAATAATTTCTTTATTTATCCAGGTCAAAAATTAAAAGTATCAGGTACTGCTAGCTCAAGTAACGCTGCGAGCAATAGTAGCCGTCCATCAACGAACTCAGGTGGCGGATCATACTATACAGTACAAGCAGG  >  USA300TCH1516_ALE/491789‑492067

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: