Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 517,732 T→C 29.7% intergenic (+5800/‑963) recR → / → speA_2 Recombination protein RecR/Arginine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE517,7320TC29.7% 60.5 / 66.4 88intergenic (+5800/‑963)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (33/28);  new base C (8/18);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 6.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

CTGACGAATACTAATCGATCGAAGACTTAATCAAAATAAATGTTTTGCGAAGCAAAATCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTGAGACCGCAAGGTCTCTTTTTTTATGTCTAAA  >  USA300TCH1516_ALE/517593‑517871
                                                                                                                                            |                                                                                                                                           
cTGACGAATACTAATCGATCGAAGACTTAATCAAAATAAATGTTTTGCGAAGCAAAATCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTc                                                                                                                                            <  2:77677/141‑1 (MQ=31)
    cGAATCCTACTCGACCGAAGACTTAGTCAAAATAAATGTTTAGCGAAGCAAAAACAGTATTACTAACAATAGAGTTTTGAATGTATAAATTACATGCATATGTCGGGTCGTTACATCAAAGAGGTCACACCAGTTTCCATg                                                                                                                                        >  2:311668/1‑141 (MQ=2)
         aCTAATCGATCGAAGACTTAATCAAAATAAATGTTTTGCGAAGCAAAATCACCTTAACTTACTGTCTAGTTTTGAATGTCTAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGaa                                                                                                                                   <  1:368850/141‑1 (MQ=255)
          cTAATCGATCGAAGACTTAATCAAAATAAATGTTTTGCGAAGCAAAATCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGCCTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAac                                                                                                                                  >  1:261525/1‑141 (MQ=17)
              tcgatcgaAGACTTAATCAAAATGAATGGTTTGCGAAGGAAAATCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATGCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAg                                                                                                                              <  2:214863/141‑1 (MQ=255)
              tcgatcgaAGACTTAATCAAAATAAATGTTTTGCGAAGCAAAATCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGCGGTCACACCTGTTTCCATGCTGAAcacaa                                                                                                                              >  2:533636/1‑140 (MQ=17)
              tcgatcgaAGACTTAATCAAAATAAATGTTTTGCGAAGCAAAATCACTTTTACTTAATAGCTATTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTAACACCTGTTTCCATGCCGAACACAg                                                                                                                              >  2:318164/1‑141 (MQ=255)
                    gaAGACTTAATCAAAATAAATGTTTTGCGAAGCAAAATCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTa                                                                                                                        >  2:324486/1‑141 (MQ=32)
                    gaAGACTTAATCAAAATAAATGTTTTGCGAAGCAAAATCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTAACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACTCAGAAGTTa                                                                                                                        >  2:419482/1‑141 (MQ=17)
                      ataCTTAATCAAAATAAATGTTTGGCGAAGTAAAATCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAg                                                                                                                      <  2:101796/139‑1 (MQ=12)
                         cTTAATCAAAATAAATGTTTTGCGAAGCAAAATCACTTTTACTTACTATCTAATTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGACGAAGACAGAAGTTACGCtt                                                                                                                   >  2:230580/1‑140 (MQ=255)
                           tAATCAACATAAATGTTTTGCGAAGCAAAATCACTTTTACTTACTATCTAGGTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCt                                                                                                                 <  1:231062/141‑1 (MQ=17)
                           tAATCAAAATAAATGTTTTGCGAAGCAAAATGACTTGTGCCTCAGACTTTGGTTTGAATGTATAAAGTACATTCACATGTCTGGTGACTATAGCAAGGAGGGCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCt                                                                                                                 <  2:453553/141‑1 (MQ=11)
                            aaTCAAAATAAATGTTTTGCGAAGCAAAATCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAATTAGGCTCAAGCCGCCTTCTTGCTGAACACAGAAGTTAAGCTCCtt                                                                                                                >  1:88518/1‑141 (MQ=11)
                              tCAAAATAAATGTTTTGCGAAGCAAACTGACCCTCAGCTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAg                                                                                                              >  2:523561/1‑141 (MQ=11)
                              tCAAAATAAATGTTTTGCGAAGCAAAATCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAg                                                                                                              >  1:360751/1‑141 (MQ=32)
                                aaaaTAAATGTTTTGCGAAGCAAAATCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCg                                                                                                            >  2:58150/1‑141 (MQ=32)
                                      aatGTTTTGCGAAGCAAAATCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCCTAGCGTCCATg                                                                                                      >  2:296153/1‑141 (MQ=17)
                                      aatGTTTTGCGAAGCAAAATCACTTTTACTGACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATg                                                                                                      <  1:199313/141‑1 (MQ=17)
                                         gTTTTGCGAAGCAAAATCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTa                                                                                                   >  1:258130/1‑141 (MQ=32)
                                              gCGAAGCAAAATCAGCTTGCGTTACTAGCTCGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGa                                                                                              <  2:194898/141‑1 (MQ=11)
                                                 aaGCAAAATCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACt                                                                                           <  1:5602/141‑1 (MQ=32)
                                                  aGCAAAATCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACtt                                                                                          >  1:582134/1‑141 (MQ=32)
                                                         tCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTGGTCGAACTTACGTTcc                                                                                   <  2:509361/141‑1 (MQ=17)
                                                           aCTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGc                                                                                 >  2:294226/1‑141 (MQ=32)
                                                           aCTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGc                                                                                 <  1:636491/141‑1 (MQ=32)
                                                                   ttactATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAg                                                                         <  2:257389/141‑1 (MQ=32)
                                                                      gtATCCAGTTTTGAATGTATACACGACATTCATATGGCTGGTGACTATAGCCAGGAGGAGTCACCTGTTTCCATGCCGAACACAGATGTTAAGCTCCTTAGCGCCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAAc                                                                      <  1:56629/140‑1 (MQ=11)
                                                                      ctATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAAc                                                                      <  2:14688/141‑1 (MQ=32)
                                                                      ctATCTAGTTTTCGATGTATAAAGTACATTCATATGGCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAAc                                                                      <  2:255173/141‑1 (MQ=255)
                                                                         tCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTCACGTACCGCCAGAGTAGAACGtt                                                                   >  2:415901/1‑141 (MQ=12)
                                                                                  gAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCaa                                                          >  1:281858/1‑141 (MQ=32)
                                                                                   aaTGTATAAATTACAGTCATATGTCGGGTGACTATAGCAAGGAGGTCACACCTGTTGCCATGCCGAACAAAGACGTAAAGTTTCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgt                                                         >  2:196383‑M2/1‑139 (MQ=255)
                                                                                    aTGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCaaaa                                                        <  1:523561/141‑1 (MQ=32)
                                                                                      gTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCaaaaaa                                                      <  1:620451/141‑1 (MQ=32)
                                                                                       tataAATTACATTCAGATGTCTGGTGACTATAGCAAGGAGGTCGCACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAAt                                                     <  2:605778/141‑1 (MQ=17)
                                                                                          aaaTTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGa                                                  <  1:531430/141‑1 (MQ=32)
                                                                                          aaaTTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGa                                                  <  1:58150/141‑1 (MQ=32)
                                                                                          aaaTTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCGGGCAgttttttaa                                                  <  2:385822‑M2/141‑10 (MQ=255)
                                                                                           aaTTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaat                                                 >  1:560040‑M2/1‑131 (MQ=255)
                                                                                           aaTTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaat                                                 >  1:195335‑M2/1‑131 (MQ=255)
                                                                                            actACAGTCATATGTCTGGTGACTATAGCAGGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGCTGGTAGTCGAACTTACGTTCCGGTAGAAGAGAACGTTGCCAGGCAAAAAATGGATg                                                <  2:232458/139‑1 (MQ=11)
                                                                                            aTTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatc                                                >  1:571091‑M2/1‑130 (MQ=255)
                                                                                            aTTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatc                                                <  2:340041‑M2/141‑12 (MQ=255)
                                                                                                cATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAATTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCa                                                           <  1:70112/126‑1 (MQ=18)
                                                                                                cATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAATTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCa                                                           >  2:70112/1‑126 (MQ=18)
                                                                                                cATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatcaaat                                            <  1:614036‑M2/141‑16 (MQ=255)
                                                                                                cATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatcaaat                                            <  1:406378‑M2/141‑16 (MQ=255)
                                                                                                  ttCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGACCGTTGCCAGGCAgttttttaatcaacttt                                          <  1:253855‑M2/141‑18 (MQ=255)
                                                                                                    cATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTACGCTCCTTAGCGTCGATGGTAGTCCAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatcaaattttg                                        <  2:107183‑M2/141‑20 (MQ=255)
                                                                                                      tatGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatcaaattttggt                                      <  1:417465‑M2/141‑22 (MQ=255)
                                                                                                           cTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCCCAAAATAGAGTCGATGAACCgacttg                                 >  2:384602/1‑136 (MQ=11)
                                                                                                           cTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTg                                 >  1:241028/1‑141 (MQ=32)
                                                                                                           cTGGTGACTATAGCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTCGAGTAGAACGTTGCCAGGCAgttttttaatcaaattttggttaaaa                                 <  2:579600‑M2/141‑27 (MQ=255)
                                                                                                           cTGGTGACTATAGCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatcaaattttggttaaaa                                 >  1:590777‑M2/1‑115 (MQ=255)
                                                                                                             ggTGACTATAGCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTCGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatcaaattttggttaaaaaa                               >  1:457934‑M2/1‑113 (MQ=255)
                                                                                                              gTGACTATATCAAGGAGGTCACATCTGTTTCTATGCAGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCATGCAAAAAATGGATGCGATGAGCCGCATTgaga                              <  2:333257/141‑1 (MQ=11)
                                                                                                              gTGACTATAGCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCATGCAgttttttaatcaaattttggttaaaaaat                              >  1:617027‑M2/1‑112 (MQ=255)
                                                                                                                gACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTGAGAcc                            >  1:140624/1‑141 (MQ=32)
                                                                                                                  cTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTCAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGAATGCGATGAGCCCCATTCAGActcc                          >  2:413749/1‑138 (MQ=11)
                                                                                                                  cTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACAAAGAAGTTAAGCTCCTTAGCGCGGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAAGGGATGCGATGAGCCGCATTGAGACCGc                          >  1:622966/1‑141 (MQ=255)
                                                                                                                  cTATAGCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatcaaattttggttaaaaaataaaa                          >  1:205386‑M2/1‑108 (MQ=255)
                                                                                                                     taGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTGAGACCGCAAg                       >  1:33802/1‑141 (MQ=32)
                                                                                                                     taGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACATTGCCAGGCAGAAACGGGAAGCGATGAGCCGCATCTAGACCGCAAg                       >  2:326912/1‑141 (MQ=11)
                                                                                                                     taGCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatcaaattttggttaaaaaataaaatgg                       <  1:638392‑M2/141‑37 (MQ=255)
                                                                                                                      aGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCTTCTATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATTAGCCGCATTGAGACCGCAAgg                      <  1:404079/141‑1 (MQ=12)
                                                                                                                      aGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTGAGACCGCAAgg                      <  1:96843/141‑1 (MQ=32)
                                                                                                                       gCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTCGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgctttttaattaacttttggttaaaaaataaaatggac                     <  2:567068‑M2/141‑39 (MQ=255)
                                                                                                                       gCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGCACGTTGCCAGGCAgttttttaatcacattttggcaaaaaaatataattgaa                     <  2:475811‑M2/141‑39 (MQ=255)
                                                                                                                       gCAAGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatcaaattttggttaaaaaataaattggac                     <  1:549665‑M2/141‑39 (MQ=255)
                                                                                                                         aaGGAGGTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgtcttttaatcaaattttggttaaaaaataaaatggacaa                   <  2:245886‑M2/141‑41 (MQ=255)
                                                                                                                          aggaggTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTGAGACCGCAAGGtctc                  <  1:324486/141‑1 (MQ=32)
                                                                                                                          aggaggTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTGAGACCGCAAGGtctc                  >  1:488441/1‑141 (MQ=32)
                                                                                                                           ggaggTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatcaaattttggttaaaaaataaaatggacaaga                 >  1:332255‑M2/1‑99 (MQ=255)
                                                                                                                            gaggTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGAGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTGAGACCGCAAGGTCTCtt                <  2:289488/141‑1 (MQ=17)
                                                                                                                            gaggTCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatcaaaatttggttaaaaaataaaatggacaagat                <  1:233486‑M2/141‑44 (MQ=255)
                                                                                                                             aggTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCCTTGAGACCGCAAGGtacattt               >  2:507516/1‑135 (MQ=12)
                                                                                                                             aggTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTGAGACCGCAAGGTCTCttt               >  2:478533/1‑141 (MQ=32)
                                                                                                                             aggTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTGAGACCGCAAGGTCTCttt               >  2:230537/1‑141 (MQ=32)
                                                                                                                              ggTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTCGCGTCCCAGGACCCACCCCCCAGGGGCCGCTAGAGTATAACGTTGCCCGGCAAAAAATGCATGCGATGAGCCGCATTGAGACCGCAAGGTCTCtttt              >  1:304620/1‑141 (MQ=2)
                                                                                                                                 tCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGCACTTACGTTCCGCTAGAGTAGAACGTTGCAAGGCACAAAATGGATGCGATGAGCCGCATTGAACCCGCAAGGTCTAttgatt            >  2:570495/1‑137 (MQ=11)
                                                                                                                                 tCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTGAGACCGCAAGGTCTCtttttt            >  2:55520/1‑141 (MQ=32)
                                                                                                                                 tCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATGAGCAGCATTAAGACCGCAAGGTCTCtttttt            >  2:306910/1‑141 (MQ=17)
                                                                                                                                 tCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatcaaattttggttcaaaaataaaatggacaagataaaa            <  1:18343‑M2/141‑48 (MQ=255)
                                                                                                                                 tCACACCTGTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACGTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatcaaattttggttaaaaaataaaatggacaagataaaa            <  2:485078‑M2/141‑48 (MQ=255)
                                                                                                                                       cTGTTTCCATGCCAAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATCAGCCGCACTGAGACCTCAATGTCTCTTTTTTTATGTc      <  1:504992/141‑1 (MQ=11)
                                                                                                                                         gTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGCACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatcaaattttggttaaaaaataaaatggacaagataaaaaaagttat    <  2:247762‑M2/141‑56 (MQ=255)
                                                                                                                                         gTTCCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAgttttttaatcaaattttggttaaaaaataaaatggacaagataaaaaaagttat    <  1:526673‑M2/141‑56 (MQ=255)
                                                                                                                                          tttCCTGGCGGACCACAGTAGTGCAGCTGCTTAGAGTCGATGGTAGTCGAACTTACGGTCCGCTCGATGAGAACGTTGCCAGGCAAATAATGGATGCGATGAGCGGCATTGAGACCGCAAGGTCTCTTTTTTTATGTCTaa   <  1:186130/141‑1 (MQ=11)
                                                                                                                                          tttCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTGAGACCGCAAGGTCTCTTTTTTTATGTCTaa   <  1:502682/141‑1 (MQ=35)
                                                                                                                                           ttCAATGCCGAACACTGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGCATGCGATGAGCCGCATTGAGACCGCAAGGTCTCTTTTTTTATGTCTaaa  <  2:106141/141‑1 (MQ=12)
                                                                                                                                            |                                                                                                                                           
CTGACGAATACTAATCGATCGAAGACTTAATCAAAATAAATGTTTTGCGAAGCAAAATCACTTTTACTTACTATCTAGTTTTGAATGTATAAATTACATTCATATGTCTGGTGACTATAGCAAGGAGGTCACACCTGTTTCCATGCCGAACACAGAAGTTAAGCTCCTTAGCGTCGATGGTAGTCGAACTTACGTTCCGCTAGAGTAGAACGTTGCCAGGCAAAAAATGGATGCGATGAGCCGCATTGAGACCGCAAGGTCTCTTTTTTTATGTCTAAA  >  USA300TCH1516_ALE/517593‑517871

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: