Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 92,987 A→G 5.2% intergenic (‑72/‑65) ricR_1 ← / → glpE Copper‑sensing transcriptional repressor RicR/Thiosulfate sulfurtransferase GlpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE92,9870AG 5.2% 1536.1 / 43.6 452intergenic (‑72/‑65)ricR_1/glpECopper‑sensing transcriptional repressor RicR/Thiosulfate sulfurtransferase GlpE
Reads supporting (aligned to +/- strand):  ref base A (212/216);  new base G (22/2);  total (234/218)
Fisher's exact test for biased strand distribution p-value = 2.93e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 472 total aligned reads displayed.

TTTTTTTATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGTTTAT  >  USA300TCH1516_ALE/92897‑93077
                                                                                          |                                                                                          
tttttttATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTaga                                                                                            <  2:5185969/91‑1 (MQ=255)
tttttttATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTaga                                                                                            <  1:2025217/91‑1 (MQ=255)
 ttttttATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAATATATATTAATAGATAAACACAAATGTGTCAATTACCCCTAGAg                                                                                           <  2:4075505‑M2/91‑2 (MQ=255)
 ttttttATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAg                                                                                           <  1:5722069‑M2/91‑2 (MQ=255)
     ttATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAggtat                                                                                       >  1:1415162‑M2/1‑86 (MQ=255)
      tATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAggtatt                                                                                      >  1:3546216‑M2/1‑85 (MQ=255)
       aTCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAggtattt                                                                                     <  2:1685527‑M2/91‑8 (MQ=255)
        tCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAggtatttg                                                                                    >  1:4934153‑M2/1‑83 (MQ=255)
        tCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATACATAAACACAAAAGTCTCCAATACACCTAAAtgtatttt                                                                                    <  2:321611‑M2/91‑9 (MQ=255)
         cATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAggtatttga                                                                                   >  2:6287023‑M2/1‑82 (MQ=255)
         cATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAggtatttga                                                                                   <  1:5878416‑M2/91‑10 (MQ=255)
         cATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAggtatttga                                                                                   <  2:6135039‑M2/91‑10 (MQ=255)
           tAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAggtatttgaca                                                                                 >  1:2411208‑M2/1‑80 (MQ=255)
           tAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAggtatttgaca                                                                                 >  2:3700177‑M2/1‑80 (MQ=255)
            taTTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCTCTATATGTATTTGACaa                                                                                >  2:227610/2‑91 (MQ=255)
              ttCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGt                                                                              <  2:1226479/91‑1 (MQ=255)
              ttCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGt                                                                              <  2:5722069/91‑1 (MQ=255)
                cATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtc                                                                            <  1:4125427/91‑1 (MQ=255)
                cATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtc                                                                            >  2:464359/1‑91 (MQ=255)
                 aTAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtcc                                                                           >  1:5109203/1‑91 (MQ=255)
                   aGTCATAGCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccat                                                                         <  1:5575847/91‑1 (MQ=255)
                   aGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccat                                                                         <  2:3599360/91‑1 (MQ=255)
                   aGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccat                                                                         <  1:4560899/91‑1 (MQ=255)
                   aGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccat                                                                         <  2:5473240/91‑1 (MQ=255)
                    gTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccatc                                                                        >  2:20733/1‑91 (MQ=255)
                    gTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccatc                                                                        >  2:4841549/1‑91 (MQ=255)
                     tCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTTTTTGACAAGTtccatcc                                                                       >  2:3935568/1‑91 (MQ=255)
                     tCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccatcc                                                                       <  1:2071475/91‑1 (MQ=255)
                     tCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccatcc                                                                       <  1:4339008/91‑1 (MQ=255)
                     tCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccatcc                                                                       >  2:2843880/1‑91 (MQ=255)
                      cATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccatcca                                                                      <  2:6937532/91‑1 (MQ=255)
                      cATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccatcca                                                                      >  2:6620859/1‑91 (MQ=255)
                       aTACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCaa                                                                     <  2:6824941/91‑1 (MQ=255)
                       aTACCTCCACTTGTAATTGAATAAAAAGATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCaa                                                                     <  2:6000203/91‑1 (MQ=255)
                        tACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAAc                                                                    <  1:371263/91‑1 (MQ=255)
                         aCCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACt                                                                   <  2:1170313/91‑1 (MQ=255)
                          ccTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTg                                                                  >  1:936809/1‑91 (MQ=255)
                            tCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGtt                                                                <  2:4929986/91‑1 (MQ=255)
                            tCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGtt                                                                >  2:745139/1‑91 (MQ=255)
                             ccACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGttt                                                               <  2:6213273/91‑1 (MQ=255)
                             ccACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGttt                                                               <  1:4828518/91‑1 (MQ=255)
                             ccACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGttt                                                               >  1:342907/1‑91 (MQ=255)
                             ccACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGttt                                                               <  2:1591479/91‑1 (MQ=255)
                              cACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTa                                                              >  1:22402/1‑91 (MQ=255)
                              cACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTa                                                              >  2:3741329/1‑91 (MQ=255)
                              cACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTa                                                              >  2:1988749/1‑91 (MQ=255)
                                cTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTaaa                                                            >  2:4451899/1‑91 (MQ=255)
                                  tttAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAAt                                                          >  1:3120192/1‑91 (MQ=255)
                                    tAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGtt                                                                >  1:5336979/1‑83 (MQ=255)
                                        tGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCt                                                    >  1:4196375/1‑91 (MQ=255)
                                        tGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCt                                                    <  1:2621582/91‑1 (MQ=255)
                                         gAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTa                                                   >  2:846900/1‑91 (MQ=255)
                                          aaTAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTAc                                                  <  2:2936117/91‑1 (MQ=255)
                                             aaaaaTATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAgg                                               <  1:1861704/91‑1 (MQ=255)
                                              aaaaTATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGt                                              >  2:109632/1‑91 (MQ=255)
                                               aaaTATATTAATAGATAAACACAAATGTGTCAAATACCTCTAGGggtgtttgacacatttgtgtttatctattaatatatttttattcaat                                             >  1:427616‑M2/1‑44 (MQ=255)
                                               aaaTATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTa                                             >  1:5925889/1‑91 (MQ=255)
                                                aaTATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTAt                                            >  1:6374391/1‑91 (MQ=255)
                                                 atatatTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTAtt                                           <  1:1011334/91‑1 (MQ=255)
                                                 atatatTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTAtt                                           <  2:5922430/91‑1 (MQ=255)
                                                  tatatTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTAttt                                          >  2:2630479/1‑91 (MQ=255)
                                                   atatTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTAtttt                                         >  1:613582/1‑91 (MQ=255)
                                                   atatTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTAtttt                                         >  1:4630074/1‑91 (MQ=255)
                                                    tatTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTAttttt                                        >  1:5090239/1‑91 (MQ=255)
                                                       tAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAgg                                     >  1:3606813/1‑91 (MQ=255)
                                                        aaTAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAggg                                    <  1:2676035/91‑1 (MQ=255)
                                                         aTAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGa                                   <  1:1829296/91‑1 (MQ=255)
                                                           aGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAgg                                 <  2:3445023/91‑1 (MQ=255)
                                                             aTAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGtt                               <  1:114260/91‑1 (MQ=255)
                                                             aTAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGtt                               <  1:2565030/91‑1 (MQ=255)
                                                             aTAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGtt                               <  1:1580048/91‑1 (MQ=255)
                                                               aaaCACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGttat                             >  1:2889727/1‑91 (MQ=255)
                                                                 acacAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGttatta                           >  2:4598153/1‑91 (MQ=255)
                                                                   acaAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATg                         <  1:3832203/91‑1 (MQ=255)
                                                                    caAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGa                        <  2:5009828/91‑1 (MQ=255)
                                                                    caAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGa                        >  1:1970656/1‑91 (MQ=255)
                                                                    caAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGa                        >  1:6943624/1‑91 (MQ=255)
                                                                      aaTGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGaaa                      <  2:6355120/91‑1 (MQ=255)
                                                                      aaTGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGaaa                      >  1:2332627/1‑91 (MQ=255)
                                                                      aaTGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGaaa                      >  1:81071/1‑91 (MQ=255)
                                                                        tgtgTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTCGGGAGGTTATTATGAAACa                    >  2:2432767/1‑91 (MQ=255)
                                                                            tCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATAc                >  2:2548906/1‑91 (MQ=255)
                                                                              aaaTACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACgg              >  1:6834568/1‑91 (MQ=255)
                                                                              aaaTACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACgg              >  1:3763713/1‑91 (MQ=255)
                                                                              aaaTACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACgg              >  2:6525754/1‑91 (MQ=255)
                                                                                aTACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGgag            <  2:4403754/91‑1 (MQ=255)
                                                                                 tcccccTAGGGGTATTTGACAAGTTCCATCCGACTGTTTAAACTACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGgaga           <  2:5384359/89‑1 (MQ=255)
                                                                                 tACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGgaga           <  1:1517933/91‑1 (MQ=255)
                                                                                  acctcTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGaa          >  2:1390331/5‑91 (MQ=255)
                                                                                   ccccTAGAGGTATTTGACAAGTTCTATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGaaa         <  1:6624432/91‑1 (MQ=255)
                                                                                    ctcTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGaaaa        >  1:2496694/3‑91 (MQ=255)
                                                                                    cccTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGaaaa        >  2:1512548/1‑91 (MQ=255)
                                                                                     ccTAGAggtatttgacacatttgtgtttatctattaatatatttttattcaattaaaagtggaggtatgactatgaattatgataaaaaaa       <  1:1344550‑M2/91‑86 (MQ=255)
                                                                                     ccTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAg       >  1:404322/1‑91 (MQ=255)
                                                                                      cTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGt      >  1:3376369/1‑91 (MQ=255)
                                                                                      cTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGt      <  2:4673214/91‑1 (MQ=255)
                                                                                       tAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGtt     >  1:3014717/1‑91 (MQ=255)
                                                                                         gagGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGTTTa   <  2:2345871/91‑1 (MQ=255)
                                                                                          agGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGTTTAt  <  1:6504915/91‑1 (MQ=255)
                                                                                          agGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGTTTAt  <  2:431933/91‑1 (MQ=255)
                                                                                          |                                                                                          
TTTTTTTATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGTTTAT  >  USA300TCH1516_ALE/92897‑93077

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: